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Figure Captions 


Figure 

1. 

Oxygen and ozone cross sections with OAO filters 
transmission functions shown. 

Figure 

2. 

Schematic view of occultation. 

Figure 

3. 

Typical occultation row data. 

Figure 

4. 

Geometry of oc cut at ion for any stellar object. 



Figure Captions Appendix IV. 

Figure 

1, 

Tangent ray 'height error as a function of satellite 
height error and altitude. 

F igure 

2. 

Tangent ray height error as a function of satellite 
latitude error. 

F igure 

3. 

Tangent point longitude error as a function of 
satellite latitude error. 

Figure 

4. 

Satellite tangent ray height error as a function 


of satellite longitude error. 


1. Introduction 


This report describes the results of a study of the stellar 
occultation technique for measuring the composition of the atmosphere. 
The work reported here was supported by NASA Contract NAS 1 9958. 

The stellar occultation technique is based upon classical 
absorption spectroscopy, utilizing spacecraft -borne photometers as the 
detector and a star as the light source. During occultation the ultra- 
violet stellar spectrum of the star shows absorption features which are 

uniquely related to various atmospheric constituents. In particular, 

o o 

the spectral region between 1200A and 2000A is related to absorption 

o o 

by molecular oxygen and the region between 2000A and 3500A is related 
primarily to absorption by ozone. This is illustrated in Figure 1 where 
the photo absorption cross sections of both molecular oxygen and ozone 
are shown. 

The intensity of starlight was monitored during the occultation 
using the Wisconsin stellar ultraviolet photometers aboard the Orbiting 
Astonomical Observatory (OAO-A2). A schematic diagram of an occul- 
tation is shown in Figure 2 where the change in intensity at a given wave- 
length is illustrated. The vertical projection of the attenuation region is 
typically 60 km. deep for molecular oxygen and 30 km. deep for ozone. 

Intensity profiles obtained during various occultations were 
analyzed by first determining the tangential column density of the absorb- 
ing gases, and then Abel inverting the column densities to obtain the 
number density profile. Errors are associated with each step in the 
inversion scheme and have been considered as an integral part of this 
study. 
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CROSS SECTION, cm 



Figure 1. Oxygen and ozone cross sections with GAO filters 
transmission functions shown. 




Most of the work carried out under this contract has been 
published in the open literature, thus, we include these papers in this 
report where appropriate. A general description of this study and the 
results were published by Hays et al. (1972) and is included as Appendix 
I of this report. 

2 * Data Procurement 

The data used in this study were obtained through the courtesy 
of the O. A. O. Wisconsin Experimenters and N. A. S. A. Goddard Space 
Flight Center. The principle difficulty associated with the data procure- 
ment centered on obtaining intensity information at high rates. It was 
possible to sample the stellar photometer accumulator at high frequency, 
but only bits 7 through i4 of this register were available in this mode of 
operation. During this procedure the photometers continued to accumu- 
late for 1.28 periods of U.54U80 seconds. Upon completion of this process 
the register was locked until a new command was received, at which time 
the cycle was repeated. Unfortunately, a finite period of time was 
required for the new command to be sent and a gap appears in the data 
every i28 samples. This caused many occupations to be rejected due to 
re -set blanks appearing in the middle of the attenuation period. 

A typical data record is shown in Figure 3, we note that the 

bottom of each record corresponds to b4 accumulated counts and the top 

to lb, 384 accumulated counts. The overflow within this sub-set is 

obvious. The data analyses utilized differences between individual 

samples and a smooth interpolation to determine the overflow and count 

accumulated during a sample period. The record labeled Channel 44 in 

o 

Figure 3 is from the 1500A filter used to measure molecular oxygen. 
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The stellar photometers have broad band filters which 
o 

typically cover a 2U0A wide region in the spectrum. The relative 
transmission of these filters are shown at the bottom of Figure 1. 

Due to the broad spectral band of these filters it was necessary to 
measure the spectrum of all stars used in the study. This was done 
prior to the occultation using the O. A. O. spectrometer. 

In addition to intensity and spectral information it was 
necessary to determine the altitude of the tangent ray. This requires 
an accurate knowledge of the spacecraft positions as well as the 
coordinates of the star. We acknowledge the considerable effort of the 
NASA Goddard Space Flight Center who carefully re -checked the satellite 
ephemeris in order to achieve the accuracy required for our study. 

3. Tangent Ray Geometry 

Throughout an occultation period, the tangent ray point follows 
a curved path above the surface of the earth. The coordinates used to 
describe this path are shown in Figure 4. In this system, the x-axis points 
toward the First Point of Aries, the y~axis is also in the equatorial plane 
pointing 90° east of Aries, and the z-axis points toward the North Pole. 

The time of scan enters into the coordinate system through the Greenwich 
hour angle of the x-axis, A . 

o 

To obtain the height of the tangent point, the angle ^5 at the 
spacecraft between rays to the star and to the center of the earth must be 
evaluated in terms of the direction numbers of the star, 

1 = cos cos 5 
m= sin<< cos £ 
n = sin £ 
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and the direction cosines of the satellite position, 

a = cos X cos jf cos A - sin X cos p sin a 
s s g s • s g 

b = s;n X cos i cos > + cos i cos 6 sin A 
s s g s s g 

c - sin w 
r s 


where ^ ascension 

(V = declination 
X = longitude 
0 ■= latitude 

Subscripts s refers to the spacecraft and subscript o will 
refer to the tangent ray point. Then utilizing the above quantities. 


where 


r = r sinv. ; 
o s v 



(1) 

6 * IT - COS 

* [ al+bm+cn j 

9 


0 = tan * 

f z <x 2 + y 2 >“*] 
t O O J o J 

9 

(2) 

X 0 = cos ’ 1 

f x (x 2 + y 2 )-*] 
Loo J o •* 

Xg. 

(3) 

If > 

o, o i X Q 1 180; 

i f y 0 * o. 

180 £ A ^ 360. 
o 


The additional formulae needed to evaluate Equations l -3 are 
x q = r g (a+1 cos $ ) 
y D = (b+m cos {J ) 


Z o = r s ^ c+n cos ? ) 


Stellar refraction is neglected in our analysis due to the low 
index of refraction at the altitudes considered in this study. 


4. Data Processing and Error Analysis 

The technique of data processing and the analysis of resulting 
errors has been reported during this study by Roble and Hays (t972). This 
paper is included as Appendix II and is a complete description of the 
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results of this study. A complete description of the computer routine 
used to evaluate the occultation results obtained in this study is included 
as Appendix III of this report. The reader is referred to these detailed 
appendices for further information on the results of our study. A 
separate discussion of influence of orbital errors is included in Appendix 
IV. 

5. Data Analysis 

A very large number of scans have been reduced in the course 
of this contract period. Table 1 summarizes these results and includes 
all orbits which have yielded at least a partial result e. g. , a molecular 
oxygen profile or an ozone profile. To provide information on quality in 
the individual channels, the following code is used 
c = processing completed 

+ = resulting number density profile, n{r Q ), acceptable 
- = n(r Q ) not acceptable 
? = n(r Q ) questionable 

Other information in the table refers to the starting time 
(GMT and local) of the scan, the location of the ray tangent point, the 
target star, and the nature of some of the difficulties. Starting times 
are only approximate, as they vary up to 2 u seconds from channel to 
channel and they must be entered to the nearest half -second for accept- 
able accuracy in height of ray tangent point. 

These results have been analyzed from the point of view of 
their geophysical information and are to be published in journals during 
the next few months. All papers are in press at this time. Detailed 
computer output containing tabular input data and results are available to 
the contracting agency on request, but are too lengthly to include in this 
report. 9 



Table I 


GMT local 


orbit 

date 

_h_r 

min 

_hr_ 

min 

Long, (deg) Lat. (deg) 

5778 

1/13/70 

5 

28 

23 

52 

-84 

48 

5779 

1/13/70 

7 

8 

23 

51 

-109 

48 

8884 

8/17/70 

4 

55 

5 

55 

15 

4 

8898 

8/18/70 

4 

19 

5 

51 

23 

9 

8899 

8/18/70 

5 

59 

5 

51 

-2 

9 

8913 

8/19/70 

5 

22 

5 

54 

8 

11 

8942 

8/21/70 

5 

55 

4 

27 

-22 

24 

8943 

8/21/70 

7 

35 

4 

23 

-48 

24 

8956 

8/22/70 

5 

18 

4 

26 

-14 

24 

8957 

8/22/70 

6 

58 

4 

22 

-39 

25 

8971 

8/23/70 

6 

20 

4 

44 

-24 

27 

8985 

8/24/70 

5 

45 

4 

13 

-23 

24 

8986 

8/24/70 

7 

25 

4 

13 

-48 

24 

9000 

8/25/70 

6 

47 

5 

17 

-32 

27 

9014 

8/26/70 

6 

10 

4 

30 

-25 

26 

10795 

12/28/70 

2 

31 

21 

55 

-69 

42 

10796 

12/28/70 

4 

12 

22 

12 

-91 

43 

10902 

1/4/71 

13 

29 

22 

38 

-223 

-16 

11795 

3/7/71 

15 

59 

0 

33 

12 8 

-34 

12178 

4/3/71 

6 

8 

22 

33 

-114 

36 

14579 

9/17/71 

5 

42 

2 

51 

-43 

26 

14580 

9/17/71 

6 

22 

1 

50 

-68 

26 


ozone 
o o 


star 

2390A 

2460A 

o? 

remarks on data loss 

r Peg 

- 


+ 

O 3 filter 

If Peg 

- 


+ 

Og filter 

Lup 



+ 

O 3 reset 

p Lup 

9 

? 

? 

Height 

p Lup 

9 

? 


Height 0 ^ enevenness 

Lup 



+ 

O 3 reset 

'T Sco 

+ 

- 

? 

O 3 data cutoff, O 2 reset 

T Sco 



+ , 


T Sco 

+ 



Og reset 

T Sco 

+ 

- 

+ 

Height 

$ Sco 

+ 

+ 

- 

O 2 reset 

T Sco 


+ 

+ 


T* Sco 


+ 

+ 


S Sco 

+ 

+ 

+- 


l Sco 

+ 

+ 



i Her 



+ 

€ Aur spectrum used 

L Her 

- 

- 

- 

£ Aur spectrum used 





Height 

And 

+ 


+ “ 


$ Ori 

+ 


+ 


i Ori 


+ 

+ 


O Oph 

+ 


- 

O 2 reset 

OQph 

? 


+ 

O 3 values high in bulge 


In addition, orbits 13161, 13187, 13189, and 13195 of June 10 and June 12’ 1971, provided high quality data 


which are awaiting accurate satellite position reports. 



6. Geophysical Results 

Interpretations of the results obtained under this contract 
have been presented in papers published in the open literature. Molecular 
oxygen results were compared with theory and previous measurements 
by Hays and Roble (1972b). Ozone results were discussed by Hays and 
Eoble (i972c), Roble and Hays (1972b) and Roble and Hays (1972c). These 
papers are included as Appendices V through VIII. 

7. Summary of Study 

The present study of the stellar occultation technique of 
recovering atmospheric composition has led to the following conclusions: 

a) Any constituent which causes significant absorption in a stellar 
spectrum can be monitored by using the ultraviolet occultation 
technique. 

b) The error associated with the inversion of occultation data limits 

< 

the altitude region for which results of high quality can be obtained. 
In most cases the high quality data results from regions where the 
intensity lies between 1U% to 90% of the unattenuated flux. Finite 
difference errors are responsible for limiting the structural detail 
which can be recovered. Orbital errors are serious and are 
reflected in the altitude scale of the recovered density profiles. 

c) Serious geophysical studies can be carried out using as principal 
data the density profiles of molecular oxygen and ozone obtained 
from stellar occultation. The following major geophysical features 
have been identified during this study: 

(1) Molecular oxygen shows a strong solar cycle variation in the 
lower thermosphere. This probably results from changep in 
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the constant pressure level surfaces at the base of the 
thermosphere. 

(2) Molecular oxygen in equatorial regions shows strong 
magnetic storm variations due to large scale dynamic 

motions resulting from polar region heating. > 

(3) Ozone in the mesopause region is controlled by wet 
chemistry and supports modern photochemical theories. 
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APPENDIX I. 


Terrestrial Atmospheric Composition from Stellar Occultations 
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Reprinted from 
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Terrestrial Atmospheric Composition from Stellar 
Occultations 

Abstract. Stellar ultraviolet light transmitted through the earth’s upper atmo- 
sphere is strongly absorbed by ozone and molecular oxygen. The stellar ultraviolet 
photometers aboard the Orbiting Astronomical Observatory ( OAO-2 ) satellite 
have measured the intensity changes of several stars during occultation of the star 
by the earth's atmosphere. From the occultation data the nighttime vertical num- 
ber density profiles of molecular oxygen at altitudes from 120 to 200 kilometers 
and of ozone at altitudes from 60 to 100 kilometers have been obtained. 


In the earth's upper atmosphere, 
stellar ultraviolet light is strongly ab- 
sorbed in the Schumann-Runge contin- 
uum of Oo and the Hartley continuum 
of O s . By monitoring the intensity of 
ultraviolet starlight in these continuum 
regions from a satellite as the star is 
occulted by the earth’s atmosphere; we 
are able to obtain information on the 
number density profile of 0 2 in the 
lower thermosphere and 0 ;< in the up- 
per mesophere (/, 2). In this report 
we describe the technique used to ob- 
tain the number density profiles from 
stellar occultation measurements. The 
data were obtained by the University 
of Wisconsin stellar photometers aboard 
the Orbiting Astronomical Observatory 
(OAO-2) satellite. We discuss here the 
inversion process and also show the 0 2 
and O a distributions determined from 
a typical occultation scan. 

In the occultation technique classical 



Fig. I . Geometry of ultraviolet stellar 
occultation. 


absorption spectroscopy is used to de- 
termine the number density profile of 
the absorbing species in the upper 
atmosphere. The star is the source of 
ultraviolet light, the OAO stellar pho- 
tometers are the detectors, and the 
atmosphere 'between them is the ab- 
sorption cell. During the occultation 
process, the ultraviolet light is selective- 
ly absorbed in spectral regions for 
which 0 2 and O a have large absorption 
cross sections. The intensity of the 
transmitted ultraviolet light is related to 
the number density along a tangential 
column of the absorbing species (N<) 
by Beer’s law 

I(\,r o) = 

MM exp ■ AMr 0 )J (i) 

where l x (\) is the unattenuated inten- 
sity above the atmosphere of the star at 
wavelength A, ^(A) is the absorption 
cross section of the tth absorbing spe- 
cies, and N,(r 0 ) is the tangential column 
number density of the ith absorbing 
species at a tangent ray height r 0 . The 
star’s spectrum is measured above the 
atmosphere by the OAO ultraviolet 
spectrometer. By also knowing the ab- 
sorption cross section, one can relate 
Ni(r 0 ) to the intensity of the trans- 
mitted ultraviolet light. In practice, Eq. 

1 must be integrated with respect to 
wavelength because of the finite pass- 
band of the ultraviolet filter. The best 




Fig. 2. (A) Normalized intensity of a star as a function of the tangent ray height 
during occultation of the star by the earth, The curves correspond to the intensity 
measured in the 1500- A and 23 90- A channels of the University of Wisconsin stellar 
photometers aboard the OAO-2 satellite [22 August 1970, 04 hours 22 minutes (local 
time). 25°N, 39°E]. (B) Number densities of 0 3 and O, as a function of tangent ray 
height determined from the occultation measurements. The solid line is the O- number 
density from the CIRA 1965 model atmosphere. 


results are thus obtained when a single 
species dominates the absorption proc- 
ess. and the OAO filters were, there- 
tore, selected accordingly. 

Once N s (r„) is known, it is a simple 
matter to invert the data and obtain the 
vertical number density profile of the 
absorbing species at the occultation tan- 
gent point, A simple geometrical argu- 
ment (Fig. I) shows that N,{r u ) along 
the ray path for a spherically stratified 
atmosphere can be written as 


A ')(/■„) 


I’ n t (r) rdr 

J (r* - 

r o 


( 2 ) 


where «,(r) is the number density of 
the ith absorbing species at radius r. 
Equation 2, the Abel integral equation 
(/, 2), is easily inverted to give the 
number density of the absorbing spe- 
cies at a tangent ray height r 


n ; ( r) 


*Lr±fS_ dro l 
dr L r„ (/■„“ — r) v ’ J 


(3) 


Thus, the stellar occultation technique 
can be used to obtain the vertical 
density profile of any absorbing atmo- 
spheric species which can be spectrally 
isolated. 

Hays and Roble (2) have calculated 
the tangential ultraviolet transmission 
of the earth’s upper atmosphere. Their 
results show that both the strong atmo- 
spheric absorption of CX in the Schu- 
mann-Runge continuum near 1500 A 
and the strong atmospheric absorption 


of 0 :} in the Hartley continuum near 
2500 A occur in spectrally isolated 
regions in which the stellar ultraviolet 
absorption is primarily due to a single 
species. The absorption cross sections 
of CX and 0 3 at these wavelengths 
have a peak around 10- 17 cm 2 . 
Therefore, we are able to determine 
the distribution of these species near 
altitudes at which N t is approximately 
10 1 ' cm - -. If a spectral region away 
from the peak cross section is 
utilized, one is able to observe 
higher tangential column number densi- 
ties or, equivalently, to measure the 
number density at lower altitudes within 
the atmosphere. The ultraviolet light in 
the wavelength interval from 1400 to 
1 600 A is absorbed primarily by CX at 
altitudes from 130 to 230 km. In the 
wavelength interval from 2400 to 2600 
A, the ultraviolet light is absorbed 
primarily by O h at altitudes from 60 
to 100 km (2). 

The normalized intensity data ob- 
tained during one of the many occulta- 
tions are shown in Fig. 2A for two stel- 
lar photometers having filters centered at 
1500 and 2390 A. The normalized in- 
tensity is obtained as a function of time; 
however, by knowing the star’s position 
and the orbital elements of (he satellite, 
we are able to relate time to the tangent 
ray height of the occulting star. Because 
of high-altitude absorption, the normal- 
ized intensity in the O-, channel at 1500 
A decays first. Then the intensity in the 


2390-A channel decays rapidly at alti- 
tudes where 0 ; , absorption becomes 
important. The excellent quality of the 
data obtained by the OAO-2 stellar 
photometers allows determination of de- 
tailed structure. The data in the 1500-A 
channel were inverted, and the results 
are >hown in Fig. 2B where the 0„ 
number density profile is shown as a 
function of height. For comparison, the 
O;, profile of the Committee on Space 
Research (COSPAR) CIRA 1965 (J) 
model atmosphere is shown also. These 
data are only a sample of the approxi- 
mately 20 stellar occultation scans that 
have been reduced thus far, hut they il- 
lustrate the quality of the data obtained 
from the occultation measurements. 

Most of the 0 ;1 scans obtained with 
the 2390-A and 2460-A filters aboard 
the OAO-2 satellite have a normalized 
intensity scan similar to the one shown 
in Fig. 2A. The intensity of the occul- 
ting star decreases as the starlight pene- 
trates into the atmosphere until a slight 
increase in the intensity curve occurs at 
a tangent ray height near 70 km. These 
data, when inverted, give the number 
density profile of 0 : , shown in Fig. 2B. 
The structure in the measured intensity 
curve is caused by a bulge in the night- 
time O a number density profile, with 
the peak occurring near 82 km and a 
minimum near 75 km for this particular 
scan. The stellar occultation measure- 
ments clearly define the structure of the 
nighttime O a profile at high altitudes 
where no previous measurements had 
been made. 

The occultation technique has also 
been used with the sun as the light 
source. The measurements, however, 
are restricted to sunrise and sunset, and 
a general review of the subject has 
been given by Link (4). 

P. B. Hays 
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Research, Boulder, Colorado 80302 

A, N. Shah 

University of Michigan. Ann Arbor 

References and Notes 

1. P. B. Hays and R. G. Roble. Planet. Space 
Sci. 16, 1 197 (1968). 

2. . /. Atmos. Sci. 25. 1141 (1968). 

.1. CO SPAR International Reference Atmosphere 
(North-Holland. Amsterdam, 1965). 

4. F. Link, Eclipse Phenomena in Astronomy 
(Sprlnger-Verlag, New York, 1969). 

5. The National Center for Atmospheric Research 
is sponsored by the National Science Foun- 
dation. 

10 January 1972 


Copyright 1972 by the American Association for the Advancement of Science 



APPENDIX II. 


A Technique for Recovering the Vertical 
Number Density Profile of Atmospheric 
Gases from Planetary Occultation Data 


17 


PRECEDING PAGE BLANK NOT FILMED 



Planet. Space Sci. 1972, Vol. 20. pp. 1727 to 1744. Per gam on Press. Printed In Northern Ireland 


A TECHNIQUE FOR RECOVERING THE VERTICAL 
NUMBER DENSITY PROFILE OF ATMOSPHERIC 
GASES FROM PLANETARY OCCULTATION DATA 

R. G. ROBLE 

National Center for Atmospheric Research,* Boulder, Colorado 80302, U.S.A 

and 

P. B. HAYS 

Department of Aerospace Engineering, University of Michigan, 

Ann Arbor, Michigan 48104, U.S.A 

{Received 29 March 1972) 

Abstract— The occultation technique of determining the properties of the atmosphere using 
absorption spectroscopy is examined. The intensity of a star, in certain atmospheric 
absorption bands, is monitored by a satellite tracking the star during occultation by the Earth’s 
atmosphere. The intensity data in certain wavelength intervals, where absorption is attributed 
to a single species, are related to the tangential column number density of the absorbing species 
through Beer’s law. The equation for the tangential column number density is the Abel 
integral equation which is inverted to obtain the number density profile of the absorbing 
species at the occultation tangent ray point. Two numerical schemes for inverting the Abel 
integral equation for signals of low intensity with statistical noise superimposed are presented; 
one for determining the number density profile of atmospheric species that decrease exponen- 
tially with height, and the second for determining the profile of constituents having a more 
complex vertical structure, such as ozone. The accuracy of retrieving the number density 
distribution from planetary occultation data is examined. A theoretical analysis of the errors 
in determining the number density from occultation data of very low signal intensity is also 
presented. The errors in retrieving the number density profile are related to the intercity 
of the source, the number of data points per scan, and the degree of data smoothing required 
before inversion. As a specific example, calculations are made of the errors in retrieving the 
molecular oxygen and ozone number density profiles from occultation intensity data in the 
Schumann-Runge continuum of molecular oxygen at 1450 A and the Hartley continuum of 
ozone at 2450 A. 

1. INTRODUCTION 

Absorption spectroscopy is an important technique for determining the composition and 
number density profiles of constituents in the upper atmosphere (Watanabe, 1958; 
Hinteregger, 1962). Most of the measurements have been made from rockets monitoring the 
absorption of the solar u.v. flux in certain atmospheric bands (Johnson et al , 1951; 
Bryam et al, 1957; Kupperian et a!., 1959; Jursaefa/., 1963, 1965; Carver et al, 1964, 
1966; Weeks and Smith, 1968; Opal and Moos, 1969; Wildman et al, 1969; Quessette, 
1970; Brannon and Hoffman, 1971). Satellites have also been used to measure the intensity 
of certain u.v. sources (i.e., Sun, Moon, stars) during occultation of the source by the 
Earth’s atmosphere (Fig. 1). The intensity data measured during occultation are then used to 
obtain information about the absorbing species in the upper atmosphere (Venkateswaran et 
al, 1961 ; Rawcliffe et al., 1963; Miller and Stewart, 1965; Thomas etal., 1965; Fesenkov, 
1967; Thomas and Norton, 1967; Norton and Wamock, 1968; Stewart and Wildman, 
1969; Hinteregger and Hall, 1969; Lockey et al, 1969; Link, 1969; Reid and Withbroe, 
1970; May, 1971; Reid, 1971; Roble and Norton, 1972; Hays et al, 1972). 

Basically, the occultation technique is similar to the classical technique of absorption 
spectroscopy. The Sun, Moon, or a star are used as a source of light and the satellite 

• The National Center for Atmospheric Research is sponsored by the National Science Foundation. 
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Fig. J. Geometry of stellar occultation. / 00 (A) is the intensity of the unattenuated 

STAR ABOVE THE ATMOSPHERE AND 7(A, r 0 ) IS THE MEASURED INTENSITY OF THE STAR AT TANGENT 

RAY HEIGHT r 0 . 

photometers as a detector, with the atmosphere between acting as the absorption cell. As 
the satellite moves in its orbit, the source ultimately is occulted by the Earth. During the 
occultation process, the intervening atmosphere absorbs progressively more of the light due 
to the strong photo-absorption features of the atmospheric gases. The occultation intensity 
data can be related to the tangential column number density of the absorbing species if the 
absorption cross-sections are known. The equation for the tangential column number 
density of the absorbing species is the Abel integral equation which is easily inverted to give 
the number density profile at the occultation tangent point. The occultation technique 
works best for a point source, such as a star, where the vertical resolution at the tangent 
point is very sharp (Hays and Roble, 1968b). For a finite-size source, such as the Sun, the 
light rays emitted from the upper and lower limb can have tangent ray heights separated by 
about 26 km when observed from normal satellite altitudes. This complicates the problem 
of relating the transmitted solar intensity to the tangential column number density (May, 
1971 ; Roble and Norton, 1972). However, once the tangential column number density is 
obtained, it again becomes easy to invert the Abel integral equation and obtain the local 
number density at the tangent ray point. For simplicity, we therefore assume that the light 
source is small, such as a star (Hays et aL, 1972) or a small portion of the solar disk 
(Hinteregger and Hall, 1969; Reid and Withbroe, 1970; Reid, 1971). 

In the present paper, we describe a technique for retrieving the number density profile of 
the absorbing species from occultation intensity data. The data-reduction procedure is 
similar to the method of computing the radial emission distribution of a cylindrical plasma 
source, as treated by the plasma physicists. The similarity exists because both techniques 
reduce to the problem of solving the Abel integral equation (Hays and Roble, 1968a, b; 
Barr, 1962). 

In practice, the numerical solution of the Abel integral equation requires special 
handling because of its sensitivity to small random errors in the data. Therefore, we 
develop two separate techniques of numerically inverting the Abel integral equation to 
obtain the number density distribution of the absorbing species. The first utilizes an ex- 
ponential approximation for atmospheric species which decrease nearly exponentially with 
altitude; the second, a quadratic approximation which works better for atmospheric 
species having a more complex vertical distribution such as atmospheric ozone (Hays et al., 
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1972). The retrieval accuracy of both schemes is examined, and a theoretical analysis of 
errors is presented which enables us to determine the accuracy of the measurement in the 
presence of random statistical errors inherent in light sources of very low intensity. The 
inversion technique and analysis of errors, although examined with respect to molecular 
oxygen and ozone, are general and can be used to determine the accuracy of any occultation 
measurement if the signal level and data-acquisition rate are known. 

2. OCCULTATION TECHNIQUE 

The intensity of a star, as measured by a telescope aboard a satellite, during occultation 
of the star by the Earth’s atmosphere (Fig. 1) is expressed as 


F*(r„) = ?- [”r(A)/(^ r 0 ) <U 0) 

where r 0 is the tangent ray height of the light ray passing tangentially to the Earth’s surface, 
(p = AQ e T 0 , A is the aperture area of the telescope, Q, is the quantum efficiency of the 
photo-multiplier, T 0 is the optical transmission of the system, T(X) is the broadening 
function or transmission of the dispersive element at wavelength X and I(X, r 0 ) is the 
intensity of the star at wavelength X and tangent ray height r 0 . 

The intensity of starlight passing through the Earth’s upper atmosphere is determined by 
considering the combined attenuation of all of the absorbing species along the ray path 
(Beer’s law). 

I(X, r 0 ) = I m (X) exp ^-2 (2) 

where 7^(2) is the stellar intensity above the atmosphere and a ( {X) is the absorption cross- 
section of the ith atmospheric species. In the case of a spherically stratified atmosphere, the 
tangential column number density A^(r 0 ) is 


N<(r 0 ) = 


f 00 nj(r)r dr 
L Vr 2 - r 0 2 


( 3 ) 


where n,(r) is the number density of the rth absorbing species at a distance r from the center 
of the Earth. 

Hays and Roble (1968b) calculated the tangential u.v. transmission for the mesosphere 
and the lower thermosphere. They show that in certain specific wavelength intervals the u.v. 
starlight is absorbed by a single species. These include a broadband region near 1500 A 
where molecular oxygen acts as a single absorber in the lower thermosphere and also a 
region near 2500 A where ozone is a single absorber in the mesosphere. In both of these 
regions, Hays and Roble (1968b) show that refractive attenuation, Rayleigh scattering, and 
absorption by other minor atmospheric constituents are negligible over most of the altitude 
range where absorption occurs. In addition, they point out that only the relative intensity is 
important for the occultation measurement. The stellar occultation data are therefore 
normalized as 

rT(X)I„(X) exp (-<r(*)JV(r 0 )) cU 

F(r a ) = ^ . (4, 

J o W)il 

The interference filter used in the occultation measurements is selected so that T(X) 
includes only a wavelength region where absorption is caused by the single species of 
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interest. By knowing the filter-broadening function, the stellar spectra above the atmos- 
phere, and the absorption cross-section, the normalized signal is then only a function of the 
tangential column number density ^F(r 0 ) = F(N(r 0 )). The intensity data can be inverted to 
give the tangential column number density N(r 0 ) = N(F(r Q )). Equation (3) for the tangential 
column number density is the Abel integral equation (c.f. Hays and Roble, 1968a, b) which 
is directly inverted to give the number density of the absorbing species at radius r 

--If (5) 

J ' s/r 0 ‘ - r s 

In the next sections, we develop two separate numerical inversion schemes; one applicable 
for determining the number density profile of atmospheric species decreasing nearly 
exponentially with altitude, such as molecular oxygen; and the second for atmospheric 
species having a more complex number density profile, such as ozone. 


3. DATA REDUCTION 

(a) Exponential form. The intensity of a star measured by photometers aboard a 
satellite during occultation is generally weak and usually obtained in digital form (Hays et 
al., 1972). At low intensities, pulse-counting systems are used with a finite integration time 
between data points. Therefore, the intensity data are obtained at a series of discrete 
tangent ray heights F(r 0 ). The normalized intensities are converted to the tangential column 
number density of the absorbing species along the ray path N(r 0 ) as described in the previous 
section. Thus, a set of discrete values of the tangential column number density is obtained as 
a function of tangent ray height during occultation. We approximate the tangential column 
number density between data points in the interval r { < r < r m by 

N(r) = «. exp (— ft.(r - r t -)) (6) 

where r t and are the i and i -f I tangent ray heights of the data points, respectively. The 
coefficients a,- and /?, are determined from a least squares fit to the data in the vicinity of the 
/th data point by minimizing the variance of Equation (6) with the data and smoothing M 
adjacent points (c.f. Roble and Norton, 1972), 


& i = 2 - «< exp (-ft(r, - r,))) 2 (7) 

where A x = / — Mj2, K 2 = i + Af/2, N } is the tangential column number density at they'th 
data point. Then by differentiating Equation (6) and replacing the integral in Equation (5) 
by a finite sum, we obtain 


7T»=l Jrt yj r * _ r f 

Expanding the denominator and integrating, we obtain 


( 8 ) 


nfr,) -if -SfiL-J (, + I ^ _ _L _i_\ M 1/2 

+ 2 (r ( + r,) 4ft (r, + r,)) \fij 

X [erf (ft(r, +1 - r,)) 1,a - erf (ft(r< - r,)) 1 ' 2 ] exp (ft(r, - r,)) 

+ 2ft(r ( 1 + r,) [(r ' +1 “ CXP ( ~ A(r <+> “ r <)) - <T< ~ '•i) 1 ' 2 ]} 
where erf(x) is the error function of x. 


( 9 ) 
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Thus, the number density of the absorbing species at radius r t is determined directly 
from Equation (9) using the calculated values of a, and 

(b) Quadratic form. When the absorbing constituent has a number density height profile 
which is not exponentially decreasing with altitude, but varies over the scan altitude, such as 
for ozone (Fig. 2), then the exponential scheme developed in the previous section may not be 


Oj NUMBER DENSITY (cm"’) 

| 0 6 | 0 T I0 B 10* I0 w 10” 



Fio. 2. The molecular oxygen and ozone number density as a function of raiGHT. 
The solid curve for O, is the number density distribution from the mean CIRA 1965 
model atmosphere and the dashed curve is an analytic approximation to the distribution , 
The 0 8 number density distribution is an analytic approximation of a night-time 

PROFILE AND IT IS GIVEN BY EQUATION (17). 


best suited for the data-reduction process. In this case, we approximate the tangential 
column number density in the altitude region z, < z < z f+1 by the quadratic relationship 

N(z) ==Az* + Bz + C (10) 

where z = r — r t and r e is the radius of the Earth, 6371 km. 

The coefficients A, B and C are again determined by minimizing the variance of 
equation (10) with the data and smoothing M adjacent points 
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Minimizing Equation (11) and solving directly for A it B i and C { we obtain 

A = - (N)(z))(M(?) - (f)(7)) - ^ ~ (z)(z)) 

(M(z 3 ) - (z a )(f))(Af(z 3 ) - (z)(z*)) - (M(z 4 ) - (z 2 )(z 2 ))(Af(7) - (z“)(f)) 

„ (M(iVz) - (z)(W)) - (M(z 3 ) - (z 3 )(z)H 

*** “ “ (12) 

(Af(z 2 ) — (z)(f» 

c ^ W - (z 2 )*i - (m 

* M 

where 

K* 

r ) = i ( v 

As in the previous section, by differentiating Equation (10) and replacing the main integral of 
equation (5) with a finite sum and then integrating between data points, we obtain 


a(r £ ) - 



r ( a ) l/a ~ W - rtf") 


+ (B, - Iog e 


(ri+i + (r? +1 - r t a ) 1/a ) ~[ 

( r t + ( r < a — r j 2 ) 1/2 ) J ^ ' 


Thus, the number density of the absorbing species as radius r t is determined from equation 
(13) using the coefficients and B { calculated for each data point. 


4. ATMOSPHERIC MODEL 

The molecular oxygen number density profile in the lower thermosphere, obtained from 
the mean CIRA 1965 model atmosphere, is shown in Fig. 2. For mathematical simplicity, 
we approximate this profile by an atmosphere with a constant scale height gradient 
(Nicolet, 1960), 

H = H, + P{r- O (14) 


where H a is the scale height of the atmosphere at radius r s , and is the vertical gradient of 
the scale height. The molecular oxygen number density is then represented as 


, . r n 

n(r) = n. 

*lH, + £(r - r s )J 


(i 


(15) 


where n a — 2 X lO^cm -3 , r 3 = 6471 km, H a = 5 km, and = 023 to give the ap- 
proximation to the CIRA 1965 molecular oxygen distribution shown in Fig. 2. The 
tangential column number density of molecular oxygen along the ray path is determined by 
inserting Equation (15) into Equation (3), expanding the denominator, and integrating to 
obtain 


w(ro) _ (ip r — ^ — nr — t — rvi 

\rj Lh, + ^(r 0 - r,)J \ H, + - r.)J \2 \d ill 


+ r fc 7 fc^)] B (l’ (l“l))) (16) 


where <3 = >3/(1 + /?) and B(x, y) is the Beta function, B(x , ;;) = y(x)y(y)ly( x + y) and 
y(.v) is the Gamma function. The normalized intensity as a function of tangent ray height is 
calculated from Equation (2) and is shown in Fig. 3 for absorption at a single wavelength 
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1450 A in the Schumann-Runge continuum where the molecular oxygen absorption co- 
efficient is approximately 2 x 10" 17 cm 2 . 

Hays et al. (1972) presented a typical night-time ozone distribution determined from 
occultation measurements made by the Orbiting Astronomical Observatory (OAO-2). The 



Fig. 3. The normalized tangential transmission of the Earth’s atmosphere at various 

WAVELENGTHS AS A FUNCTION OF HEIGHT. 

ozone number density profile given in their paper can be reasonably represented as a sum of 
exponential and Gaussian functions 

n(r) = n b exp (— (r - r b )jH b ) + n e exp ( — ((r — r^/tf,) 2 ) (17) 

where n b = 5 x 10 10 cm” 3 , r b = 50 km, H b = 4-34 km, n e = 8 x 10 8 cm- 3 , r c = 83 km, 
and H e = 5 km. The ozone number density distribution obtained from equation (17) is 
shown in Fig. 2. 

The coefficients used in this study give an ozone number density slightly higher than those 
for a typical night-time ozone distribution determined by Hays et al. (1972). This difference, 
however, is not important because the profiles are similar and the main features of the ozone 
analysis are preserved in either case. 

The total columnar number density along a ray path is obtained by inserting equation 
(17) into equation (3), expanding the denominator, and integrating to give 

N(r 0 ) — 2n b exp (rJfft)r 0 J? 1 (r 0 /J7,) + 2n 0 exp (— ((r„ — r c )fH t f ) 

X (2r 0 r 1,s J yii) exp ((r„ - r,) 8 /(2/f/)) 

X - rj) + 

X exp((r 0 - r.)7(2H t ! ))D_ 1/a ((2) 1 ' !! /f o - , (r 0 - r e ))J (18) 
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where /^(x) is the modified Bessel function of the first kind, y(x) is the Gamma function, 
and D_ v (x) is the parabolic cylinder function. The normalized intensity as a function of 
tangent ray height is obtained from Equation (18) and Equation (2), and is shown in Fig. 3 
for absorption at a single wavelength in the Hartley continuum of ozone at 2500 A. At this 
wavelength, ozone has an absorption cross-section of about 10~ 17 cm -2 . 

5. SAMPLING ERROR 

The accuracy in retrieving the number density profile with the data-reduction scheme 
depends upon (a) the data-acquisition rate or number of data points per scan, (b) the 
statistical noise inherent in a signal of low intensity, (c) the smoothing parameter M used in 
the inversion process, (d) knowledge of the satellite position during occultation, and (e) 
departures from spherical symmetry. The first three items are examined in this paper. Item 
(d) depends on the satellite tracking network and the degree of accuracy of the orbital 



Fig. 4. The error in retrieving the analytic molecular oxygen number density as a 

FUNCTION OF HEIGHT FOR VARIOUS DATA-ACQUISITION INTERVALS. 


elements. The intensity of the star is measured as a function of time, but by knowing the 
star’s position and the satellite orbital elements, we can relate the intensity to the tangent ray 
height of the occulting star. Thus, errors in satellite time and position appear as errors in 
the tangent ray height and will not be discussed in this paper. The last item may be im- 
portant for certain atmospheric species which change rapidly with time, as perhaps during 
sunrise or sunset; however, in this analysis, spherical symmetry is assumed. 

The error in retrieving the molecular oxygen number density profile using the ex- 
ponential data-reduction scheme is plotted as a function of height in Fig. 4 for various 
data-acquisition rates. The tangential column number density for the analytic molecular 
oxygen distribution is calculated at the various altitudes using Equation (16). These data are 
inverted using the exponential data-inversion scheme, and the retrieved molecular oxygen 
number density is compared to the number density given by the analytic model. The errors 
in retrieving the number density are then plotted in Fig. 4. The results show that the errors 
are smallest when a large number of data points per occultation scan are used in the 
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inversion process. The smoothing parameter H is also a factor in the retrieval accuracy. 
The smoothing parameter is primarily used to smooth the analytic representation of the 
tangential column number density through noisy data points and its importance will be 
evident in the next section. But for noise-free data, an increase in the smoothing parameter 
decreases the retrieval accuracy for a number density profile having a slight curvature. 

At high altitudes, the retrieval accuracy would decrease due to the practical requirement 
of truncating the inversion integral at an upper boundary. It is difficult to convert the 
normalized intensity to tangential column number density for the upper part of the scan 
where the normalized intensity approaches 1*0, especially in the presence of statistical noise. 

100 
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-B -4 O 4 8 12 16 -20 24 

ERROR IN RETRIEVING O s NUMBER DENSITY 1%) 

Flo. 5. The error in retrieving the analytic ozone number density as a function of 

HEIGHT FOR VARIOUS DATA- ACQUISITION INTERVALS. 

To reduce this error and extend the upper limit of the occultation scan, it is convenient to 
utilize a model atmosphere and calculate the tangential column number beyond the upper- 
most data point for use in the inversion integral. 

The error in retrieving the model ozone number density profile using the analytic 
expression for the tangential column number density given by Equation (18) and also using 
the quadratic data-reduction scheme is shown in Fig. 5. The retrieval is good for a high 
data-acquisition rate, but as the data-acquisition rate decreases, the accuracy in defining the 
vertical structure and in particular the ozone bulge near 83 km also decreases. The lowest 
data-acquisition rate considered in Fig. 5 show's considerable error in retrieving the ozone 
bulge. These results were obtained using a smoothing parameter of M = 2. As the smooth- 
ing parameter is increased, the accuracy of retrieving the vertical structure depends critically 
on the data-acquisition rate. For a particular application, there is a tradeoff between the 
data-acquisition rate and the smoothing parameter. 

For simplicity, the results presented here assume absorption at a single wavelength. 
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However, the technique which has been developed is general and can be applied to any 
broadband filter by including a wavelength integration. 

6. THEORETICAL ANALYSIS OF STATISTICAL ERRORS 
For occultation scans of very low intensity, the statistical noise which is superimposed 
upon the basic signal intensity becomes an appreciable part of the signal. The errors in 
retrieving the number density profile of the absorbing species from occultation data having 
random statistical noise are considered in this section. 

(a) Molecular oxygen. The solid curve in Fig. 6 shows the calculated normalized 
intensity variation of a star during occultation using Equations (16) and (2). The peak signal 



Fio. 6. An occultation scan in the Schumann-Runge continuum of molecular oxygen 
AT 1450 A WITH RANDOMLY SCATTERED STATISTICAL NOISE SUPERIMPOSED. THE STANDARD 
DEVIATION OF THE NOISE IS PROPORTIONAL TO THE SQUARE ROOT OF THE BASIC SIGNAL LEVEL. 

level in this example is approximately 1000 counts sec -1 and the calculations are made at a 
single wavelength of 1450 A. Also shown in Fig. 6 is the basic scan with a superimposed 
Poisson noise distribution having a standard deviation proportional to the square root of the 
signal level. The molecular oxygen number density profiles shown in Fig. 7 for the two scans 
given in Fig. 6 have been retrieved using the data-reduction scheme described in the previous 
sections. Other statistical noise distributions yield an ensemble of retrieved molecular 
oxygen number density profiles which can be statistically analyzed. 

If we assume that the tangential column number density is N — N 0 + AN where N 0 is 
the basic tangential column number density without statistical noise and AN is the tangential 
column number density error due to the statistical noise, then in the approximation given by 
Equation (6) we expand a and as 

a = a 0 + Aa 

+ A$. (19) 

Inserting these relationships into Equation (7), expanding the exponential, and subtracting 
the basic unscattered state, we obtain 

K t 

V = 2 [AN, - Aa, exp (-/3 0 ,(z, - z,)) + a 0 ,AA(z, - z,) exp (-/?„,(*, - z ( ))f (20) 

i 



DENSITY PROFILE OF ATMOSPHERIC GASES 


1737 



Fig. 7. The retrieved molecular oxygen number density profile for the occultation 

SCAN WITH RANDOMLY SCATTERED NOISE SHOWN IN FlG. 6. TlIE SOLID CURVE GIVES THE ANALYTIC 
MOLECULAR OXYGEN NUMBER DENSITY DISTRIBUTION. 

j 

The coefficient errors A a* and are calculated in the least squares sense from the 
relationships 


ddf d&S 

dAotj 3A/?, 


( 21 ) 


and are given in the Appendix. 

The intensity of the star, considering single wavelength absorption, is calculated from 
Beer’s law 


Ft = F{X, r 0 ) = FM exp (-<r(A)W(r 0 )) (22) 


where i r w (2) is the signal measured above the atmosphere in counts sec -1 . The linear 
variation in tangential column number density due to a signal with statistical noise F = 
F q + AF is obtained by introducing the value of N into Equation (22) and expanding the 
exponential to obtain 


AN = 


1 AF 
<j(A) F(X, r 0 ) 


(23) 


If a Poisson random noise with a standard deviation for A F proportional to the square 
root of the basic signal ( F 0 ) m is introduced into the expressions for the coefficients Act, and 
A(i f , the mean square errors for the coefficients (Aa*) and (A/9,-) are determined using the 
propagation of statistical errors theorem and they are presented in the Appendix. 

The errors in the number density of the absorbing species at the tangent ray point are 
determined from the integral 


An(r t ) = 


1 f 09 [d ANjdr ] dr 
v h Vr 2 — r ( 2 


(24) 


The tangential column number density of the absorbing species is determined by 
inserting Equation (19) into Equation (6). The noise in the number density of the absorbing 
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species at the tangent ray point is obtained by expanding the exponential, subtracting the 
basic state, and differentiating the tangential column number density error A N to give 


= - I (ftiAai 


7 T i~l 


fri + 1 

» + 


exp(-^ 0i (r - r<)) 


vV- 


dr 


Integrating, we obtain 


^ dr. (25) 

7Ti-=i *h< % / r 2 _ r a 


An(r,) = -2(R,A«, + S ( AA) 


7M=.l 


(26) 


where i?, and S { are coefficients given in the Appendix. The standard deviation of the 
number density error at the tangent ray point due to statistical noise is calculated using the 
propagation of statistical errors theorem 

<An(r,)> = T, 

where 


(27) 


T, = L 


'00 — 1 1/2 

I (R,X, + SM* 

Li=I 


and and Y { are given in the Appendix. 

(b) Ozone. The errors in the ozone number density caused by statistical noise super- 
imposed upon the basic signal are determined from the quadratic data-reduction scheme. 
The signal intensity variations are related to the tangential column number density variations 
through Equation (23) and the number density error of ozone at the tangent ray point is 
obtained from Equation (24). The tangential column number density with superimposed 
variations caused by statistical noise is N = A 0 -f AN. The coefficients in Equation (10) are 
assumed to have a variation A t = A oi -f A A ir B t = B 0i + A B {y and C< = C 0i + A C { . 
Introducing these variations into Equation (11) and subtracting the basic state, we have 


Kt 


= I (A N, - AA iZ * ~ A B i z f - AC,) 2 . 


i-Ki 


(28) 


Solving for the coefficients A A 0 A and AC,-, we obtain relationships similar to Equation 
(12) with the exception that N f is replaced by A N } 


A A t = 2 Qi AN, 

i<*Ki 

AB, = f Q* AN, 

}^K i 


(29) 


where Q } and Q { * are given in the Appendix. The number density error at the tangent ray 
point is obtained by inserting Equation (29) into Equation (24) and integrating to give 

i 00 r Ka A F "1 

Mr,) = - ~ I I (( 2W U - 2 r.WJQt + ^,8,*) . ^ (30) 

TTi=lLf-Kl CfF. i, J 


where 


Wu = (r? +1 - r,T 2 - (r, 2 - r, 2 ) 1 ' 2 
/ '•i+i + Oja - r, 2 ) 1 /2 

' r i + (r j 2 ~r,r i 


= log, 


)■ 
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The standard deviation of the number density error at the tangent ray point is calculated 
using the propagation of statistical errors theorem and assuming that the standard deviation 
of the noise is proportional to the square root of the basic signal 

<A rt(r I )) = £,.i^ 1/2 (31) 

where 

E, = -{I I [ [(2 W u - 2 r.WWCi + *] ;(~) l, T) 7 ■ 

i-ATil <*'r 0j / J ) 

♦ 

7. ANALYSIS OF STATISTICAL ERRORS 


The error in retrieving the molecular oxygen number density from occultation scans 
having a random statistical noise superimposed upon the basic signal is defined as 


e(r,) = 


(A»(r t )) 


(32) 


The numerator is the standard deviation of the retrieved number density, given by 


equation (27) for molecular oxygen and equation (31) for ozone, and n(r t ) is the number 
density at radius r v The standard deviation of the retrieved number density is inversely 
proportional to the square root of the unattenuated signal level for both molecular oxygen 
and ozone. Therefore, in Fig. 8 r^r,))- 1 is plotted as a function of altitude for various 
smoothing parameters M . The error in retrieving the molecular oxygen number density for 



Pig g. Theoretical analysis of errors in retrieving the molecular oxygen number 

DENSITY FROM OCCULTATION SCANS WITH RANDOM STATISTICAL NOISE SUPERIMPOSED. THE 
RETRIEVAL ERROR IS OBTAINED BY DIVIDING THE PARAMETER T^flO,))- 1 BY F*' 1 , THE SQUARE 
ROOT OF THE UNATTENUATED BASIC SIGNAL LEVEL. THE CIRCLES REPRESENT THE RESULTS OF A 
NUMERICAL ANALYSIS OF ERRORS FOR M = 4. 
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any signal level is obtained by dividing the parameter r,(n(r,))- 1 by F l £\ An increase in the 
smoothing parameter M decreases the error in retrieving the molecular oxygen number 
density. However, as discussed earlier, the accuracy of retrieving the basic molecular 
oxygen number density profile from unscattered occultation scans decreases as the smoothing 
parameter M increases. Therefore, there exists a value of M which gives a minimum in the 
overall error in retrieving the molecular oxygen number density profile at a given data- 
acquisition interval. The results of a numerical analysis of errors are shown in Fig. 8 for 100 
occultation scans with different randomly scattered noise distributions. There is general 



Fio. 9. Theoretical analysis of errors in retrieving the ozone number density from 

OCCULTATION SCANS WITH RANDOM STATISTICAL NOISE SUPERIMPOSED. THE RETRIEVAL ERROR 
IS OBTAINED BY DIVIDING THE PARAMETER («(r j))" 1 BY F"*, THE SQUARE ROOT OF THE UNATTENU- 
ATED BASIC SIGNAL LEVEL, THE CIRCLES REPRESENT THE RESULTS OF A NUMERICAL ANALYSIS OF 

ERRORS FOR M — 4. 

agreement between the results of the theoretical and numerical analysis of errors. Also 
shown in the figure are arrows indicating the altitudes of the 0-1 and 0-9 normalized intensity 
points. The retrieval of the molecular oxygen number density is best for data lying between 
these normalized intensity ratios. 

The standard deviation of the retrieved ozone number density is also inversely pro- 
portional to the square root of the unattenuated signal level. Therefore, in Fig. 9 the 
parameter E t (n(r{)) 1 is plotted as a function of height for various smoothing parameters. 
In general, the error in retrieving the ozone bulge and the number density below about 
70 km is low. However, in the region of the ozone minimum near 75 km the retrieval error 
increases. The data-acquisition interval in this case is 1 km and an increase in the smoothing 
parameter M decreases the error in retrieving the ozone number density. As the data- 
acquisition interval increases, the error in retrieving the ozone number density increases. 
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For a given data-acquisition interval, there also is a value of the smoothing parameter M 
which gives a minimum in the overall error in retrieving the ozone number density. The 
circles in Fig. 9 again are the results of a numerical analysis of the errors for 100 ozone 
occultation scans with different randomly scattered noise distributions. The general agree- 
ment between the theoretical and numerical analysis of errors enables us to theoretically 
calculate the error coefficients and thus determine accuracy of an occultation scan for any 
signal intensity. 

8. DISCUSSION 

We have presented a numerical technique for retrieving the number density profile of 
absorbing atmospheric gases from planetary occultation data having random statistical 
noise superimposed upon the basic signal. Because of the sensitivity of the Abel integral 
equation to random noise, some smoothing is required before evaluating the inversion 
integral. An exponential technique was developed for atmospheric species whose number 
density profile decreases exponentially with height, whereas the quadratic scheme works 
better for atmospheric species having a more complex vertical profile. These techniques have 
been used to reduce the OAO-2 stellar occultation data and obtain the night-time molecular 
oxygen and ozone distribution in the upper atmosphere (Hays et at . , 1972). With a modi- 
fication to account for the finite size of the solar disk, these techniques have also been used 
to reduce the SOLRAD-8 solar occultation data (Roble and Norton, 1972). The retrieval 
errors are determined by knowing the approximate signal level and acquisition rate of the 
occultation data. For each specific case, there is a minimum in the retrieval error determined 
by the degree of smoothing for the appropriate data-acquisition rate and basic signal level. 
Generally, the best accuracy is achieved between the normalized intensity limits of 0- 1-0*9 
where the derivative in the numerator of the Abel integral equation changes rapidly. Figure 2 
shows that by selecting the appropriate spectrally isolated wavelength interval, the number 
density distribution at practically any altitude may be obtained in the region where the 
tangential column number density is approximately equal to the inverse of the effective 
cross-section. 

The analysis was performed for single wavelength absorption, whereas measurements 
are made for some finite passband determined by the characteristics of the monitoring 
spectrometer or photometer. In this case, an additional integration over wavelength is 
required in the analysis, but it does not fundamentally alter the results as long as the 
absorption process is spectrally isolated. Actually, a wavelength interval selected to cover a 
continuum region where the absorption cross-section is smoothly increasing or decreasing 
allows deepest penetration into the atmosphere. The altitude range of the occultation data is 
greater in this case than for the case of a constant cross-section over the selected wavelength 
interval. 
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APPENDIX 

Several of the coefficients discussed in the text are presented in this Appendix. 
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Appendix III. a. 

This Appendix gives some details of data processing which is 

employed to obtain molecular oxygen profiles and ozone profiles from 

occultation data. The details are divided into three parts, or processes, 

dealing with (1) the computation of tangential column number density for 

molecular oxygen, N(0 9 , r ); (2) the computation of tangential column 

z o 

number density for ozone, N(Og, r Q ); and (3) the retrieval of number density 
of either gas. Computer programs with accompanying notes are included. 
Process 1 

To obtain N(O q , r ), the raw modulus 256 occultation data on 
z o 

the strip chart or digital printout are converted to normalized signal inten- 
sity by a program filed in "PLOT". The input to this program is the 
sequence of counts taken at an interval of J sec spacecraft time. The 
measure of intensity which is used is the difference between successive 
counts. The output consists of a graph of normalized signal intensity v. 
time and the corresponding punched cards. 

Following the conversion of the data to normalized intensity, 
a program filed in "OAOPR" is activated to complete the process. First, 
the ray tangent point location is computed. Then an interpolation table is 
formed consisting of normalized signal intensity v. tangential column 
number density. The table is then used to find N(0 2 , r^) corresponding to 
the normalized, measured signal intensity. 

Process 2 

To obtain N(O n , r ), the raw occultation data on the strip chart 

O O 

or digital printout are converted to normalized signal intensity by the "PLOT" 

programs and the "OAOPR" program is activated as for Og. However, in 

this case corrections to N(O q , r ) are made for the effects of the 0„ Herzberg 

z o z 

bands and Rayleigh scattering upon ozone. 
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Process 3 


Once the tangential column number density, N, is computed 
l'or all m data levels, a parabolic fit is forced upon N(h) and the smoothed 
values are then inverted by the finite differencing scheme. 


t \ 2d 
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J 7T dr 


f -jr- <r 2 - 

J. r o 0 


r 2 ) N{ r ) dr? 
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Details of the smoothing technique and inversion for various species 
is given in the attached paper on inversion techniques. 
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SSIG SAFI?' C 

= 510 

•MAURICE GRAVES 1 

**LAS? SIGNON * AS : 09:36.59 12-29*7 1 

USER "SAFE 

» SIGNED ON AT 21:18.13 ON 12-29-71 

SLIST PLOT 



1 


DIMENSION AY (150) , Y { 150) 

2 


DIMENSION BY (150) 

2. OS 


DIMENSION YI NT (150) 

2.1 


READ (5,86 1) (BY (K) , K = 1 , 8 0) 

2.15 


PRINT 861, (BY (K) ,K= 1, 60) 

3 


READ (5, 894) N , MAX 

3.5 


PRINT 894, N,MAX 

4 


READ 893, (AY (I) ,1=1, N) 

4. 1 


LMT= 5 

4.2 


PRINT 893, (AY (I) ,1=1, N) 

5 


DO 882 K= 1 , LMT 

6 

882 

Y (K) = AY (K) 

6. 1 


JA=LMT-1 

6. 2 

• 

DO 1 K= 1 , JA 

6.3 

1 

YINT (K) = Y (K) -Y (K+1) j 

7 


J=N-LMT 

7.1 


WRITE (6 ,861 ) (BY (K) ,K=1,80) 

a 


WRITE (6,899) 

9 

899 

FORMAT (20X, 8HEXPECTED, 9X , 6HACTU AL, 14 X , TOHDIFFERENCE, 10X, 

10 


18HRECORDED, 12X,9UINTEMSITY,/, 22X, SHV ALUE, 9X, 5HVALUE, 35X , 

1 1 


25HVAEUE//) 

11.1 


PZ = 0. 1 5 * M A X 

12 


DO 883 K = 1,J 

12. 1 


Y(K)=AY{K) 

12.2 


DO 2 M P = 1 , J A 

12.3 

2 

Y (K+MP) =Y (K+MP-1) -YINT (K+MP-1) 

12.41 


YB= ( Y (K) + Y (K+1) + Y { K+2) + Y (K + 3) + Y(K + 4))*0.2 

12.42 


AI= (Y (K) -2. * Y (K+1) + 2. + Y (K+3)-Y (K + 4) ) / (- 12. ) 

12.43 


31= ( 4 . * Y (K) + Y (K+ 1) + Y (K+3) + 4.*Y (K + 4) -10.*Yfl) / 1 4. 

12.44 


CI= (Y (K) -8. *Y (K+1) +8 • * Y (K + 3) - Y (K+4) ) / 12. 

12.45 


B= (-9 • * AI) +BI 

12.46 


C= (27 . * AI) - (6. *81) +CI 

12.47 


D=- (27. *AI) +9. *31-3. *CI-2.*BI + YB 

12.5 


P 1=2 1 6* AI + 3 6*B+6*C + D 

12.52 


L=K + 2 

13 


A= (Y (L) -2. *Y (L+1) +Y (L+2) ) /2. 

14 


B= (- 5. * Y (L) + 8.*Y (L+1)-3.*Y (L + 2) ) /2. 

15 


C= 3 . * Y (L) -3 • * Y (L+1) + Y (L + 2) 

20 


P= 1 6 . *A* 4. *B + C 

20.2 


P= (P + P 1 ) /2. 

2 V 

, 

X=AY (K + LMT) 

22 


L=0 

22. 1 


PRINT 3,K,J,X,P,P1 

22.2 

3 

FORMAT (T2,2I5,3F20. 6) 

23 

886 

IF (A8S (X-P) -160.) 884,885,885 

24 

884 

IF ( (P-X) .GT. 100.) X=X+MAX 

24.2 


PW=0 .08*YINT (K+LMT-2) 

24. 4 


IF (PW .LT. PZ) PW=PZ 

24.42 


IF (PW .GT. 160.) PW = 160 * 

24.5 


IF ( (Y ( K+LMT-2) -Y (K+LMT-1) ) .GT. (Y (K+LMT- 1) -X + PW) ) X=X-MAX 

24.8 


Y (K + LMT) =X 

25 


B=X-P 

26 


C=Y (K+LMT-1) -Y (K+LMT) 

26.1 


YINT (K+LMT- 1 ) =C 

27 


WRITE (6,895) K,P,X,B,AY (K+LMT) ,C 

28 

895 

FORMAT (15, 5F2Q. 0) 

40 



29 

663 

CON riMJJS 

3 0 


J=N-1 

32. Ob 


WRITE (7, 861) (BY (K) ,K=1,80) 

3 J* 1 


W ill T E (7,865) (YINT(K) ,K-1, J) 

33 


WRITE (0,898) 

34 

696 

F OEM AT ( 1 1 GRAPE OF INTENSITY 

33 


EE AO (5 , 8 oO) DY,CY,EY 

3o 

6 60 

FORMAT ( 3 A 1 ) 

3 o . 1 


FEINT do 1 , (dY (K) ,K=1,80) 

37 


DO 692 K = 1, 1 2 

33 

692 

li Y (i\ ) = 1 0 * K 

3 7 


WRITE (6, 661) (BY { K) ,K=1,12) 

40 

661 

FORMAT (/3 X , 1 2F10. 0) 

4 1 


DO 690 K=1, 150 

42 

6 90 

BY (K) -C Y 

43 


DO 691 K= 1 , 2 5 

44 

691 

BY (5* (K- 1 ) + 1) - EY 

45 


WHITE (6,861) (BY (K) ,K=1, 120) 

40 


BY (1) =DY 

47 


DO 670 K-2, 1 50 

4 d 

670 

BY (K) = BY (K- 1 ) 

43 


N = N- 1 

43.1 


J B= 1 

49.2 


DO 741 J=1,N 

4 9.3 


IF ( Y I NT ( J) . GT. YINT (JB) ) JB=J 

4 9.4 

741 

CONTINUE 

50 


DO 88d K-1 , N 

51 


BY (1) — C Y 

52 


J A= K/5 

53 


J A= J A*5 

54 


IF (JA . EQ. K) BY (1) =EY 

54. 1 


L= YI NT (K) *100. /YINT (J3) 

5 o 


BY (L)=EY 

57 


PRINT 861, (BY (J) ,J=1,12Q) ,L 

5d 

661 

FORMAT ( 1 20 A 1 ,15) 

59 


BY (i-) =BY (L+1) 

60 

dad 

CONTINUE 

6 3 

S o5 

FORM AT ( 12F6. 0) 

64 

693 

FORMAT (14F5.0) 

65 

69 4 

FORMAT (215) 

6 b 

667 

CALL EXIT 

67 


RETURN 

66 

685 

L* L+ 1 

69 


X = AY (K + LHT) -MAX*L 

69. 1 


PRINT 4 ,K,L, MAX,X 

69.2 

4 

FORMAT (T2, 315, F20. 6) 

70 


IF(L-400) 886,887,887 

7 1 


END 


END OF FILE 


TIME'//) 
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Part 1 : Line numbers 1-6 

In this part, an informative comment card is read (line 2. 1), the 
number of data values (N) and the maximum counts on calibration (MAX) are 
read (line 3) f the raw data (AY) are read (line 4), and the first 5 values are 
set in Y (line 6). These five values must be monotonic estimates of values 
on the curve of actual (not normalized) intensity. 

Part 2 ; Line numbers 6. 1-12, 5) 

In this part, an allowable departure of a predicted; value of 
intensity from the actual value is established. This departure is ±0.15 , 
(line 11. 1). Then the coefficients (AI, B, C, D) for cubic curve -fitting 
with five points are calculated (lines 12-12, 47), and a sixth point is predict- 
ed (line 12. 5). 

Part 3: Line numbers 12. 52 -20. 2 

In this part, a parabolic curve -fitting is carried out for the third, 
fourth, and fifth points, yielding a new prediction of the sixth point (line 20). 
An average of the two predictions becomes the final estimate (line 20. 2) of 
the sixth point. This procedure continues within loop 883 until the entire 
set of data points is exhausted. 

Part 4: Line numbers 21-29 

In this part, the modulo conversion is accomplished, based upon 

some intuitive conditions regarding signal intensity. The resulting signal 

* 

intensities (C) are written out (line 27). 

Part 5: Line numbers 30-71 

In this part, the signal intensities are punched (line 32. 1) for the 

b v 

data input to the final processing program. The intensity values are also 
plotted (line 57) after normalization (line 54, 1). 
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iMi a :j t 1 1 >: !* a y a sol to least square in version for occult atign data 

1 ij 1 ;i F J C.ti *w\ i*i FOR S LIT I NC STAR ONLY.*** 1 ^ 
i .ii'i-it I I i'l l. A >■ '' ( A — H , 0"* 1 1 2 — L ) 

J Oil Iv, ( 1 50) , ALAMO ( 150) 

COMMON L (TOG) , \ N (200) , A (TOO) , U (TOO; , C (200) , ALT , W LL ( 1 00) ,FILT (100) 
1FIO < 100) ,SG ( 100) ,SG02 ( 100) ,SGM (100) ,CCM (2 00) , NP,N?X 
U i i" N:> i C N DA (TOO) , F O c ( T 0 0 ) , T NO ( 20 0 ) 

LIMiASlGN CF (TOO) , COT (TOO) ,r:LR(20C) 

DINE NSaGN t\ { T O 0 ) , FN (TOO) ,1N (20 0) , I*' MO ( 200) ,ZP (200) 

DIMENSION DEN (200) , Z N ( 2 0 0 ) 

DIMENSION 0(20) ,T (150) 

DIMENSION 1(0 (200) 

DIMENSION UD (20 0) , UZ (200) 

DIMENSION UT (200) 


1 C, £ * ****** * *>* + «*****£ * v £ v * * * **♦**####*;*+#******♦*****#****#*#** ** ** * * 


16 C* 


FART-1 


17 

16 

1 9 
20 

2 1 
2 2 
2 I 


C* 

C* THIS PART CALCULATES TANGENT RAY HEIGHT POINT* INPUT IS 5 VALUES 
C* OF SPACECRAFT LONG. LAX. AND ALT. IN 10 SEC. GMT INTERVAL, 5 VALUES 
C* OF GUA OF ARIES IN 10MIN. GMT INTERVAL, DEC. AND R. A. OF OCCULT- * 
C* IMG STAR, THE STARTING TIKE IN HR* , MIN., AND SEC., AND NUMBER OF 
C+ DATA POINTS NEEDED 


2 4 £************** ^^t******^****** ***##**##***«#***#********#**#********** 


2 5 


2 6 
27 

y?') 

JO 
3 1 

32 

33 
3 4 
3 5 
3o 
J 7 
3 ti 

3 9 

4 0 

4 1 

4 2 

43 

44 
43 
4 i 
4 7 
4 S 
4 9 

50 

51 

52 

53 


56 

57 
50 


Cl--; U/37. 2957755131 
WRITE (6, 109) 

109 FORMAT (//10X ( * S/C POSITION, LONG. LAT . , AND ALT. IN 10 SEC. • 

1 'STEPS AND' ,/,10X, ' GiiA. OF ARIES IN TIN MINUTE STEPS') 

2i-.AJ 11, Y1,Y2,f3,Y4, YS 

tRIMT 1 1, Y 1 ,Y2,Y3,Y4,Y5 

CALL S ET ( Y 1 , Y 2 , Y J , Y 4 , Y 5 , A 4 , li4 , C 4 , D 4) 

READ 11 , Y 1 , Y 2, Y 3, Y4, YS 

PRINT 1 1, Y 1 ,Y2,Y3,Y4,Y5 

CALL SET ( Y 1 ,Y2,Y3,Y4 / Y5,A2,B2,C2,D2) 

a E A D 1*1, Y 1 , i 2 , 7 3 , Y 4 , Y 5 

PRINT 11, Y1,Y2,Y3,Y4,Y5 

CALL SET (Y1, Y 2 , Y 3, Y 4 , Y 5 , A 3, B3 ,C3 , D 3) 

READ il,Y1,Y2,Y3,Y4,Y5 

PRINT 1 1, Y1,Y2, Y3,Y4,Y5 

CALL SET (Y1 , Y 2 , Y 3 , Y 4 r Y 5 , A 5 , 135 , C5, D5) 

2 i-AD 11,Y1, Y2,Y3,Y4, Y5 
PRINT 71, i 1 , Y2, 7 3 

71 FORMAT (///5X, 1 STAR DECLINATION* ,?5. 1 , * DEGREES • , F5. 1 , ' MINUTES', 
IF 7. 2 , 'SECONDS'/) 

Li:C=Y1+ ( Y2/6 0 . ) + ( Y 3/3600. ) 

READ 11, Y 1, Y2, Y3, Y4, YS 
PRINT 72, Y1,Y2,Y3 

72 FORMAT (5a, ' STAR RIGHT ASCENTION * , F5. 1 , ' HOURS ' ,F5. 1, ' MINUTES', 
If 7.2, 'SECONDS'/) 

R A- (Y 1 + ( Y 2/6 0 . ) * (Y3/3 60 0.) ) *15.00 
C HC= DEC/ 57. 2957U 
PA = 3 A/57 . 29573 

C TRIGONOMETRIC COMPUTATIONS START HERE. 

DCL = jCCS (DEC) *DCGS (RA) 
dCM = DCO S (DEC) +DSIN ( R A ) 

DC ii = DSIii (DEC) 

READ 612,T1,T2,T3,NP 


43 



59 

612 

FOR MAT ( 3 F 1 0 . 5,15) 

60 ’ 


PRINT 641, T 1 , T 2 , T 3 

6 1 

641 

FORMAT (5X, ' STARTING TIME ’,1 

6 2 


1 1 AND TIME IS INCREASING') 

o 3 


T 1 = T 1 + (T 2+ (T 3/ 6 0 • ) ) / 60 * 

b 4 


KPP=NP 

65 


1 = 0 

6 6 

73 

IF (T2 .LT. 10.) GO TO 74 

67 


T 2=T 2- 10. 

68 


GO TO 73 

69 

74 

L=T 3/10 • 

7 0 


T X=T2+ (L* 10./60. ) 

71 

75 

IP (T3 .LT. 10.) GO TO 7 

72 


IJ=T3-10. . 

7 3 


GO TO 75 

74 

7 

3=1 .+ (13/10.) 

75 

c* 

SPACECRAFT POSITIONS AND GUA < 

76 

C* 

GMT TIME HI USING THE FUNCTIC 

77 

c* 

Y = A*T *T + B*T + C 

78 


PHI S =C AL ( S , A 2 , E 2 , C2 , D2 ) *C1 

79 


ALT= C AL (S,A3 ,B3,C3,D 3) 

80 


ALAMS=CAL (S, A4,B4,C4,D4) *C1 

8 1 


IF (ALA MS .LT. 0.) ALAMS=6. 

82 


X = 3 . + ( ( E X + (T 3/60.) ) /10.) 

63 


ALAMG=C AL (X, ,V5 , B5 , C5 , D5) *C 1 

84 


1 = 1+ 1 

85 


M=NP-I+ 1 

8 6 


P A= DCO S (ALA MS) 

67 


I- 3 = D C 0 S (PHIS) 

68 


PC=DCOS (ALAMG) 


89 

90 
9 1 

92 

93 

94 
9 5 


OP ARIES ARE CALCULATED AT DIFFERENT 
CAL WHICH USES THE CURVE 


13 


D=DSIN (A LA M3) 

Z = L>3IN (ALA MU) 
DCA=PA*FB*PC-D*PB*E 
DCd=D* PB^PC + P A*PB*E 
DCC=D8IN (PHIS) 

CALL HEIGHT (PHIS , RS A T , RK ) 
X3 = RSAI*iJCA 


9 b 
97 
• 9 8 
99 
100 
10 1 
10 2 
10 3 

104 

105 

10 b 
0 7 

108 

109 

1 10 
1 1 1 

1 1 2 
1 1 3 
,;i4 
\ 1 5 
1 10 

117 

118 


Y S- R SAT* DC 3 
Z S=RS AT + DCC 

HS^DSORT (XF*X3+YS*YS+ZS#ZS) 
ARG=DCA+DCL+CCD*DCM+DCC*DCN 
B ETA = 3* 14 159 2 0 536-DA RCOS (AEG) 

W=DC03 (BETA ) 

XO=XS+ RSAT*W*DCL 
YO = I 3+ESAT* W*DCtt 

IF ( YO. LT. 0, . AND. A KG. GE. 0 • ) BETA- 6. 293186-BETA 
IF (YO.LT.0. . AN D. AUG. LT. 0 - ) BETA= 3* 141593+BETA 
ZO=Z 3 + R S AT *VI *DCN 
ROO^DSylU' (XO* XC + Y 0* YC^ZO*ZO) 

PUIO (ft) = D aTA H (ZO/DSQKT (XO* XO+ YO* YO) ) 

ALAMO (M)=DAUCCS (XO/DSQrtT (XO*XQ+YO*YQ) ) -ALAMG 
L>P-DSIH (BETA) 

D C L S = X 0/ (R3AT*£F) 

DCH8=Y0/ (KSAT*QF) 

BCN8=ZO/ (PSAT*C?) 

CP Ea-DC L* DCL S + BCM *DCMS+ DC N*BC NS 
CALL HEIGHT (P HIO (M) , RP , RR) 
ii 0 - F\CC“ SR 
ZP (rf) - HO* 1OOG00. 

I; (M) =ROC*100000. 



SO (M) = U ;<* 100000. 

PRj.O(M) =rulO(M) *57.246 
A L A i'i L (N } = A LA MO (M) *57.20 5 
2 51= ( 1 . 5 7 0 7 0 0 3 2 6 8- 13 L 1 A ) /C 1 
V3=T3+ { 62.5 15/1 16. 34; 

IF ( I . LT . UP) GO TO 7 
FORMAT (5110.5) 

# 4 * * * * * * * * * 4e >.<1 4 *c 4 * ^ * >fe <. 4 * <c ’fc 4 + £ 4*4* ******** ***************** 4**444 

2 A ST-2 


C* Tills PART TAXES AS INPUT TUI’ FILTER TRANSMISSION, CROSS SECTION 
C* AN J STAR SPECTRUM. THEN A TABLE OF TOTAL INTEGRATED NUMBER DENSITY 
C* TM VS. NOR MALI SEC SIGN A3. F NO , IS COMPUTED 

£****** V ********** V ******** V * * * * £* + **$*****:************£************* * 

TX= AL AMO (NPJ/15.0 
T l=T 1 + TX 

IF ( i’i . LE. 0.) T1=T 1+24.0 

IF ( T 1 . GT. 2 4.0) . T1»T 1-24. C 

L 1 = T 1 > 

T 1= ( T 1 - L T ) *00. 

Lz= 2 1 

11= (T1-L 2) *60 . 

NT A 8 = 106 

C INPUT DATA 

C FILTH DATA IN STEPS CF 20 A. 

READ 3 5 1 , FX, FY ,10 2 
3 51 FORMAT (2F5. 0 , 15) 

KPX = 1.0 0 1+ (FY-I'X) /20 . 

DO 302 1-1, MPX 
302 u J.L ( I) - F X + (I- 1) *20. 

IF (1 02 . NE. 2) GO TC 1260 
M P X '= M P X + 1 3 

HEAD 1 20 1 , .(WLL* (I) , 1=25, MPX) 

1201 FORMAT (8F10. 3) 

1260 R L A i> 300, (FILL (I) ,1 = 1, MPX) 

C CROSS-SECTION DATA 

fi^AD 6G0, (SG (I) ,1=1, MPX) 

C STAR SPECTRUM IN STEPS OF 20A. 

READ 30 3, (FIO (I) ,1=1 ,MPX) 

300 FORMAT (1 1F5. 3; 

30 3 r 0 R MAT (1115.0) 

60 0 FORMAT (6 El 0. 2) 

S U N = 0 . 

DO 101 J= 2 , MPa 

101 S U w = S U N + ( (FI LT (J) *F1C (J) +FILT (J-1) *FIO (J-1) )/2. *(WLL(J) -WLL(J-I) 
C Tills SECTION MAKES A TABLE OF TOTAL INTEGRATED NUMBER 
C DENSITY VS. NORMALISED SIGNAL. 

DO dOO IHC=1 , NTAIi 


READ (3, 12 16) TT 
1216 FORMAT (5E16.6) 

IF ( IC 2 .ME. 2) GO 


1261 


126 1 


1000 


READ (2,121b) (SG (I) , 1=25, MPX) 

C 0 N T I iJ (.J E 
IN (I ilC)"=TT 
S U M = 0 . 

DO 1000 J=2, MPX 

S 0 M = S U M * ( (FILL ( J ) *F lo (J) *DEX? (-SG (J) *TT) +FILT (J-1) #FIO(J-1) * 
1 J E X P (-SG (J-1) *TT) ) / 2m) * ( W L L (J) -*ILL (J-1) ) 

F NR = SUM/ S1JN 

IF (FNR . LT. 1. D-20) r'NR=1, C-20 


45 



F N 0 (I HO) = FN P 
CON TIN UF 




PABT-3 


C* INPUT TO THIS PART IS THE MEASURED SIGNAL- THEN NORMALISED SIGNAL 
C* EOF' IS DSHD TO DETER MI UH TOTAL INTEGRATED NUMBER DENSITY TNO 
C* nY INTERPOLATION IN Till’ TABLE COMPUTED IN PART-2 ***<1, 

<J**+****#V* <****** <«******** *'* *+**+******************************** 

a^AO loci, («(i) ,1 = 1,10) 

l’RIWT 1001, (C(I) ,1 = 1,10) . 

1001 FORMAT ( 1 C Ati ) 

1002 FORMAT (13) 

C OCCULTA I1C I* DATA 

.READ 750, (DA (1) ,1=VNPP) 

760 FORMAT (12F6.0) 

DMIA=D A ( 1) ! 

DO 657 1=1, NPP 

IF (DA (I) -LT- DM1A) DMlA=DA (I) ' 

ti5 7 CONTINUE 

DO H61 1=1, NPP 
661 DA (1 ) =D A (1) -DMIA 
. NMXY=NPP-4 

suc=o. 

DO 761 I=NMXY,NPP 

761 SUC=SUC+CA (I) 

D AO = SUC/ 5 « 

DO 762 1-1, NPP 

762 F OF (I) =CA (I) /DAO 

WRITE (6 ,100 4) 

1004 FORMAT (//'MODIFIED INPUT DATA*/) 

PRINT 1005, (EA(I) ,1=1, NPP) 

100 5 FORMAT (1GF8.C) 

C THIS SECTION CALCULATES TOTAL INTEGRATED NUMBER DENSITY; 

C USING MEASURED SIGNAL AND PREVIOUSLY COMPUTED TABLE. 

DO 611 J = 1 , N T A 3 
L (J) ~D LOG (F NO (J) ) 

XN (J) =DLOG (TN (J) ) 

611 CONTINUE 

609 FORMAT (T 5 , 2D 15. 6) * 

M = 4 

N P = N T A B 
CALL PA BAB (M) 

DO 6 30 L= 1 , NPP 

IF (POP (I) -LE. 0.) FOP (L) = 1* D-5 

DO 7.02 I=1,NTAB ^ 

IF (FOF (L) -GE. FNQ(I) . AND. FQF (L) .LT.FNO (1+1)) GO TO 613 

792 CONTINUE • 

613 J=I 

X = DLOG (BOF (L) ) 

1NO(L)=A{0) *X*X+3(J) *X + C(J) 

TNO (L) -D2XP (TNO (L) ) 

6 1 4 FORMAT (T5, 2D 15. 6) 

630 CONTINUE 

IF (102. NE„ 0) GO TO 2501 

********************************************************** ******** 
C* PART* 3 A * 

C* IN THIS PART, THE OZONE COLUMN DENSITY IS CORRECTED FOR 
C* HER ZB ERG AND RAYLEIGH SCATTERING * 

c * ************************************************* ****************** 


FOF (L) = 1 • D-5 

AND. FOE (L) .LT.FNO (1+1) ) GO TO 613 


792 

613 


2 0 0 D 


3 _> } 


i\ E A J , (C3CJ (X) , I - 1 ,MPX) 

FORMAT ( / c 1 0 . 2 ) 


3 40 


L'O 20 0 5 I = 1 , M P X 

3.4 1 


SG.. u ) --= 3 . 7 5 1 J- 23* ( (1 59 0. »- (I- 1) *20. J * 1. D-4) **(-4.05) 

'♦ *. 

200 5 

CONTINUE 

3 

5 14 

FORMAT ( 1 a , o S 1 5. 4) 

24 4 


W 2 LTH (5,5 1 5) 

2 5 

3 1 5 

FORMAT (//5X, • WAVE J//X/ FILTER* , 10 X, 'STAR IN7. » ,5X, 

2 4 it 


TOi X-ScCTION 1 , IX , *02 X- SECT! ON • , 3.C, 1 HA YLEIGH X-S ECTION • //) 

24 / 


DO 377 1=1, MPa 

2 4 0 


«RI2E (t. ,8 14 J WLL (I) , MLT (I) , F IO ( I ) , SG ( I) , SG0 2 (I) , SGM (I) 

! 4 ) 

5/7 

CONTINUE 

2 ’>0 


READ 2000, ITEST, EPSC 

2s 1 


R E AD 20 0 0 , M Z 2 , I’l O , 0Z Z 

25 2 

2000 

FORMAT ( J! 5 , 2 E 1 2 . 4 ) 

2:> 3 


W R IT E (5, 300 1) ITEST, EVSC , M Z Z , Z ZO , D ZZ 

2 5 4 

300 1 

FORMAT (/T4 , 1 ITEST - ',14, * EPSC = 1 , F7. 4, • MZZ = ',15, 

2 5 5 


1' HZO - ' , i) 1 0 . 3 , 1 DZZ - ' , D 1 0.3/) 

2 5 u 


READ 2009 , (CCM (I) ,1= 1 , HZZ) » 

25 7 

2009 

FORMAT (T2 , 7 E9 . 2) 

2 5 a 


DO 2500 1 1 ■= 1 , N P p 

25 3 


I = NP P* 1 -II 

25 0 


it t=o 

25 1 


K K= (Si? (I) -ZZO) /DZZ+ 1 

2o 2 


I F { K K . L E . 1) KK-1 

2 f> .3 


Cu2 (I) =CCM (KK) * D E X P ( (ZP (I) - (ZZO + DZZ* ( KK — 1 > ) ) * CLOG (CCM (KK+1) 

2 5 4 


1/CCM (KK) ) /DZZ) 

2 5 5 


C M ( I ) = CC2 (I) /0. 2095 

Ho o 


S U N = F A T (0.0, CC2 (NPP) ,CM(NPF)) 

) 7 


FATO=JLOt; (FOP (I) ) +SUN 

c o 5 


IF (FATO.LE. 1. OE-5) 7 ATO= 1.0 £-5 

2o 5 


C 1- i N 0 (I) 

27 0 


FI --FAT (Cl ,C02 (I) ,CM (I) ) 

271 


C 2= C 'i *0.7 5 

272 


IF (C 2 . LE. 1. 0L10) C2= 1 . 0D10 

27 j 

2002 

r 2 - FAT (C 2 , CO 2 (I) ,CM (I) ) 

2 1 4 


IF (DABS (r 2-F 1 ) . LE. 1 . OS- 10) GO TO 200 4 

27 b 


CJ=C1+ ( F A TO - FI J * (C2-C1) / (F2-F1) 

27 (; 


IF (DABS (C3) .GE. 1. 0EJ5) GC TO 2004 

27 7 


I TT= I TT+ 1 

. 27 y 


IF (DABS { (C3-C2) /C2) .LE.EPSC) GO TO 2004 

2 7 9 


IF ( (iTT-ITEST) .GT. 1) GO TC 2004 

2b 0 


IF (DADS ( (F 1-FA TO) / (F2-FAT0) ) • LE. 1.0) GO TO 2003 

2 3 1 


F 1 = F 2 

25 2 


0 1=02 


253 
2 a n 
2 a 5 
\n a 
2(7 
2a a 
2 a 'j 
20 0 
2 7 i 
25 2 
j 

>?4 

20 5 
29 5 
20 7 


2003 C 2 = C J 

GO 20 2002 

2004 CONTINUE 


*!'!0 (1) =C 3 

IF (C3.LT.0-0J T NO (I) =0. 0 
2 500 CON TIN U E 

C ***** * ****** ********** ** ********** ** ** *** *** ********** *************** 

C* PART- 4 

C* 


C* 

C* 

C* 

c* 

c* 


7HIL PACT TITS A PAR A 30 LA THROUGH TANGENT RAY HEIGHT ZP COMPUTED 
IN PART-1 VS. TOTAL INTEGRATED NUMB ER DENSITY T NO COMPUTED IN p ART 
~3. THEN USING LEAST SQUARE COEFFICIENTS, INVERSION IS CARRIED ON 
To COMPUTE THE NUMBER DENSITY. THIS NUflB&K DENSITY IS, THEN, USED 
IN THE CALCULATION OF MOLECULAR OXYGEN TEMPERATURE. 


— - ’ - — w W W* ft * V U « 

c* ********************************* j^***^^ ******************** ** * 
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24 8 

2501 

29 'j 
30 0 


30 1 

1211 

302 


’30 3 
30 4 
305 

C 

30 6 
307 

21 

30 8 


109 

123 

310 
3 1 1 
312 
31 3 

12 

314 

315 


316 

317 
316 
319 


320 

321 

4 

322 
32 3 
324 

601 

325 

326 

327 

23 

326 


329 

330 

331 

332 

333 

334 
33 5 
33c 


337 

336 

339 

150 

34 0 

577 

341 

342 

22 

343 

978 

344 
_ 34 5 


34 6 


347 

348 
345 

350 

351 

352 
■ 353 

'354 

9 11 

355 
3 Db 
357 



TIME* ,16, Mi OUR,' ,16, 'MINUTE, 


N 1'Mlt-L 

DO 1211 1-1 , fU’P , 

2(LJ=Z?(I) 

XN ( I) = 4 I MC (I) 

Tais*SKCTIO«“oE'f AIKS XIU COEFFICIENTS FOB LEAST SQUARE PARABOLA 
K-a 

WRITE (6,21) M 

FORMAT (' 1* , V SMOOTHING WITH M=',I3//) 

WRITE (6,123) 

CALL PA 11 A 13 ( IT } ri . ipirtmel//( 

FORMAT (' IGNORE TEMP. IN OZONE RESULTS //} 

FORM AT (5 X , 6 E 1 2. 4) 

N 1 = N P- 1 
DO 4 J- 1 , N 1 
DEN (J) - 0.0 
DO 5 1=0, N1 

DEN W=DEN (J) ♦ U-'tDSaRX **2-B (J) **2) -DSQB.T (R CI)**2-B (j)**2) 

1 + (it ( I) -i.*A(I.)*RO (I) )* 0LCG(IR{I+1)+1)SQRT(R(I + 1) **2-B (J) **2) ) / 

2 (U(i) +DSCFT(R (I) **2-E (J) **2) ) ) 

DKN (J)=-DEN (J)/PI ■ 

CONTINUE 

W RITE (b ,60 1 ) LI , L2 ,T 1 
FORMAT (5X , * LOCAL STARTING 
1F7.3, ' SEC CNDS. ') 

W FORHA'T # (//^0X r ' HEIGHT * ,9 X, 1 NUMBER * ,9 X,* TEMP. * , 10X, 1 LAT * * * 

1 ' LONG. • ,10X, 'COLUMN' ,10X, • NOBM ALIZED' , 1 0X , ' 02 COLUMN DEN^ITI , 
2,1 IX, 'IN CM.',8X,'DEN5ITZ',aX,'IN DEG. K. • , bX , IN DliG. ,8X, 

3<'iM DEG. ■, BX, -DENSITf , 9 X , ' INI ENSIT* • , 1 IX , 'DENSITY',/) 

a i = n i - 1 

DO 977 J= 1 , N 1 

T (J) = 0. 0 

DO 150 I = J , N 1 

A2-DEN (I*1)/(a(I+l) **2) 

A 1 = D F. N (I) / (3 U)**2) 

A 1 = D ABS l A 1) 

T (J) =T (J ) + ( A 2** A 1) /DLGG (A 2/A 1 ) * (Z (1+ 1 ) - Z (I) ) 

1 (J) = (3.8033D-4)*T(J)/DEN (J) 

/ WHITE (6 \ '12) ( z\j) ,DEii [J) ,T (J) ,PHlO<J) , ALAMO (J) ,XN (J) ,FOF(J) ,C02 
FORMAT (4X, 6 El 5. 6, El 8. 6, E 19.6) 

STOP , 

DO 911 J = 1 , N 1 
UO(J)-DEN (J) 

UZ (J) =Z (J) . .. 

UT (Jj =T (J) 

IF ( U D ( J } 

(J Y ill H=7 . 0 
U 7 2X-0. 5 
U XMIN=40. 0 E 05 
UDX = 20 . 0 E 05 

IF (102 . EQ. 0) GO TO 1531 

CALL PLTOFS (100.0,100.0, U-XHIN , UDX , 1 . 5 , 1 . 0) 

UPX=1G0. 

C\LL PA/IS M . 5, 1. , 1 6HTHHP. IN DEG. K. ,-16,10. ,0., UPX ,UPDX,-1. ) 
CALL mWtl.I.l/lSUALTIIUDE IN CM. , 1 5 .8 . ,90. , UXMIN ,0DX. -0. 5) 


LF,. 0.0) UD (J) =1.02 01 
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5 5 3 


v. A.. 1. P‘.»R 1 2 \ ) • 5 , 1 . 0 , 0 . 5 , 0 . 9 , 2 0 , 1b) 

5 > 9 


C A L n I L I, N L ( U I ( 1 ) , U Z ( 1 ) , N ! , 1,1, 1,2, 1 ) 

8b 0 


CALL P B \ 8 B ( 1 • 8 , 9 . 5 , 0 » 1,0(1) ,0*0, 8 w) 

i ;> 1 


l F k M .HQ. 4) GO TO 13 80 

w »2 


CALL PGYMB ^ 1 . 5,9. 25,-0. T, 3lW'l- il ,0. 0,3) 

'% ; 


G v * TO 1 ‘ > 3 2 

_c> 4 

1 i 0 

C At. L P b Y N li £ 1 . 5 , 9 . 2 5 , - 0 . 1 , 3 11 K-- 4 , 0. 0 , 3 ) 

5 » > ! > 

\ 83 2 

CA LL i'l. VEND 

.in (V 

153 1 

CONTINUE 

Jo / 


CALL PLiOFS ( U Y MI N , U Y La , U a cl IN , U DX , 1.5, 1. 0) 

it, 3 


CALL PLGAXS ( 1. 5, 1. 0, 14 tl NUMBER DEN SI l’ Y 14 ,10, 0 , 0 . Q , U Y Hi H , -U Y 3 X) 

3 o') 


CALL L'AXiS (1 . 5, 1. 0, 1 5HALVITUDE LiJ CM. , 15, 8. 0, 90. 0, UXHIN ,UDX,-0. 

370 


CA,.L PLGGRD ( 1. 5, 1 .0, 1 .0/UYEX, 10-0,8.0,0.0) 

3 7 1 


C a L L PG 810(1.5,1.0,10.0,0.5,1,15) 

3 7 2 


CALL PIT LCS (2) 

J 7 3 


CALL BLINK (UJ (1) , UZ ( 1) , N 1 , 1 , 1 , 1 , 2 , 1 ) 

•37 4 


CALL FIT RFC 

37 5 


CALL PS YMB (1 . 5,9. 5,-0. 1 , Q (1) ,0 .0, 80) 

37 o 


jlF ( M .EC., 4) GO TO 912 

37 7 


CALL PSY Ho (1.5,5.25, -0.1 , 3 H M - 8 , 0 , 0,3) 

37 d 


GU TO 913 

37 9 

912 

CALL PSYKB ( 1 . 5,5., 25,-0. 1 ,3Hft-4 ,0. 0,3) 

3 6 0 

9 1 3 

CALL P .LI END 

38 1 

6 

C O N T I K U L 

3B 2 


CALL EXIT 

38 3 


END 


(J ********* ******* ************************************** ** ************4 

J-R-5 C* PART-* 5 


it- t: 



i t* 8 
39 0 
>9 1 
392 

39 3 
39 4 
39 5 
5 V o 

39 7 
j 9 o 
.39 3 
4 1 ' 0 

40 1 
4 0 2 
40 3 
40 4 
40 5 
40 c 
4 0 7 

40 5 
4 1 0 
4 1 1 
4 1 2 


C* 

c* SUBROUTINE SLOT TAKES 5 EQUIDISTANT VALUES AND FITS A LEAST SQUARE 
C* PARABOLA Y=F*X*X + y*X + k . COEFFICIENTS P,Q AND H ARE RETURNED 

(J * * * * >< * * * * * * V * * ** * ***V****'*«: *********** ********************** ***** * ** ' 

SUBROUTINE SET (V 1 , Y2 ,Y J, 4’ 4 , Y5 ,P,Q , R,S) 

IMPLICIT' REALMS { A-il, O-Z) 

COMMON l (200) , a N (200) , A (200) ,D (200) ,C (200) , ALT , WLL ( 1 00) ,FILT(10( 
1FIO(100) , SC. (3 00) , SGO 2 ( 1 00 ) ,SGK(100) ,CCS (200) ,NP,KPX 

?■ ( (2. + Y 1) - Y2- £2. *Y J) - Y4 + (2.* if 5) ) /14, 

QM (- 146. *Y1) + (4D. *Y2) + < 120. J) + (74.*Y4) - <92.*Y5) ) /140. 

R = ( ( Lb, *Y1) - (5(>.*Y J) - (4 2.*Y4) + (4 2.*Y5) ) /70. 

RETURN ' 
ii N D 

ij ******************************************** ************************ 

C* PART- 6 

C * 

C* SUBROUTINE HEIGHT CALCULATES TilE RADIUS ‘OF EARTH RR, AND DISTANT 
C* R P OF A PCI NT i'POri CENTER Of EARTH; WHEN LATTITUDE PHI, AND ALTII* 
C * ALT, OF THE POINT ARE KNOWN. 

^ ****** **************** ** ********************** ******************** ** 
SUBROUTINE. HEIGHT (Pli.l, HP, SB) 

IE PL 1 SIT K F. A L * 8 ( A - H , C - Z ) 

CC.1SGN Z (200) ,XS(200) , A (200) , 3 (200) ,C (200) , ALT , WLL (1 00) , FILT (1Q f 
IF IO (100 ) ,3G (100) , SCOP (100; ,St3M (100) ,CCM (200) , NP,MPX 
REAL*8 KAPPA 


Vi 4 

415 EC--. 0818202 

416 K=6378.387 

417 uENOFi = DSQRT ( 1.- (EC *EC*DSI N (PHI) + DSIN (PHI) ) ) 
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ACO$=DcCS (PHI) * K/ D E N C M 

as ln = r* { i tc^fcc) *us in crux ) /dlnom 

K tt - i) S C K 1 l H C 0 S *RCO S + R FIN* BSI N) 

PHIP = DARSJ'N (K3IN/KR) 

KAPPA=3. 1415926536 

is ?- DSL'KT (KR *NK+ALT*AIT-2. *EL*ALT*DCOS (KAPPA) ) 

RETURN 

END 

****** ** * * * * * * * * * * * * * v **************** ******** ******************** ** * ♦ 

C* PART- 7 

C * 

C * SUBROUTINE PALAU 1 ; ITS TilE LEAST SOU ARE PARABOLA THROUGH M POINTS 
C* AT A TILE IN THE CURVE OP Z VS. XN. THE PROCESS IS CONTINUED FOR 
C* NP POINTS AND COEFFICIENTS A, 13 AND C ARE RETURNED TO MAIN PROGRAM 
C* NOTE THAT XN=A*Z*Z + 6*Z + C 

(J**** ******** * fr* ******************************** ************************ 

SU5R0UTI NE' PAR AH (M) 

IMPLICIT RZAL*8 (A-H, C-Z) 

COMMON Z (200) ,XN (200) ,A (200) ,B (200) ,C (20C) , ALT f WLL (100) ,FILT(100) 
1FI0( 100) , SG .( 100) , SG0 2 (100) , SG M (100) ,CCM (200) , NP,MPX 
L = &/ 2 ’ 

DO 2 1=1, HP 
IZ-NP+1-I. 

K = 1 

IF (I.GT. L. AND. I. LT.IS) K=.I-L 

IF (I.GE.IZ) K-NP-M 

25=0. 

Z 213 = 0 • 

235=0* 

245=0. 

XN5-0V 
Z N B = 0 . 

Z2NB-0.. 

K 2=K + M 
DO 3 J = K,I<2 
ZB=ZB + Z(J) 

Z2B=22B+Z (J) **2 * 

Z3b=Z3D + Z (J) ** 3 

Z4B=Z4B+Z (J) * *4 ' * 

2 H 13 = Z N 5+2 (J) *XN (J) 

2 2N B=Z 2N 3 + XN (J) * (Z ( J) **2) 

3 X N8= XN9+XN (J) 

J = I 
M 1 = M 
rt=ii + I 

A (J) = ( ( M * Z N B - X N B * Z B ) * ( M* Z3 B-Z B* Z 213) - ( M*Z 2N3-XNB*Z 2B) * (M*Z2B-ZB**; 
1) / ( (M*23L 3 ~Z23*ZB) * ( M*Z3 B-ZB* Z2B) - ( M* Z4B-Z 2B**2) * ( M*Z2 B-ZB** 2) ) . 
C=A (J) , 

*3 (J) '* ( (M *ZM 3- ZB* XNB) - ( h* 23B-Z 3*Z2B) *D) / (M*Z2B-ZB*23) 

Z = 5 (J) 

C (J) = (XNE-D*Z2B-2*ZB)/M 
M = 1 1 

2 CONTINUE . 

• RETURN 
END 

FUNCTION CAL (X,A # B,C # D) 

IMPLICIT REAL*8 (A-H r O-Z) 

CAL* (A*X**2) + (B*X) *C 
RETURN 
END 
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Part l: 


Line numbers 1-125 

In this part, the spacecraft position values are entered for use in 
interpolation. There are five such values for longitude, latitude, and altitude 
furnished by GSFC several weeks after the occultation. The first of these 
five values corresponds to the time (given at a 10-second interval) which is 
less than but nearest to the starting time. The Greenwich hour angle of Aries 
is then read in for five times spaced ten minutes apart and centered upon the 
time less than but nearest to the starting time. Star position information is 

i . 

also entered, followed by the starting time and the number of data points. The 
position of the data point is computed from the above inputs, using subroutines 
HEIGHT and SET and function CAL. Program comments for Part I give 
additional details. 

Notation 

DEC = declination of star, in radians 
RA = right ascension of star, in radians 
NP, NPP - number of data points in scan 

PHIS = latitude of subsatellite point in radians 
ALT = altitude of spacecraft, in cm 
ALAMS = longitude of subsatellite point, in radians 
ALAMG - Greenwich hour angle of Aries, in radians 
RSAT, RP = distance of subsatellite point from Earth's center, in cm 
ROO, R = distance of data point from Earth's center, in cm 
HO, ZP = height of data point, in cm 

RO, RR = radius of Earth along a line to the data point, in cm 
PHIO = latitude of data poins, in degrees 
ALAMO = longitude of data point, in degrees 
Tl, T2, T3, TX, X, S = temporal variables 
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Part 2: 


Line numbers 126-180 

In this part, the filter transmission values, the appropriate absorp- 
tion cross section, and the star spectrum are entered for the pertinent wave 
band at a 20 A interval. A normalization factor for the signal intensity is com-' 
puted and a table of normalized signal intensity v, tangential column density 
(integrated number density) is formed. The latter values are already stored 
in file CROSSDAT. 


Notation 


'' ■, 

NTAB 

= 

number of tabular points Q 

FX 

= 

beginning wavelength of band^ in A 

FY 

= 

final wavelength of band, in A 

102 

= 

switch for gaseous species and absorption cross section 
(liozone, 2 : 0 2 , 3 : 0 2 with Hudson's cross section) 

MPX 

= 

number of wavelength intervals in band 

WLL 


wavelength 

FILT 

' = 

normalized value of filter transmission 

SG 

= 

2 

absorption cross section, in cm 

FIO 

= 

spectral star intensity, in digital voltmeter 
counts/time interval 

SUN 

= . 

normalizing factor, in counts/time interval 

TT, TN 

= 

2 

tabular tangential column density, in cm 

SUM 

= 

filtered signal intensity integrated over wavelength, 
in counts/time interval 

FNR, FNO 

s 

tabular, normalized signal intensity corresponding to TT 


Part 3: Line numbers 181-233 

In this part, a comment card is read which identifies the occultation 
data. The data are then entered, the background noise (minimum signal) is 
removed, and the data are normalized. Using (m+1) smoothing points, a loga- 
rithmic interpolation is performed in subroutine PARAB to get the least squares 
coefficients for the column densities corresponding to the "modified input data. M 
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Notation 


DA 

DMIA 
FOF 
TNO 
A, B, C 
M 


= star intensity observed during occultation, 
in counts/time interval 

• i 

= minimum value of DA 

= normalized, measured signal intensity 

. -2 

= column density along stellar ray path m cm 
= least squares coefficients for a parabolic fit 
= one less than the number of smoothing points taken 
for a least squares fit 


Part 3 A : Line numbers 234-288 (for ozone only) 

In this part, a correction to the ozone column density is made for 
the removal of Herzberg and Rayleigh scattering. To accomplish this, a 
functional inversion of F(r) is necessary, where 


F(r) = 


J Q TU)i 0 (M exp[- <To 2 (X)N o 2 (r| -I'M (X)N M (r) - ^o 3 (X)N o 3 < r) ] dX 


I T(\)I (\)dX 

/ 0 


o 


F 

T 

I 

o 

v i 


N. 

l 


normalized signal intensity 
filter transmission 
spectral stellar intensity 

th 

cross-sectional value of i species or of molecular 

2 

scattering (M), in cm j 

- th 

tangential column density of i species 

(Nm = total for all species), in cm ^ 


The desired quantity is N 0 (r). To carry out this operation, Kj Q J,\) is read 
in for the ozone band being considered e. g. , 1800-3000 A, with the Herzberg 


continuum being included. \f ' (\) is then computed for the same wavelength 
region from Rayleigh's formula, 




3. 981 x 10 _28 (lCf 4 \f 4 ‘ 05 om 2 


(\ in A) 
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Some representative values of 0 2 tangential column density are read in at a 

5 -km interval from 10 km to 310 km, and exponential interpolation is employed 

to find this quantity at each data point height. It is then converted to total 

tangential column density for all species and the final value of is reached 

through Newton's inerative method. In this method, an initial guess of N M is 

made for the value of N^ o) and function "FAT" produces a signal F x which is 

depleted by passage through an atmosphere with tangential column densities N^, 

N and N. . and through the filter. The second guess reduces nJ” by 

one quarter, and if the resulting signal F g is close enough to F x we have the 

N which is sought. Otherwise, a new is computed by the formula 
Oq 3 


N i = N j ' 2 + (F -F. i XN ^" 1 - N->' 2 )/(F, 
o, o, O 3-1 °3 °3 1 



where F is the sum of the "FAT" functional evaluation for the case where 
o , 

N =0 (actually a new normalizing factor) and the log of the signal corres- 

O3 ^ j — ^ j-1 | 

ponding to N m . A second tolerance test is provided by / (N^- )/ | 

< £ and a limit is placed on the number of iterations. The data points are 
taken from top to bottom and certain safeguards are included for computational 
instability encountered at the lowest levels. 


Notation 

SG02 

SGM 

1TEST 

ITT 

EPSC 

MZZ 

zzo 

DZZ 


= 0 2 cross section ( in cm 2 

s . Rayleigh scattering cross section in crn 

= maximum number of iterations 
= iteration counter 
= tolerance 

= number of values in interpolation table for 0 2 column density 
= base height of interpolation table, in cm 
s vertical interval in table, in cm 
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Notation (cont'd) 


tabular 0 9 column density, in cm 


-2 


interpolated 0 9 column density (NC^), in cm 


total column density (N^), in cm 


-2 


COM 

co 9 

CM 

SUN = new normalizing factor 

FATO(F q ), CKN^' 2 ), C2(N 0 ^" 1 ), C3(N 03 j ), Fl(F j . 1 ), F2(F ) 
= parameters in Newton’s method 


-2 


Part 4: Line numbers 289-383 

Using subroutine PARAB, the coefficients are obtained for fitting a curve 
of total tangential column density against tangent ray height, and the mathe- 
matical inversion is performed to get number density. Please refer to the 
notes on the error analysis for details on this procedure as well as the 
computation of temperature. 

i 

A STOP statement (line 342) negates CALCOMP plotting. When STOP is 
removed, the plotter produces number density graphs for two smoothing values. 
When molecular oxygen is being investigated, it draws two additional graphs 
for temperature. 

Notation 

Z = height, in cm 

_ o 

XN = tangential column density, in cm 

- 3 

DEN= number density, in cm 

T - molecular temperature, in °K 

Parts 5-7: Line numbers 384-499 

Subroutine SET returns the coefficients after fitting a least squares 
parabola through five equidistant points, whereas subroutine PARAB fits 
any number of points which need not be equidistant. Additional information 
is available in the program comments and in notes by Mr. A. Shah. 
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Part 8 


Line 

500 

501 

502 

503 


504 

505 

506 

507 
508-513 
514-524 
525-530 
531-539 
540-541 


Line numbers 500-541 

Sample data are listed for an ozone run. They are: 

Longitude of spacecraft at a 10-sec interval, beginning with the 
value nearest to, but less than, the starting time of the scan, in g 

Latitude of spacecraft, .ditto, in deg 
Altitude of spacecraft, ditto, in km 

Greenwich hour angle of Aries at a 10-min interval, centered around 
the value which is nearest to, but less than, the starting time of 

scan, in deg 

Declination of star in deg, min, sec 

Right ascension of star in hr, min, sec 

Starting time in hr, min, sec and number of data points 

Beginning and ending wavelengths, in A, and species identifier 

S2F5 filter 

2 

Ozone cross section, in cm 

Stellar spectral values of intensity over the same wavelength intervals 
Molecular oxygen cross section over the same wavelength intervals 
Control values for corrections in Part 3A 
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Appendix nib. 

This Appendix explains the procedure of finding a random number 
from a Poisson distribution as an additive scattering term. The procedure 
is made necessary by the fact that a Gaussian error distribution with a 
variance equal to the signal intensity causes negative values to occur at 
times when the scattered signal is superposed upon a basic signal as low as 
5 counts /&t. 

The Poisson distribution gives the probability of the number of occur- 
ences of an event within a given time interval, e. g. , 1 sec, and it is quite 
pertinent to the present case. A Poisson distribution, as described by the 
expression 

p(n) = exp [-t] 

has the values shown in figure when t = 10, 1, and 0.1. n must be an 
integer, but t, the mean signal intensity, is a real number. Since we seek 
a random number from a Poisson distribution to be added to the signal as 
random noise, an inversion of the formula is desired. In lieu of this 
difficult inversion, we can search for a random number, y, O^y-^1, which 

will permit an n to be generated such that p(n)<y< p(n+l). When the inequal- 

* 

ity succeeds, n is accepted as the simulated signal intensity. 

The procedure is employed successfully when the signal is between 10 ^ 
and 10 counts /At. At lower values of the basic signal, the computer search 
for Poisson random variables becomes too time-consuming so the signal is 
then left intact. 
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Appendix IIIc. 

This Appendix uses a simple flow chart and copies of error 
analysis programs with notes to give the details of the numerical error 
analysis for each of the species studied. Both single wavelength and 
broad waveband cases are considered. 

The first flow chart applies to molecular oxygen, which is 
assumed to have a monotonically decreasing number density with increas- 
ing height. Following this figure is a program, "THREE 11 , which handles 
the braodband filter case. Explanatory notes may be found at the end of 
the program. 

The second flow chart applies to ozone, which is assumed to 
have a Gaussian -peaked feature superimposed upon a monotonically 
decreasing number density with increasing height. Following this figure 
is a program "OZNUM", which applies to the broadband filter. Explana- 
tory notes for this program terminate the Appendix. 
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F ( 1 ) -;v- ( - { . 9 , 1 1 * V’i + 1 U 3 8 78 ) ) * ( WL 0**5 ) / t ( WL <1 )*1 .L >~3 ) **5 ) * ( DEX 

l \ ... ^.i/ {v:i ..<’* 1 F ) ) -1 . > / ( ill XP (1 . A 3 8 / ( W L ( I ) *1.0-8* TF ) ) - 1. ) * ( 1 . 0- 7 ) 

not i )=!• Ili( i }*mn 
301 F'T n^’UF 

c 

c C ' <« 1 ' i’ l J 1 1 Nl I ! ■ M A 1. [ / A r I UN I At T D l< 

C . 

SI if 1 . “A 1 . 

!)•( to just* 

>0? b i JN --- SUN «• ( { i 1 LI ( J ) *F m £ J ) f F I L T ( J- 1 )*FI0< J- LH/2. *{ WL ( J)-WL(J-l) ) ) 

WRI I f ( h , 303 ) SUN 
808 i [)9WAT( F ID. A ) 

;i,i=n 0 .n 1 ,/ ( N n — 1 J 

0 

r H, O'* A TABU; OF COLUMN DENSITY! TT) VERSUS SIGNAL INTFNSITY 

C ( SU v\ ) EASE I) I J pn N A N T NT L G R A T I ON 0 F TH E FORMULA FOR T T WHICH 

C ASSUMES SPHERICAL STRATIFICATION AMD AN ISOTHERMAL ATMOSPHERE 

C 

00 800 I HO = 1 , HP 
SUM^G. 

!*’.' ! - R E + HU+ f T HO- 1 ) * HD 
ZO-RO-RF 

IM I HO. F tJ . 182) XNi : = 5 * i.H 3 
IF ( I HO . HO. 187) XNO = 3 • I) 1 3 
T r ( IMO. EQ. 1 q/ U XNO1.013 
no 701 1-13 8 f 210 

X K = { 1-184 )/2. 

701 I F f IHH. EO. I) X NO = 1 .013/1 . 0 1 * * X K 

T T = X NO*ni: X P ( - { z O- 8C . 0 8 ) t SH 1 *D SQR T ( 1.570 79* SH*RO ) 

1 * ( 1 . + 8 .*SH/ ( 8 . * R 0 } - 1 *3 . / 3 2 . * ( $ H / R 3 ) * *2 + 3 1 5 . / 2 1 6 . * ( SH/R0)**3) 

IF! IHJ. CO. 211) TT = 1. D-1C 
T i 1UDI=TT 
. Oil i rv,o J-?_ ♦ 3? 

1CO0 s U- v T SUN + ( { F I L H J ) *F 10! J 1 *DC XP(-SG(J)*TT)+F I LT !,J-1 ) *F I 0! J -1 ) * 

1 DFXPf -SO! J-I ) *TT ) )/2. ) *( WL< J )-WL ( J- 1 ) ) 

F NO! I HU = SUM/ SUN 
FSUH! I!G ) = SU M 
A 0 0 C ONT 1 MJ E 

1 F ( NN . f: 0. 2 ) GO TO 30? 

C 

C LOOP FUR NT RANDOM ERRORS IN SIGNAL INTENSITY AND VARIOUS 

C. DATA PUI Ml SPAC 1 N GS 

C 

DO R .11 l- 1 , J I / MAX 
Di: <>0 2 I H= I , f !!‘P ,NH 1 

FUM IH) = FNU ( IK) 

W ! M F L ( 0 t 5 0 1 ) IHtlUMIH) t FNf!! I H) » F 0 F ( IH) 

■SOI FORMAT tT 2, IS, 80 1 5 . 6 ) 

60? CONTINUE 

GO TO 306 ■ 

80? CONTINUE' 

I f ( MM .NE .MMM) GO TG 803 
C 

C LOU P FOR ACCRUING ERRORS IM 0? DENSITY RESULTING FROM 

C RANDOMLY SCATTERED SIGNAL INTENSITY, FOR MANY CASFS, 

C WITH A SINGLE DATA POINT SPACING 

C 

DO 18 JIX-I.MM 
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117 

• 503 

CUNT INUF 

i 3 3 

C- 


11'. 

C 

si.-;i up for gauss ssp at 

' ? 9 

(■r «■. * 

0 

I X= 203+ l 00* ( J I x-l ) 

1 ? 2 


I F { MM . N H . M ) I X - 2 C 1 

\ j > ; 


J\=TX 

i 4 


A' ! = U. 

12 5 


P(1) = 0., • 

1 2 6 


IF ( MM . f 0 • MMM ) GU T G 80 4 

" 27 

c 


' ? c> 

c 

L f 0 p F 0 V. P A N 0 0 'A L P P C R S 

L?° 

c 

DATA POINT S PACINGS 
1 

130 

c 


131 


DO 26 J IY= 1 1 J I YMAX 

132 

8 04 

DO 310 lM=lfNPPfNH 

133 


SUM.= FSUIM IH) 

134 


IF l SUM. LT. 10. )G0 TO 9 

1 3 5 


Sl.) = 0 SOP T ( SUN ) 

136 


C ALL GAUSS ( I X * S l) ♦ AM t V ) 


137- 
1 3* 
1 39 
1*0 
l '4 1 

l'<? 

143 
1 4 4 
1 4 4 

1 * 1 7 

1 4 3 
]'m O 

1 4 0 
131 
1 4 3 
1 33 
1 54 
i 3 ^ 
1 56 
157 
1 5 3 
15° 
i 67 
■ 17 1 
16? 
/■ "• 3 
> /. 


503 I F ( FOF ( I H ) • 0 [* * 1 • )FUF(IH)=1 .-1.0-5 

I r- i FOF ( I H 1 . LT • f NO l l ) 1F0H IH)=FN0( 1) + 1.0-12 

M- 100 

. GLJ TO 809 

o i f { SUM .LT . . 1 ) FO F { IH) = FNO( IH ) 

IF ( SUM. LT. . 1 ) N*0 • i 

IF(SUM.LT..nGO TH 809 

■ DFTFRMINF A RANDOM NUMBER FROM a POISSON DISTRIBUTION 

W ML: N • 1 < S I GN AL I N T EN S l T Y < 10 • 

N‘1AX= SUM +5. 

|- =DLXP{ - SUM 1 
P ( 2 ) = $UM* F 

2 fiSUM = H ( 2 ) 

C AL L RAN DU ( JX f JY , Y F LJ 
DO 1 N- 2 » N-M A X 
I F( N. HO. 2.) GO TO 1 0 
U“ SUM** ( N-l ) * F . 
n=i. 

Il=N-2 ' 

-■> J r l » i I 

3 • C= (N-J ) * D 

GSUM=.CSUM+U/ 0 ( 

P (M) =GSUM ' 

10 IFOM.N) .GE.YFL .AND.P(N-l) . LE.YFLJGG TO 12 

IF (MM.NE.MMMJWR IT E(6, 91 )N, P (N-l ) , P (N) i'Y'FL 


165 

91 

FORMAT ( I. 5 , 3 F 1 2 . 5) 

166 

1 

CONT INUF 

167 

13 

J X ~ J Y 

1 6 H 


GO TO 2 

169 

12 

JX = J Y 

170 


FOF ( IH) =(N-2 ) /SUN 

1 71 


IF (FOF( IH) •LT,i.D-lO)FOF( |H) 

^72 

.« C. 9 

I f { MM , NC ♦ M H M R I T L ( 6 » 3 0 5 ) IF 

^73 

30 5 

FORMAT ( I 5 f 2015.O, l5 t Fl G. 5) 

174 

3 1 9 

CONTINUE 

175 

3 0 6 

HD P = H D P * 2 . 

1 7o 


T r ( MM . HO .MMM ) Hf)P = 2.05 


£4 


177 

17* 

17*5 

' F 0 

W*i 

187 
1 F 3 
184 
1 84 
l 6 4 
- N 3 7 

v i 88 

r)0 
I'U 
1.0? 
19 8 
1 94 
1 9 4 
1°6 
1*7 

1 <5 8 
1°9 

2 0 C 
201 
212 

203 

204 
\-!3 

W--6 
2 07 
2( 3 
20° 
2 10 
?11 
212 
2 13 
214 
2 1 ’3 


2 16 


217 

2 13 

219 

220 
221 
2 27 

/' ;> 3 
".. r 2 4 
! 2 4 

228 
227 
2? H 
2 7 9 
7 ? 

: i 



2 24 
239 
22 6 


NP|!= ( NPH-1 > /2+E 
l F ( J I Y. EG* 1) NPH = 181 
I F { J I Z.E0.1)NPH*181 
I F t MM * tO • MMM ) NPH=9 1 

C COMPUTE 02 NUMBER DENSITY FROM AN EXPONENTIAL MODEL 

C ATMOSPHERE FOR LATER COMPARISON WITH RETRIEVED VALUES 

C 

00 80 1 11-10*1 * NPH 
ROs-KE + HO-M IHC)-1)*HDP 
CALL ATMOD(RD) 

‘ ONIJMI ) ( 1 I I 0 ) = X P ( 2 1 > 

RHTO(IHO)*RO 
ZKM( I HO ) = RHT 0 ( I HO 1-RE 
801 CONTINUE 
C 

C INTERPOLATE IN THE TABLE TO FIND COLUMN DENS I TY < TNO ) FROM 

L A SIGNAL INTENSITY(FOF) WHICH IS A SIMULATION OF THE 

C MEASURED SIGNAL * 

C 

DO 902 J = 1 »NP 

Z( J ) ■= 0 L 0 G ( F N 0 ( J ) 1 
XNl J1 =DLOG ITWWI J1 ) 

902 CONTINUE 
M = 8 

1 F < NN.EQ* I) M=2 
CALL PARAB(M) 

DO 630 L=1 *NPP,NH 
on 612 1*1, NP 

I F( EOF { L 1 .GE .FNO ( I) .AND. FOF( L ) *LT ,FNO( 1+ 1 > )G0 TO 613 

612 CONTINUE 

613 J = I 

X s PLuG ( F OF ( L ) 1 

TNU < L ) = A ( J) * X* X + 0 ( J ) * X + C ( J 1 

TND ( L ) = DE XP ( TNO ( L ) ) 

630 CONTINUE 
C 

C INVERT COLUMN DENSITY(TNO) TO GET 02 NUMBER DENSITY(DEN) 

C 

DO 700 J * 1 , N P H 
7 ( J ) * ZKM { J ) 

K- J< £ NH— ( NH — 1-J 
700 XN(J)=TNQ(K) 

M=8 

I F ( NN • EO* I 1 M = 2 
CALL PAP. AD(M ) 

N L = N P H - 1 
DU 4 J= ItNI 
DENI J ) “0 • 

DO S I - J , N 1 
0*2. *A{ I 1 

5. OEN( J)=DEN< J ) +U* ( USQRT ( RHTO( 1 + 1 ) **2-RHT0< J ) **2) -DSQR1 ( RH TO ( 1 )** 

1RH T0( J) **2)1 + (B (I>-2 .*A(I )*RE K0L0G1 ( RHTOU + 1 1+DSORT (RHTOt 1 + 1) 
. 2-RHTO(J)**2) ) / 1 RHTO ( I ) + DSQP. T( HHTO( I 1 ** 2-RHTO ( J ) ** 2 ) ) ) 

DEN(J)=-DENtJJ/PI«2, 

C 

C COMPUTE ERROR ANC ACCRUE ERROR SUMS 

C 

DERR ( J ) = ( DEN ( J ) -fJNUMD ( J ) ) /ONUMD ( J 1 
DERRSQi J >=UERRSGU)+DERR ( J)*OERR ( J ) 
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7 .5 7 


f '' ** 

7 ^ 

4 

."<•0 

W ' 1 

•f 4 : 

.88 

2 A 3 
? A /- 

7 A 

2 A 8 

7 5 

7 6 

2 c> 

'? A 7 

2 3 

. * , 8 

C 

7 A 9 

c 

p 5 

c 

281 
2 8? 
? r '3 
2 *3 A 


2*5 


2 86 


2 87 
? 8 R 


2 l >9 

7 

26 0 


?6i 


2 62 


2 6 3 

22 

2 64 

6 

N> r . 

v.; 

(. i ; ' 


? f', 9 

2 8 

2 70 


27 1 

lb 

2 72 


2 77 

8 

2 7 A 


278 

79 

7 76' 

82 

77 7 


' 278 

C 

7 7° 

c 

'» 0 < ' 
£_ - 1 

28 i 

2 p 7 

t; 

A 4 

c 

7 P 8 

81 

2F6 

80 

2 67 
2 F 8 

?v’ : 

a 3 

2 ° : - 

c 

261 

r- 

U 

'? 

c 

V>3 

n 

2 96 

c 

7 9 9 
206 

r 


0 r: R A ( J ) - l? r K R A C J ) + O f R Is ( J ) 

CUNT 1MUC 

UKI T f (o;,P r ') JIZ , J IY f J I X 
FORMA TT / 3 I *> ) 

1 r t M M . L 0 . n M M ) GO T 0 IB 
W is I Tt (< ,?4>NRH 

FORMAT (//26H NUMB PR OF DATA POINTS =tI3/) 

WRI Tt -to ,26 ) NP 

FORMAT ( /2°H NUMBER OF TABULAR POINTS = » 1 3/ / ) 

W» IT E ( 6 , 2 3 ) 

FORMAT ( //15X,1HZ,14X , 1 HT , 1 A X , 3HN02 , 12X , 4HCN02 , i*X , AHDFRK / / ) 

COMPUTE TEMPERATURE FROM A FORM OF THE HYDROSTATIC LAW 

Ni = Nl-l 
DU 6 J = l ♦ N I 

r=o. 

DO 7 I=J , N1 ' 

A,?= DEN ( I + 1 ) / t RHTfl ( I + 1 >**2 ) 

A 1 = 0 EN( n/(KHTO( I >**2) 

A1=DABS CAL) 

A 7= DAB S ( A 2) 

T = T+-CA?-A1)/OLOG(A2/A1) j MZKMU + 1)-ZKM( 1) ) 

T=(3.80 38F-'t)*T/DEN( J) 

T - T’-RHT 0 ( I ).#RHTO( I ) 

T F { MM.NE .MMM > WRIT E (6,22 JZKMIJ ), T,DENt J ) ,ONUMu< J)iDtRR(J) 
FORMAT ( 10X , 5DI5.6) 

CONTINUE ' 

Nil- NH*2 

I F ( MM. EQ. MM.M-) NH -2 
I F ( NN • t'Q *1)00 TO 0 
I F ( mm ,£0* MMM IGO TO 18 
CONTINUE • 

IF ( MN .F'0.2* AMD. MM .NE .MMM )G0 TO 83 
CON T INUE 

I F- ( MM . tQ .MMM ) GO TO 79 
CONTINUE 

I F ( MO .NE . MMM I GO TO 80 
WR I T E ( 6 , 8 2 ) 

FORMAT (//OX , 1HZ , I OX , iOHMEAN E RROR , 7X , 5 HS I GMA , 10 X, 5HSUMSQ, 

1 1 1 X , 3M SUM/ ) * 

COMPUTE STATISTICS F ROM ACCRUED SUMS 

DO 80 J-lfNi 

0MEAN1 J ) =DERRA ( J) / (MM*1 . ) 

S I GMAU ) = DSOR T ( DERRSOt J ) /( MM*1. ) -DM E ANT J)*DMEAN( J ) ) 

WRITE (6 ,81) ZKM( J } , OME AN ( J ) , S I GM A ( J ) , DERRSQ ( J ) , DERR A ( J ) 
FORMAT ( 5D15. 6) 

CO NT I MU E 
' CALL EXIT 
FND 

SUBROUTINE PARAB(M) 

PAR AD FINDS THE SECOND ORDER LEAST SQUARES COEFFICIENTS 
A, 8, AND C FOR INPUT VECTORS Z AND XN » USING M+l SMOOTHING 
POINTS 

IMPLICIT RFAL*8(4-H,0-Z) 

COMMON F SUB ( 400) 
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Z(«00) ,XN(900),M4 00),B(400) , C ( 400 ) , XP ( A ) , FNC t 400 ) , 
lFUM<r00».TKW<400> t TN<mOO),HO,XNO,RE,NP . 

L*M/2 

on 2 i=l»NP 
1 1 = NP + 1-L 

K = 1 1 ' 

Il : ( l . GT » L • AND* J • LT • 1 Z ) K = I -L 

1F(I.GE,U)K = NP-M 

ZR-O. I 

77 8 = 0. 

Z3R=0.‘- 
Z4r3 = 0. 

XNR-O. 

ZNB“0. 

Z?Nii*0. 

K 2-K+M 
DO 3 J=K*K2 

ZB= ZOW ( J J •> 

72h*72tl + 7 ( J ) ** 2 
ZiB=Z3B+7 t J J **3 
Z4R=Z4U+Z{ J ) **4 
/Nrt= ZNB+Z ( J > #XN ( J ) 

Z2NB = 72NN+XN ( J )*( Z( J ) **2 ) 

XNR=XNH+XNt J) 

J=I 
M 1 = M 

Y-( m*Z 3U-Z2I5*ZH )■* (M *Z3lW H*Z 2B )-( M* Z4B- 1 23**2 >*t M* Z 2B-Z *2 ) 

I F ( DABS ( Y ) • LT - 1 .0-4 0 ) A ( J J = 0 . 

lWuSU^-Z2B5zB};(M*23P-ZD*Z2B)-(M.Z4B-Z26**2J*<!1*Z26-ZB**2n 

0= A ( J ) 

YY= M* 12 t3 — / D**2 

I F ( HAB S ( YYJ . LT. 1 - 0-40 ) BtJ)=0. 

I F ( DABS ( Y Y ) *LT • 1 .0- 40 ) GO TO 11 

Bi J)=( (M*ZNB-ZB*XNB)-(’M*Z3B-ZB*Z2B)*0)/CM*Z2B-ZQ*ZB) 

L 

C ( J )= (XNE~0*Z2B-E*ZB )/M 
M = Ml 

CONTI N UC 
RETURN 
e no" 

SUBROUTINE ATMOO(RW) 

ATMOfl FINDS THE 02 NUMBER DENSITY AT ALTITUDE ZW 
IMPLICIT REAL ^8 ( A-H * 0-Z ) 

CrMMfiN Z (400) ,XN( 40 C ) T A { 4 00 > ,8(400} ,C< 4CC) ?XP<4> ,FNf)(400) * 

1 FQF (40U ) , TWW (400 ) t T NO ( 400 ) « HQ, XNG* R E, NP 
PE =6371.05 
Z/! = PW-RE 

XP( 2) =B.72013*DEXP(-(ZW-80.D5 >710.05) 

RETURN 

2.00E-10 .5, COE- 19 1 . 50E-1.8 2.30E-18 8. 00 E- 18 l.iOE-17 

1.90E-1? 1.41F-17 1.42E-17 1.30E-17 

l.OOE-17 8.20E-18 7.00E-18 6.00E-18 4.50E-18 3.bOE-18 

2.50E-18 1.80E-18 1.20F-18 B.00E-19 3.50E-19 3.50E-19 


*. 7 ?.CVr-l*j 1*00 t- I s ? H.00P-20 b.00L-2.0 3. 00020 

3 •' * r l . 300? 0 7* OOP-.? 1 

v o , 0 i\ 0 » 0 i' 0 • 1 ? 0 . 160 * ? 30 .320 • 6 6 C *600 • 7 0 0 • 9 ? 0 • 9 *3 0 

^ ( O’ •‘-001,000 .9 00 .9 no ,970 .960 .8 80 .7 BO .66 0 . BOO .600 

^3*1 .3 20 .260 .200 .170 .110 .080 .060 .060 .030 .020 

no ct rut: 
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Part I: Line numbers 1-62 

In this part, the initial parameter values are set and the 
input data are read. 

Notation 


SH 

VM 

WLO 

TF 

FOO. 


RE 
HO/ HP 
XNO 
NP 


NPP, NPH 
HDP 

NN 

NH 

JIYMAX 

JIZMAX 


0 2 scale height for an isothermal atmosphere, in cm 
visual magnitude of star 

reference wavelength for stellar spectrum, in cm 

Effective blackbody stellar temperature, in °K 

adjusting factor for obtaining a desired fictitious count 

rate at the top of the scan * 

average radius of the earth, in cm 

base height of the O profile, in cm 

£ 

0 2 number density at the base height, in molecules /cm 

number of points in a formulated table of normalized signal 

v. column density. Note: column density = mass of species 
2 

in a 1 cm channel centered on the stellar ray, also called 
the total integrated number density, 
number of data points in the profile 

halfstepsize for computing O densities in an exponential 
atmosphere, in cm 

switch for introducing random error into the signal intensity 
stepsize for computing O 2 number densities, in number of tabular 
data points 

limiting value of the number of data point spacings used, with 
random error in signal 

limiting value of the number of data point spacings used, 
without random error in signal 
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JDC, JIY, JIZ ^counters 


MM, MMM 


SG 


WL 

FILT 


= switch for selecting one of three paths (see program 
comments) 

2 

s' o absorption cross sectional values in cm ; 

the values are listed at the end of the program 
= wavelength in cm 

= normalized values of S4F1 filter transmission; the values 


FIO 


are listed after the SG values at the end of the program 
= stellar spectral intensity, in digital voltmeter counts/ sec 
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Part XI: Line numbers 63-96 


In this part, a normalization factor for the signal intensity is 
computed and a table of normalized signal intensity v. column density (see 
NP under Notation, Part 1) is formed. In making the table, arbitrary 
values of 0 2 number density are selected which cover the expected range of 
values of the species profile. Then the expression for the column density, 

W - / y r is integrated numerically by an expansion 

of the modified Bessel function of the first kind, assuming an isothermal 


atmosphere: 








) 


n Q » number density 

r « distance to the earth's center 
o 

H = scale height 

= modified Bessel function of the first kind 
Superscript * refers to the reference height, which is 80 Km in this case. 
Then the signal intensity corresponding to N is computed from the expression. 



T = filter transmission 

,\ = wavelength 

0" =• species cross section 

I oo * inte g rated Stellar spectral intensity without atmospheric . 
attenuation 

Complete notes on the derivation of the formula for N,are included on the pages 
preceeding Part III. There are some minor differences in notation. 
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Notation 


SUN = integration factor, or the integrated stellar spectral intensity, 

in counts /sec 

HD = stepsize for altitude of the ray's tangent point used in the 
calculation of the column densities, in cm 
IHO » counter and altitude setter; in lines 81-89 it is used to arbi- 
trarily reduce the column density to near zero at the upper 
end of the scan 

SUM * filtered signal intensity integrated over wavelength, in 
counts /sec 

RO = distance from the center of the earth to the ray's tangent 
point, in cm 

ZQ = altitude of the ray's tangent point, in cm 
TWW - column density for the table, computed at altitude ZQ, in cm 
FNO * integrated, normalized signal intensity corresponding to TT 
FSUB = storage vector for the integrated, unnormalized signal 
intensity corresponding to TT, in counts /sec 



Part III: Line numbers 97-174 


In this part, the simulated signal intensity is obtained. If no 
random error is to be admitted, lines 98-108 select certain tabular values 
of the integrated, normalized signal intensity as the simulated signal intensity. 
If random error is to be introduced, there are two possibilities: (1) scatter 
the signal values with a set of random numbers from a Gaussian distribution 
with its mean = 0 and its variance = the signal intensity; (2) find positive 
integral values of the signal from a Poisson distribution having its parameter 
= signal intensity. 

A Poisson distribution is characterized by the probability function, 

P(n) = 0, n not an integer 
T n -I 

p (n) = ■ ■■ n e , , n = 0, 1, 2 , . . . 

I = mean signal intensity, the Poisson parameter in this case 
To find an n which belongs to a Poisson distribution, the pro- 
cedure employed (lines 123, 125, and 145-171) is this: 

1. Seek a random number (YFL) between 0 and 1 which permits a random 
variable, n, to be found such that P(n) < YFL P(n+1). 

a. The search begins with n = 0 and proceeds with increasing integral 
values. This permits the probabilities to be cumulated, i. e. , 

I 3 e “* 

p (3) + p (2) + p (l). 

b. An upper limit on n of I + 5 is imposed for computing economy. 

c. If n ~ 0 is the outcome, substitute a small positive value for 0 (n= 1 0 1 ) 
to allow subsequent. computation to continue. 
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jjc ' 

2. Assign n as the value of the signal 

- V 

In practice, the Poisson random variables were sought when 
0. 1 < I < 10. At the lower limit, P(n I 1) = 0. 12, and the search procedure 
is quite time-consuming at lesser values of I. Below 0. 1, the tabular 
number is accepted as the simulated signal, but a different smoothing pro- 
cedure in the interpolation subroutine produces a number density unlike 
that for the no-scatter case. At the upper limit and beyond, the Gaussian 
distribution yields a real number which is added to the signal in a rapid 

i 

computer operation. 

Notation 

JIZ = loop counter for no -scattering case 
FOF * simulated (measured) signal intensity, in counts/sec 
JIX = loop counter for any quantity of scattering cases with a single 
data-point spacing 

IX 55 GAUSS subroutine argument; it must be changed manually 
at the start of a run if new random numbers are desired; 
it is an odd, positive integer 

JX = RANDU subroutine argument; it must be changed manually 
before each entry to avoid repetition 
AM - GAUSS subroutine argument: the mean of the distribution 

P = probability of a Poisson random variable equalling the integral 
value of the subscript minus 1 

JIY = loop counter for one scattering case with variable data- 
point spacing 

SD = GAUSS subroutine argument: the standard deviation of the 
distribution 

* n-2 in the program corresponds to the n discussed here. 
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V 


N 

NMAX 

U 

D 


GAUSS subroutine argument: the value returned 

indicator of the program path taken and the Poisson integer 

upper limit on the Poisson integer 

numerator in the Poisson formula 

denominator in the Poisson formula 



Part IV: Line numbers 175-238 


In this part/ C> 2 number densities from an exponential 
atmosphere pinned at 80 Km are computed (lines 175-191) for later 
comparison with retrieved number densities. The former values are 

furnished by subroutine ATMOD. 

In order to carry out a mathematical inversion leading to 
retrieved number density, the column density must first be obtained. 

This is achieved through a logarithmic interpolation scheme as follows. 

1. The table (Part III) is converted to the natural logarithms of signal 

intensity (FNO) and column density (TNO). 

2. Using subroutine FARAB, a least squares parabola is fitted to M + 1 

values of ln(FNO). This is done at each tabular value of ln(FNO) to 
obtain the coefficients (A, B, C). 

3. The simulated signal (FOF) is located between tabular values of FNO, 

its logarithm is found, and the column density (N) is computed by 
= exp j^ATn I)^ + B^ In I = CL j . 

Tests with various numbers of smoothing points indicate that 


3(M s 2) is the optimum number for unscattered signal intensities and 
9 (M = 8) is optimum for randomly scattered data. 

The mathematical inversion is represented by 

<-/ 




and the error analysis may be expected to invert N(r Q ) according to an 
expression such as 

. . * 7 *' f 

W,- J) 

However, as may be seen in the ensuing notes, the random noise being 
superposed is 

Jr 

* 
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2 

and smoothing of /IN is necessary, i. e., .dN = A(r-r Q ) + B(r-r o ) + C. 
The resulting formula is 

The coefficients A. and are again supplied by subroutine 
PARAB, the arguments of the table becoming altitude (ZKM) and column 
density (TNO). > 

The error, (n Q - n Q ) / n Q is finally computed at each altitude 
level and its sums and squares of sums are accrued for multiple cases 
(lines 232 -238). 

Notation 



ONUMD, XP(2) - 0 2 number density from the exponential model atmosphere, 

-3 


RHTO, RO 
ZKM 
Z 
XN 


in cm 

distance of the ray tangent point from the earth’s center, in cm 
altitude of the ray tangent point, in cm 
abscissa of the table 
ordinate of the table 


M 

X 

DEN 
DERR 
DERRSQ 
DERRA 
A, B, C, 


one less than the number of smoothing points taken for a 
least squares fit in subroutine PARAB 

entry value for a logarithmic interpolation 

-3 

retrieved number density, in cm 
error in number density 
sum of squares of DERRs 
sum of DERRs 

least squares coefficients for a parabolic fit 
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Part V; Line numbers 239-288 

In this part, the molecular temperature (T) is computed from 
the expression 



which is derived from the hydrostatic approximation (see the following notes). 

If multiple cases are being handled, the mean and standard 

deviation are then computed for the errors at each altitude level. 

> 

Notation 

T * molecular temperature, in °K 
MM, MMM = number of cases for error statistics 
DMEAN * mean error in number density 
SIGMA * one standard deviation of number density error 
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Part VI: Line numbers 353-361 


In this part, the 0 2 scattering cross sectional values are listed 
o , 

at an interval of 20 A (lines 353-358). 

The S4F1 filter transmission values, adjusted to a maximum 
value of 1. 0, are also listed at the same wavelength interval (lines 359-361). 
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USER 
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1 


AVI ; i = 10 r- -SO •KAUHICi CHAVIS' 

SUNOS WAS: 1o:25.!i 0*3-06-7 1 

«»S A 1 1 " SICK I C CN AT 15:32.47 C N 07-12-7 1 

C2NUM 

c a zone numerical erece amuysis 




2 

t 


- 


lEFLlCI.T HiAL*8 (A-ll,C-Z) 

4 


DIMENSION SG (100) , W L (10C) ,210 (100) , FI IT (100) ,CNUMD(400) ,nHTQ(40 

c 


ItEN (400) ,DEhH (400) ,ZM (4CC) ,1N0X (400) ,P (16) , C E R UA (40 0 ) , 

5.1 


21 EH 8 5 6 (4 0 0) ,LKE AN (400) ,5 IGF A (4 00) 

6 


CCMrtCN Z (40 0) , X N (400) , A ( 4C C) ,8 (40 0) ,C (400) ,XP <4) ,FNO (400) , 

7 


IF OF (400) , TW W (4 CO) ,TNC |40 0) ,HC ,XNO, EE, N£ 

7.1 


LATA CE F R A, D EHRSQ/6C C# 0. / 

fc 


Sli = 4 . 34E 5 



X NCC 2 = 8 . E 8 

1C 


2 UO J = 8 3 . D5 

1 1 


SEGO=5. 15 

12 


V ft- 1 . 

13 


W LO = 5000 . C-8 ’ 

14 


IE *1 9000 . 

15 


ECO= . 98 1 C 16 

U 


11=3.141592653589793 

18 


F £ = fc 27 1 . C 5 

15 


H C = 3 0 . D 5 

20 


li P= ilC 

2 1 


> NO = 5 . D 1 C 

2 2 


KE= 17 1 

2 - 


NFP= 141 

24 


MH=NEE 

■ c 
<: - 


H L P= 5 • D4 

26 


NN=2 

2 7 


Ml = 1 

2 7, 1 


J 1Y = C 

2 7.2 


*3IZ=0 

28 

C 


29 

c 

MIEN FK = F M, TilE PROGRAM COMPUTES MEAN ANE STANDARD DEVIATION 

30 

c 

Cl PERCENTAGE ERROR AT HEP SPACING FOB KM CASES. WHEN 

3 1 

c 

FK .NE. EM, THE EROGRAf CCMFUTES PERCENTAGE ERROR AND TEMP. 

j 2 

c 

AT VARIOUS SPACINGS KITH RANDOM ERROR (NN=2) CE WITHOUT 

J J 

c 

RANDCF E E ROE (N*=1) . 

34 

c 


35 


FF=10 

3 6 


ft KM = 10 

3 7 


IE (MK.EC .EMM) WIITE (6 ,20) £ E 

38 

20 

FORMAT (//21H STATISTICS EOR MM=,I3/) 

4 C 

c 


4 1 

c 

INPUT 02 CROSS SEC1ICN (SG) , fcA VE1ENGTH I NTERV AL ( WI) , FITTER 

42 

c 

TRANSMISSION (FILT) 

4 3 

c 


44 


F EAD (5, 6 C C) (SG (I) ,1=1,61) 

45 

600 

F CRM AT ( 6 E 10. 2) 

4C 


DC 300 1=1,61 

47 

300 

KI (I) = 18 C C . + (I- 1) * 2 C . 

4 t 

C 

Si ;V 

49 

C 

EILTER &&A TRANSMISSION 

,5C 

c 


5 1 


READ (5, 601) ( J ILT { 1) ,1 = 1,61) 

5 2 

C < 

601 

f 

FORMAT (11F5.3) 

54 

c 

COMPUTE E1ACKECDY STAR SEICTBUM AND ADJUST TO GET DESIRED 
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C MAXIMUM COUNT BATE 

C 

CC 301 2=1,61 

F10(I)=CEXP (- (.92 1 * V M+ 19 . 367 5) ) * (WLO + *5) / ( (WL (I) * 1 . D- 8) ** 5) * (OBX 
11. 438/ (WLC+TF) ) -1.)/ (CEXP ( 1 • 4 38/ ( WL (I) *1.D-8*TF) ) -1.) +( 1.B-7) 

EIO (I) = F IC (X) *f 00 
301 CONTINUE 

C ... 

C COMPUTE NCBMAIIZATICB F ACT C F 

C 

S C N = C • 

£03020*5,61 ■' 

3 02 S UN -SU N + ( (FILT (0) *EIC (J) ♦ 1.3 IT (J-1) *FIO (J-1) )/5.* (WL <J)-WL (J-1).) ) 
HBITE (6,303) SLN 
303 FCRPAT (/ C 1 5 . 6/) 

EC=17C.C5/ (NP-1) 

C 

C FCBrt A TABLE CF COLUMN DENSITY (TT ) VERSUS SIGNAL INTENSITY 

C (SUM) EASED UPON AN INTEGRATION OF THE FOFFUIA FOB TT UUICH 

C ASSUMES SPHERICAL ST EATI F I C AT IGN AND AN ISOTHERMAL ATMOSPHERE 

C 

' DC 8C0 IHC = 1 ,NP 
• S U M = 0 • 

EC=KE+HC+ (I H C- 1 ) * H C 
Z0=5C-BI 

TT=XNC*EEXP (- (2Q-50. C5)/SE) ♦DSQRT (1. 57079+SH+EO) 

1* (1. *3.*SH/(e.*RO)-15./32«* (SH/RO) *+ 2+3 15 ./2 1 6 . * (SH/RC) **3) 

C 1 =2 ./ (SHCO+SHCC) 

VA=Z£-ZWCC 

VE=VA*1.4142/SHC0 

IF (CAES (VE) ■ G1 . 12.6) TTG^C. 

IF (DAES (*vB) .Gl.12.8) GO TC 14 
CYL1=PCFN 5,VB) 

C YL2=PCF l (- 1 .5 , V E) 

TTG* (XNCC3 + DEXP (- ( VA/SHOO) ** 2) /DSQRT ( 2 . *R C) ) * (.8862269/(01 
1**,75) *DE XP ( VA*VA/ (2 . +3HCC*ShOO) ) +CYL2+EO/ (Cl ** .2 5) * 1 .7724538* 

2 1 EXP |VA*V A/ (2.* SHOO* SHOO) )*CYL1) 

14 TNOX (I HC) =TT+TTG 

IF (IHC.GT.162) TT= 10.*+ ( ( KP-1 BC ) /5. ) 

2 F ( I EC. F£ • 17 1) IT* 1. E- 10 
T WW (IHC) = TT 
IC 10CC J = 2, 61 

1000 J£UH=SUH + ( (Fill (J) *FIC (J)*CEXP (-SG(J) *TWW (I HO) ) + FILT (J- 1) * 

1 F IC (J-1) *CEXP (-SG (J- 1)*TWW (IHC) ) ) /2. ) * { WL (J) -WL (J-1) ) 

f BC (IHC) “SUM/SUN ' , 

WRITE (6,504) IHC, TWW (IHO) ,FBC (IHO) ,TTG 
504 FORMAT (T2,I£,4D15. 6) V 1 

SCO CCNTIBUE 

IF (NN.EC.5) GO TO SCO 

c ' 

C LOOP FOE BO FA EDOM ECHOES IN SIGNAL INTENSITY AND VARIOUS 

C I AT A POINT SPACING* / 

C 

jizmax=£ 

DC 8 J12=1,JI2KAX 
900 DO 8 C 1 1HO* l,NPP,Nb 
SUM=0. 

CC 100 1 J = 2, 61 

1001 3U M = S UM + ( (Fill ( J) * FIG (J) *CFX£ (-SG(J) *TNOX (IHC) ) + FILT (J-1)* 

IF 10 (J-1) *DEXP (-SG (J-1)*TNCX (I HO) ) )/2.) * (WL(J) -WL(J-I) ) 
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IU 

1 13 

1 14 
1 15 

~17 
1 18 
1 IS 
1*10 
1i 1 
122 


1214 
12 £ 
126 
127 
12b 

129 

130 
132 
134 
1 j 6 
136,1 

136 

137 
136 

139 

140 
14 1 



i<m 

14 5 
14c 
147 
1 4 1 
149 
15C 

15 1 

152 

153 

154 

155 
16 * 
157 
15t 
159 
,6C 

16 1 
162 
163 
It 4 

165 

166 
167 


1 7 C 

171 

172 


501 

fcCl 


c 

c 

c 

c 

c 

803 

c 

c 

c 


c 

c 

c 

c 


604 


503 


9 


C 

C 

c 

c 


2 



i CF (110) = £ UM/S UN 

h«IT2 |6,5C1) IRC ,TNOX (I U C ) , 1C F ( IHO) 

FORMAT (T 2,1 6 , 2D 1 5 . 6) 

CONTINUE 

IF (NN. Et. 1) GO 10 ece 
IF (MK. NE .MMM)GC TC 603 

LOOP FCIv ACCRUING E RFC PS It C3 DENSITY RESULTING FROM 
RANDOMLY SCATTERED SIGNAL INTENSITY, FCK MANY CASES, 
WITH A SINGLE DATA F Cl NT StACING 

CC 1b JI >=1 , MK 
CONTINUE 

SET Ut FCli GAUSS SSP AND R AND C SSP 

IX = 2C3* ICC* (JIX-1) 

IF ( M K • N E • K M M ) I X = 2 0 1 

oX^IX 1 

AK=0. 
l (1) =0, 

IF (MK.EC.KMM)GC TC 804 

LOOP FOR RANDOM ERRORS .IN SIGNAL INTENSITY AND VARIOUS 
I AT A PC INI SPACINGS 

J I Y M A X = 1 

CC 2 6 01 W , JIVMAX 
CO 810 1H=1,NPP,NH 
SUM^FCF |IK) *SUN 
IE (SUM. LT .10.) GO TO 9 
5D = DSCRT (SUM) 

CALL GAUSS (IX, SC, AM, V) 

Jr CF (III) = (SUH+V)/SUN 

IE (FCF (IH) . GE. 1.) FOE (IH) =1 .*1 .D-5 

IE (FCF (IH) . LT.FNO (1 ) )FOF (IR) = FNO(1) ♦ 1. D - 1 2 

N= IOC 

GC TC 8CS 

IF (SUB. II. • 1) FCF(IH) =FNO (IH) 

IF (SUM. IT.. 1 ) K = 0 
IE (SUM . LI . . 1) GO TO 8C9 

DETERMINE A RANDOM NUMBER FRCK A POISSON DISTRIBUTION 
MIEN .KSIGNAI I NT E SSITY < 10 . 

NKAX=SUK+5. 

I - DEX F (-SUM) 

E U)-SUK*F 

C-SUM = F (2) • 

CALL RANDU (JX,0Y,YFL) 

DC 1 N = 2 , Nil AX 
IF (N.EC.2)GC TC 10 
U = SU M * * (N-1) * F 
D = 1 . 

I I = N-2 
DC 3 J=1 f II 
£ = ( N-J ) * I 
GSUi*J = GSUK+U/D 
l (N) =GSUK 

IF (P (N) .GE. YFI.AND.F (N-1) . 1E.YFL) GO TO 12 
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17 5 
17 k 
175 

17 6 

W 77 

~7t 

179 

1 8 C 
18 1 
187 
183 
‘6 <4 

185 
18 6 
186.1 
18 7 

18 S 

16 8.1 

189 

ISC 

19 1 
197 

195 
19<J 
15 5 

196 

157 

158 
199 

w 0 0 
2C 1 
202 
2C3 

20 6 
205 
2 C 6 
2C7 
2C 6 
209 
2 1C 
211 
2 12 
213 
2 Hi 
2 15 
2 16 
217 
. 16 
219 
2 2 C 
221 


223 

226 

225 

**27 



809 
6C5 
8 10 
6 C6 


C 

C 

C 

c 

c 



612 

613 


630 

C 

c 

c 


7 CO 


CONTINUE 
8X = J Y 
CC TC 2 
JX=JY 

ECE (IK) = (N-2)/£UN 

IF (FCF (IH > • LI. 1 . D-10 ) EOF (IK) = 1. D- 10 
IF (MM . N 8 .MMM) WRITE ( 6, 8C5) IH, FOF (IH) ,FNO (IH) , N , YFL 
FORMAT (I5,2D15.6 ,15 ,F 10,5) 

CONTINUE 
li EF = iJLF*2. 

IE (MM. HCP 35 . 5C 5 

NFH = <KFH-1)/2+1 
IE (JlY.Et* 1) NPH=161 
IF (JI2.EC-1) NtH = 141 
1 E (MK.E.C .KfiH) NPH=26 1 

COMPUTE C 2 NUMEER DENSITY EACH AN EXPONENTIAL ECDEL 
ATMOSPHERE WITH A SUPERPOSED GAUSSIAN CURVE ECR LATER 
CCMAFJSCN WITH RETRIEVED VALUES 1 

EC 8C7 IHC=1,NFH 
RC=RE+HC+ |IHC-1) *B£E 
CALL ATKCE (SO) 

CKUME (IHC)=XP(2) 

F H1C (I HC ) =30 

ZKM (1HC) - RHTQ (I HO) -RE 

CCNIINUE 

I NT ERIC I ATE IK THE TAELE TC FIND COLUMN DENSITY (TNO) FROM 
A SIGNAL INTENSITY (FOE) UEICH IS A SIMULATION CF THE 
MEASURED SIGNAI 

EC 502 C=1 ,NF 
Z (J) =ELCG (FNO (J) ) 

XN(J) =DICG (TWk(J) ) 

CCNIINUE 
K = 8 

IF (NN. EC* 1) M = 2 

CALL EAHAP(H) 

EC 630 I=1,NEE,NH 

IF (FC fc (I) .LT.INO(I)) EOF (L) »f NC (1) *1.D-12 
DC 612 I = 1 , N F 

1 E (FOE (L) .OE.FNO(I) • AN D. EGE (I) .LT .PNG (I* 1 ) ) GC TO 613 
CONTINUE 

J = I 

X = 0 LCG (F CF (I ) ) 

TNO(L) = A (J) ♦ X*X + B(J) ♦X+'C (J) 

TNO (I) = DE XE (TNC (L) ) 

CONTINUE 

INVERT COLUMN DENSITY (TNC) TC GET 02 NUMBER DENSITY(DEN) 

EC 7 CO J= 1,NPH 

2 (J)=ZKM (J) 

F=J*KH- (NH- 1) 

X N (J)=TKC (K) 

K = 8 

IF (NN.EC *1) M-2 
CALL FABAB(M) 

M = NFH-1 
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2 J 1 


*. j d. 

' 

* 23 3 


2 J 4 


:js 


W 2 J c 

5 

2 3 7 


2 38 


239 


240 

c 

24 1 

c 

142 

c 

243 
243.1 

244 
246 

4 

2 4 c 


247 

85 

2 4 6 


249 


2 6 C 

24 

25 1 


252 

25 

25 3 


254 

23 

255 


25c 


257 


2 5 8 

2 2 

59 

6 

'ftc 

29 

26 1 


2 6 2 


263 


2t 4 

26 

266 


2 1 fc 

18 

267 


2o8 

8 

266 


270 

79 

.• 27 1 

82 

272 


273 

C 

27 4 

C 

27 5 

c 

27 6 


;7 7 


278 


276 


2 8 C 

81 

28 1 

80 

2 6 2 

8 3 

26 3 


284 


r. 

C 


C 

267 

c. 

28 6 

c 

26 6 

c 


C C 4 J S 1 , Ml 
DEN (J) ~C . . 

I C 5 1 = J , N 1 
11 = 2. *A<I) 

IT = EAES(MiTO(l) tiTG (J)**2) 

DHU (J) =11 £ A ( J) ♦ D*(ESQEr(fihTC(i«'1) **2-RllTQ (J) * + 2) - DSQRT (DT) ) 
\ * (U (I) -2 . *A (1) + RE) *D10G ( (1UTC ( I > 1 ) «-DS CRT (RfJTO (I* 1) +*2 
2-KUTO (J) **2) ) / (RHTC (I) +DSQET (DT) ) ) 

DEN (J) --DEU <J)/P1*2, 


COMPUTE E F R 0 B AND ACCRUE IffCf SUMS 


D£iiR (JJ = iUEN (J) -ONUKE (J) ) / C RU K D (J) 

E E I<8 S Q (J) = DERRSQ(J) -* EERR ( J } ♦DERR ( J) 
CEREA(J)=D£RRA(J)+EEFR(J) 

CONTINUE 

hbXT£ (6,85) JI2,JIY,JIX 
FORMAT {/ 3l5) 

IF (Mft.Et .PiIKJGC TC 18 
w RITE (f; , 2 4) NPE 

FORMAT (//26li NUMEEfi OF £ AT A POINTS =,I3/) 

WRITE (t, 2 5) NP 

FORMAT ( / 2 9 ii NUMEER OF TABULAR POINTS =,I3/y) 

* R I T E (6,23) 

FORMAT (// 18X, IfiZ, 1 3X , 3U61C J , 1 2X , 4HC NO 3 , 1 2X ,4fiDE SR//) 

N 1 - N 1 - 1 

EG 6 J = 1 t N 1 

IF (MM. NE . MNM) V, E IT E ( 6 , 22) ZK ft ( J) ,DEN (J) ,CNUMD ( J) -DERR ( J) 

FORMAT (1CX,4D15. 6) 

CONTINUE 
Nri = Nil * 2 

IF (ME . EQ . EMM ) NH=1 
1 F ( N N . E C . 1 ) G C TO 8 
J i (fiM. EO .MMM) GC TO 18 
CONTI NUE 

IF (NN.EQ.2. ANC.MM.NE.KMM) GO TO 83 

CONTINUE 

IF (MK.EQ.NMM)GC TO 7 9 
CONTINUE 

IF ( M ft • NE - KMM)GC TC 80 
WRITE (6,8 2) 

- 1HZ,10X,10HKEAK ERROR, 7 X , 5HSIGM A , 10X,5HSUMSQ, 

1 1 'X^ 


COMPUTE STATISTICS FROM ACCRUED SUMS 
LC 8C J= 1,N 1 

DMEAN(J)=EERRAlJ)/(Mft*1.) 

S IGMA (J) = ESgiiT (DERRSQ (J) / ( ft ti+ 1 . } -D ME A N ( J) *DKEA N (D>) 

K FOREAT # (5I15 K 6) ^ #DKfAN(J) '~ IGMA (J) ,DEERSQ(J) , DERRA ( J) 
CONTINUE 
CALL E JIT 
END 

SUBROUTINE PAEAB (ft) 


PARA8 FINDS Tb E SECCNE CRCffi I EAST SQUARES COEFFICIENTS 
ICINIS Nr C J0E INfUT V£CICiS Z AMD XN ' M + 1 SMOOTHING 
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2^0 

291 ' 

292 
29 i 



296 

297 
296 
29 9 
3CC 
'.0 1 


3 C 2 
30 3 


30 1 


305 
3CC 
307 
3C6 
309 
3 1 C 
31 1 
312 
312 
314 
215 
316 
217 
318 



321 
2 2 2 
323 
224 

325 
3 26 
327 

326 


329 

330 
32 1 
32 2 


3 


1 


11 


2 



336 
) 2 7 
338 
3 3 S 
340 
34 1 
34 2 
342 
344 



347 

346 

349 


C 

C 

C 

c 


IMPLICIT E E A L * 6 ( A - 11 # G- Z } , 

CC.VU.N 2(000) ,XN(400) ,A(4CC) ,£ (400) ,C(400) , XP (4) ,FNO(400) , 

1 f Cl (4 CO) /INK {4 CO) ,'INC (400) ,HC ,XNO, RE # NP 


L=M/ 2 

DC 2 1 ,NP 


12-NP* 1-i 
h = 1 

If II.G1.L.AHD.I.L1.J2) it -I- 1 
II (1.GB.I2) K=NE-M 
ZE = 0. 

2 2 D = 0 • 

Z 3 B- 0 . 

24B = 0 . 

X N B= 0 * 

2 NB=0. 

Z2NB=C. 

N2=K+E 
EC 3 J=K,K2 
2 E = Z B+ Z (2) 

Z 2E-Z 2 E **Z (J) **2 
Z3B=Z3E + Z (J) * *3 
Z 4 B~Z4 E ■+ Z (J) **4 
ZNB=ZNE+Z (J) *XN (J) 

Z2NB=Z2 tE+XM (J) * (Z (J)**2) 
XNE = XNE*XN (J) 

0=1 


r. = Mf 1 i ^ " 

Y = (M*Z'3I“Z2E*ZE) * (M*Z3B-ZE*Z2B)-{M*Z4E-Z2B**2) * ( M*Z2 B~ZB**2) 
IF(EAES(Y)*LT*1* D-40 ) A (0) =0, 

II (I AES (Y) .IT - 1* D-40) GO TC 1 

A (J)= ( (R*ZND-XNti*ZB) * ( tt* 2 3E -ZB+Z2B) - (M*Z2 NB-XK E*Z2 B) * (M + Z2B-ZB* 
1}/ ( (M*Z3E«Z2B*ZB)* (W*Z3 E-2 E*Z2B) - (M*Z4fi-Z 2E**2) * (M«Z2B-ZB**2) ) 

E = A(J) 

Y Y=H*22E-2B**2 

If (LAES (YY) • LT • KE-4C) B(J)=0. 

IF (D AE 3 ( Y Y) . II - 1 • D-4 0) GC 1C 11.* 

E (J) = ( (M*ZNB-ZE*XNE)-(N*Z3E-ZE*Z2B)*D)/(H*Z2B-ZB*ZB) 

E = B (J) 

C (J) = (XNE-D*Z2E-E*ZB)/M 
fc = M 1 

CONTINUE 

BETUBN 

END 

SUBROUTINE ATHCD(RK) 


ATW0BE3SDS TEE 03 NUrtUEE DENSITY AT ALTITUDE ZW FOB AN 
ISOTHEKK AI ATHCSPEERE WITK A SOPEBPOSED GAUSSIAN CURVE 

IMPLICIT BKAI*8 (A-H,C-Z) 

COMMON Z (400) ,XN (400) ,A (4CC) ,B (400) ,C (400) ,XP(4) ,FNO(400) , 
1FC.F (400) f TViW (400) ,1NC (40 0) , EC ,XNO , RE ,NP 
E i=6 37 1 . C5 
ZW-EK-B E 

XP (2)=5.D10*DEJP (- (Z fc-50. C5)/4. 34 D5) ♦8.E8+DEXP {- ( (ZW-83.D5) 

1/5- 05) **<) •* 

RETUBN 

END 

SUBROUTINE GAF.tA (XX/GX,IEE) 

IMPLICIT EEAL*6 (A-H,0-Z) 
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J f: C 

3 5 1 4 
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35 1 
55 3 
358 
J 5 S 
3(>C 
i fc 1 
50 2 

36 i 
JO 4 
365 
306 
JO i 
3o b 
5 c 5 

570 

571 

37 2 


6 


10 

15 


50 

60 

62 

u 4 
70 
80 


1 10 


120 

130 


E FILE 

w 


It ( X X- 57 • ) 6,6,4 

2 Hi = 2 

GX = 1.D75 • 

b E T UK K 
X -XX 

EiU\ = i.on-o 
2 Ki = 0 
G X- 1 . 0 

1 r ( X. ~ 2 • C ) 50, 5C, 15 

IF (X -2 . C ) 110,110,15 
X^X-1.0 

gx na* x 

G C TO 1 C 

IF (X-1 •) 60,120,1 10 
It (X-HUO 62, i 2 , 8C 
Y=FLt AT |IDI NT (X) ) -X 
IE {I AES (Y) - ERR) 150, 130,(4 
IF (U-Y-E.KR) 130,130 ,70 

It (X-1.) 80, 60, 1 10 » 

G X =G X / X 
X=X+1. 

GC TG 7C 
Y=X- 1. 

G Y=1 ,0+Y* (-0.5771017+Y* £0. 9856540 +Y* (-0.876 42 18 + Y* (0.83 28212+ 
1Y* (-0. 566 4729+Y* (C. 2 5 482 05 + Y* (-0. 0514 9930) ) ) ) ) ) ) 

G X=G X*G Y 
C ETUKN 
I L U = 1 
K ITU K N 
E bD 

FUNCTION ICFN (V, X) 

IMPLICIT tiJAL*e {A-a,C-Z) 

ECRTPI = 1. 772 455 65C SC 551 
EC FT 2=1 .4142 1 356237 309 

f AG 1 = 2. ** (-. 5*V-. 5) +Sv>RTPI*DEXP (- .25*X*X) 

CALL GAMA (*5*V+1.,GAK1,IEf) 

CALL CAHEA £. 5*V + . 5,GAM2,JEB) 

E AC 2=C FtiG (. 5*V + „ 5, . 5, . 5*X*X) / GAMl-S^KT2*X*CFHG £.5*V+ 

1 1 • , 1 . 5 ,.5+X*X)/GAK2 
PCFN = t AC 1*FAC2 
EETUFK 
END 

FUNCTION C F H G ( A , C , X ) 

IMPLICIT EEAL*6 £A-H,0-Z) 

EC«=1. 

T 1:30= 1 . 

E = 1 . 

TEEM=TEBM*X* (A + S- 1.)/{S*(C + S-1,)) 

EUB=EIK+T£BK 
S=S + 1 . 

It 13UK.EC.0. ) G C TC 1 

2 t (UAES (TERM/SOW) .GE..0CCC1) GO TO 1 
CFIIG = £UK 

I ITU IN 
END 
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iL 1ST TtilP 


1 

7.8CD-19 

. 7.SCD-19 

7.00C-19 

6.5CD-19 

5.800-19 

5. 20D- 19 

2 

4.7CE-19 

4. 30D- 19 

3.800-19 

2. 4 CO- 19 

3.200-19 

3.20D-19 

3 

3.30D-19 

3. 6 CD- 19 

4.30E-19 

5. 20 D- 19 

6. 60D-19 

8. SOD- 19 

4 

1. 101-18 

1.-4CD- 16 

1.60C-18 

2.2CD-18 

2.70O-1 8 

3.20D-18 

c 

3.9CD-18 

4.SCD-10 

5 . 20 E - 1 8 

5.80D-18 

6. 60C-18 

7. 40D- 18 

6 

8. 00 L- 18 

9.0GD- 18 

1.00D-17 

1 .010-17 

1 .10D-17 

1 .20D-17 

7 

1.200-17 

1.150-17 

1 . 1GE-17 

1. 080-17 

1. OSD-17 

1.00D- 17 

6 

9. 80 L- 18 

9. COD- 16 

8. 200-18 

7. 6 CD- 18 

6. BCD-18 

6.00D-18 

5 

5.0 CD -T8. 

4. 5CD- 18 

4.00E-18 

3.40D-18 

2. 70 D- 18 

2. 300- 18 

10 

1 . 90 E- 18 

1.5CD-10 

1. 200- 18 

5. CCD- 19 

7.000-19 

5.40D-19 

11 

3. 20C -18 






12 

.020 .021 

.022 .024 

.026 .030 

.035 .040 

.048 .055 

.060 

13 

.070 .080 

. C 9 C .100 

.130 .160 

.200 .260 

. 360 . C 20 

.520 

14 

• b 20 . 740 

.840 .940 

• 99C 1 .0 C ( 

.960 .940 

.860 .640 

.800 

15 

.780 .740 

.680 .640 

.600 .500 

.400 .360 

.300 .240 

.200 

16. 

. 170 . 140 

.120 .100 

*.08C • C60 

.CSC .042 

.040 .037 

.035 

n 

.0 33 • C 30 

.028 .025 

.023 .020 


t 
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Part I: Line numbers 1-61 


In this part, the initial parameter values are set and the input 
data are read. 

Notation 


sri 

XN003 

ZWOO 

SHOO 

VM 

WLO 

TF 

FOO 

RE 

HO, HP 

XNO 

NP 


NPP, NPH 
HDP 


NN 

NH 


Og scale height for an isothermal atmosphere, in cm 
0 3 number density at Gaussian peak, in molecules/ cm 3 
altitude of Gaussian peak, in cm 

half -depth of Gaussian layer, in cm 

! > 

visual magnitude of star 

reference wavelength for stellar spectrum, in cm 

effective blackbody stellar temperature, in °K 

adjusting factor for obtaining a desired fictitious count 

rate at the top of the scan 

mean radius of the earth, in cm 

base height of the O^ profile, in cm 

O s number density at the base height in molecules/cm 

number of points in a formulated table of normalized signal 

intensity v. column density. Note: column density = mass of 
2 

species in a 1 cm channel centered on the stellar ray, also 
called the total, or integrated, number density 
number of data points in the profile 

halfstepsize for computing densities in an exponential 
atmosphere, in cm 

switch for introducing random error into the signal intensity 
stepsize for computing number density, in number of data 
point intervals skipped 
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MM, MMM = test for selecting one of three paths (see program comments) 

2 

SG * 0 3 absorption cross sectional values, in cm ; the values 


' WL 
FILT 

FIO 


are listed at the end of the program 
= wavelength, in cm 

= normalized values of S2F5 filter transmission; the values 
are listed at the end of the program 
s stellar spectral intensity, in digital voltmeter counts /time 


interval 



Part II: Line numbers 61. 1-99 


A table of normalized signal intensity v. column density is 
formed in the same way as for C> 2 (p. 10 of O 2 notes). The reference 
height is now 50 Km. 1 

In addition, a representation of the mesospheric peak is 
obtained by use of the Weber cylinder function. To avoid data roughness 
at the altitudinal limits of the superposed Gaussian curve, the cylinder 
function is forced between 40 Km and 130 Km. 

Notation 

SUN = integration factor, or the integrated stellar spectral intensity; 

in counts /time interval 

HD - stepsize for altitude of the ray's tangent point used in the 

calculation of column density, in cm 

IHO = counter and altitude setter; in lines 90-91 it is used to arbitratily 

reduce the column density to near zero at the upper end of the scan 
RO = distance from the center of the earth to the ray' s tangent point, 

in cm 

&Q = altitude of the ray's tangent point, in cm 

TT, TWW = column density for the table, computed at altitude ZQ, in cm 2 
VB = parameter for Weber cylinder function 

Cl, CYL1, CYL2 = parameters for column density in Gaussian superimposed 

curve 

TTG . = column density in Gaussian contribution, in molecules/cm 

o 

TNOX = total column density, in molecules /cm 

SUM = filtered signal intensity integrated over wavelength, in counts/ 

time interval 

FNO = integrated, normalized signal intensity corresponding to TT 

in the table 
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Part III; Line numbers 100 - 181 

This part is essentially unchanged from the Og program* 

A signal is computed from 

X = f -r(x> x„iaK , 

which has the Gaussian superimposed column densities incorporated in N . 

If the signal is scattered, a Poisson distribution is available at the lower end 
of the scan (0. 1 < I < 10.) If the signal intensity is less than 0. 1, the scattering 

is bypassed. 

Notation 

JIZ, JIZMAX - loop counter and its limiting value for no scattering case 
FOF = Simulated (measured) signal intensity, in counts/time interval 

J1X 3 loop counter for MM scattering cases with a specific data 

point spacing 

IX, JX, AM, P, SD, V, N, NMAX, U, D = GAUSS and RANDU subroutine 
arguments; Poisson parameters as in Og program 
JIY, JIYMAX * loop counter and its limiting value for a scattering case 
with variable data -point spacing 
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Part I V: Line numbers 182 - 245 

The procedure and notation are identical to those used in the 0 0 

a 

error analysis. However, subroutine ATMOD yields the model Og through 
solution of the expression 

™ (0 3 ) = rr? £XP [-«-£*)/„] + ~ 

where 

3 > 

n * number density, in molecules /cm 

n* = number density at a reference altitude, in molecules /cm 3 
h = altitude, in cm 
h* = reference altitude, in cm 
H = scale height of O^, in cm 

o 

n’ = 0 3 number density at Gaussian peak, in molecules /cm 
h’ altitude of Gaussian peak, in cm 
H’ - half -depth of Gaussian layer, in cm 

Numerical values are inserted into ATMOD for the above 
parameters. 
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Part V: Line numbers 246 - 283 


Temperature computations are omitted. 

If multiple cases are being handled, the mean and standard 
deviation are computed for the error at each altitude level. 

Notation 

MM, MMM = number of cases for error statistics 

DMEAN = mean error in number density, in percent 

SIGMA = one standard deviation of number density error, in percent 
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Part VI: Line numbers 1-17 


O 3 cross sectional values are listed row by row at an interval 
o 

of 20 A (lines 1 - 11 ). 

S2F5 filter transmission values, adjusted to a maximum value 
of 1 . 0 , are also listed at the same wavelength interval (lines 12 - 17). 



Appendix IV. Orbital Errors 

The effects of the uncertainty in satellite position have been 
traced by imposing small input errors upon the satellite coordinates and 
mapping the error departures of the ray tangent point. The symbols used 
are as in Section 3, Section 3, with the addition of star -superscripts to 
denote, reference, or zero -error quantities. A special case is assumed, 
for convenience, which places the spacecraft and the stellar ray in the 
equatorial plane, and the source star at the First Point of Aries in the 
same plane. 

Tangent ray errors are shown in Figures 1-4 which have 
distance scales on the abscissae to supplement degree indicators of longi- 
tude and latitude. They disclose that altitude errors in tangent ray point 
nearly match vertical satellite displacements; they amount to a few meters 
when a satellite is displaced latitudinally 10 Km; and they reach magni- 
tudes of a few Km when a satellite is displaced longitudinally 10 Km. 

Thus, as expected the tangent ray heights are quite sensitive to the longi- 
tudinal accuracy of spacecraft positioning. 
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Figure 2. Tangent ray height error as a function 
of satellite latitude error. 
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Figure 4 


Satellite tangent ray height error as a 
function of satellite longitude error. 
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Abstract— Stellar ultraviolet light near 1500 A is attenuated in the Earth’s upper atmosphere 
due to strong absorption in the Schumann-Runge continuum of molecular oxygen. The 
intensity of stars in the Schumann-Runge continuum region has been monitored by the 
University of Wisconsin stellar photometers aboard the OAO-2 satellite during occultation 
of the star by the Earth's atmosphere. These data have been used to determine the molecular 
oxygen number density profile at the occultation tangent point. The results of 14 stellar 
occultations obtained in low and middle latitudes are presented giving the night-time vertical 
number density profile of molecular oxygen in the 140-200 km region. In general, the measured 
molecular oxygen number density is about a factor of 2 lower than the number densities 
predicted by the C1RA 1965 model. Also, the number density at a given height appears to 
decrease with decreasing solar activity. Measurements taken at low latitudes during the 
August 1970 geomagnetic storm showed a decrease in the molecular oxygen number density 
at a given height several days after the peak of the storm followed by a slow recovery to 
pre-storm densities. 

r 1. INTRODUCTION 

The first attempt at determining the molecular oxygen concentration in the upper 
atmosphere by u.v. absorption spectroscopy was made on a V-2 rocket experiment in 1949 
(Friedman et al., 1951). The molecular oxygen concentration was obtained from the 
solar u.v. absorption measurements made at various altitudes by a spectrometer aboard the 
rocket. Since then, numerous rocket flights have been made to examine the molecular 
oxygen distribution and its variations, in addition to determining the u.v. spectra of the 
Sun (Byram el a!., 1955; Kupperian et a /., 1959; Jursa et al., 1963; Hall et a/., 1963; 
Hinteregger et al ., 1965; Week and Smith, 1968; Opal and Moos, 1969; Quessette, 1970; 
Brannon and Hoffman, 1971). 

The molecular oxygen distribution in the altitude range 100-200 km has also been deter- 
mined from mass spectroscopic measurements (Schaefer and Nichols, 1964; Nier et a!., 
1964; Hedin and Nier, 1966; Schaefer, 1968; Krankowsky et al 1968; von Zahn and 
Gross, 1969). In addition to molecular oxygen, the other major and minor constituents of 
the upper atmosphere are likewise determined as a function of altitude during the rocket 
flight. 

More recently, satellites have been used to determine the properties of the upper at- 
mosphere from u.v. absorption measurements. The satellites have monitored the attenua- 
tion of solar u.v. radiation in various isolated wavelength intervals during occultation 
at orbital sunrise and sunset. These data have been used to retrieve the neutral air density 
(Thomas et al., 1965; Kreplin, 1965; Venables, 1967; Landin et al . , 1965; Landin et al., 
1967) and the molecular oxygen distribution in the lower thermosphere (Thomas and Nor- 
ton, 1967; Norton and Warnock, 1968; Link, 1969; Stewart and Wildman, 1969; Lockey 
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et a/., 1969; Reid and Withbroe, 1970; May, 1971; Roble and Norton, 1972). These 
measurements, however, are limited only to sunrise and sunset. 

Hays and Roble (1968a) suggested that u.v. stars may be used as source for occultation 
measurements to determine the night-time distribution of molecular oxygen and ozone in 
the lower thermosphere and upper mesosphere. Their calculations showed that in the spec- 
tral region near 1500 A attenuation is primarily due to absorption by molecular oxygen, and 
that Rayleigh scattering and absorption by other minor constituents can be neglected. 

During the past few years, we have used the Orbiting Astronomical Observatory (OAO- 
2) to obtain u.v. stellar occultation data in various spectral intervals. These data have been 
used to obtain the night-time molecular oxygen number density and neutral gas temperature 
in the lower thermosphere from about 140-200 km. The purpose of this paper is to describe 
the experimental technique and discuss the results which were obtained during quiet and 
disturbed geomagnetic conditions. 

2. EXPERIMENTAL TECHNIQUE 

The general details of the stellar occultation technique (Fig. 1) have been described by 
Hays and Roble (1968a, b), Hays et al. (1972) and Roble and Hays (1972). Here we 
describe the specifics of the occultation measurements made by the OAO-2 satellite. 



The OAO-2 satellite has one 16 in. dia u.v. telescope, four 8 in. dia u.v. telescopes, and 
an u.v. spectrometer having a resolution of approximately 5 A. The University of Wis- 
consin optical package, used in conjunction with the telescopes, consists of a series of u.v. 
filters which are used for stellar photometry. Filter (4-1) shown in Fig. 2 has a broadband 
transmission function centered near 1450 A in the Schumann continuum region of molecular 
oxygen. This filter was used to obtain the molecular oxygen distributions from the stellar 
occultation measurements. The detection systems of the u.v. telescopes have a variable time 
integration range. The data shown in Fig. 4, which is typical of the high data rate occulta- 
tion scans, has approximately a 1-5 km altitude resolution between data points. 

Prior to occultation of the star by the Earth, an unattenuated u.v. spectrum of the 
star is obtained from the spectrometer. As an example, the spectrum for a typical bright 
u.v. star is shown in Fig. 3 for the wavelength region of the molecular oxygen u.v. filter. 
Although only a measurement of the relative change of the u.v. stellar intensity is required 
during occultation (Hays and Roble, 1968a), the spectral distribution of the star’s energy 
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Fia. 2. Absorption cross sections for molecular oxygen and ozone in the spectral 
region extending from 1000 TO 3500 A. Dashed curves are the O.A.O. stellar photom- 
eter filter transmission curves used in this study. 



flux is required to determine the transmission through the broadband filter. Thus, a detailed 
stellar spectrum covering the spectral passband of the filter is necessary for data reduction. 

A schematic diagram of an occultation of a star by the Earth is shown in Fig. 1. The 
satellite acquires the star in its telescopes prior to occultation above the absorbing atmos- 
phere. As the satellite moves in its orbit, the source is ultimately occulted by the Earth’s 
atmosphere and the geometry at two satellite positions is illustrated schematically in Fig. 1 . 





HEIGHT, km 
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Fio. 4c. Temperature vs. height deduced from the molecular oxygen number density 

DATA SHOWN IN FlG. 4t> ASSUMING DIFFUSIVE EQUILIBRIUM. 


The intensity of the star is measured as a function of time during occultation by the Earth. 
By knowing the satellite position as a function of time and the star’s position, we can relate 
the star’s intensity to the tangent ray height of the star during occultation. By also knowing 
the position of the star, satellite, and the time, we obtain from geometry the geographic 
position of the tangent ray point. 

The intensity data are normalized to the intensity of the star above the atmosphere and 
the normalized data are related to the tangential column number density of the absorbing 
species, here O a , by the integral relationships expressing Beer’s law 


= (J o ' 


7W.W exp (-eo.WA'o.M HWM dA 


!(£ 


f 


where from the geometry in Fig. 1 , assuming a spherically stratified atmosphere 



n o t ( r ) r dr 
Vr 2 - r 4 * * 


( 1 ) 

( 2 ) 


Here 7(r,) is the photometer count rate at tangent ray height r t and I m is the photometer 
count rate of the stellar signal above the atmosphere. T(X) is the filter transmission and I m (X) 
is the intensity of the star at wavelength X above the atmosphere determined by the spectrom- 
eter; u 0l (2) is the absorption cross-section of molecular oxygen shown in Fig. 2 (Ditchburn 
and Young, 1962; Hudson et al t 1969; Ackerman, 1970) and N 0l {r t ) is the tangential 
column number density of molecular oxygen at tangent ray radius r , ; and n 0j (r) * s the local 
number density at the radius r. 

Molecular oxygen is the sole absorbing species in the wavelength interval determined by 
filter 4-1, therefore the normalized intensity data are converted to tangential column num- 
ber density data using Equation (I). The basic measurements are thus converted to data 
giving the tangential column number density of molecular oxygen as a function of tangent 
ray height. Equation (2) for the tangential column number density is the Abel integral 
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equation (cf. Hays and Roble, 1968a, b) which is inverted to obtain the local number 
density at the tangent ray point assuming spherical symmetry, 


27 r dr Ur r ( V r t 2 — r 2 l 


(3) 


The actual numerical techniques for accomplishing this inversion are described in detail 
in a paper by Roble and Hays (1972). 

The neutral gas temperature is derived by assuming diffusive equilibrium above 130 km 
for the 0 2 number density. Thus, 

To ‘ w= 5wI” s(r>0 - (r ' )dr ' (4) 

where m 0a is the mass of molecular oxygen, k is Boltzmann’s constant, and g is the ac- 
celeration of gravity. 


3. RESULTS 

A relatively large number of absorption profiles were obtained during the time period 
extending from January 1970 to August 1971. Of these, command errors and tracking 
errors have caused serious question about many profiles. The final result of carefully 
selecting only those profiles where altitude is accurately known and where there are no 
command interrupts during the occultation yields 14 useful occultations. Of these, ap- 
proximately 6 are in geomagnetically quiet periods and one series of 8 profiles was obtained 
during the magnetic storm of August 1970. Almost all of these profiles were obtained in low 
latitudes due to the low inclination of the orbital plane of the OAO-2 satellite. In addition 
to the molecular oxygen occultation scans, measurements were also made with a filter located 
in the Hartley continuum of ozone. (Filters 2-5 and 3-2 shown in Fig. 2). A single night- 
time ozone number density profile was presented by Hays et al. (1973) and the results of the 
other scans are presented by Hays and Roble (1973). 


(a) Molecular oxygen profiles during geomagnetically quiet periods 

The geographic position, date, and local time of all of the scans analyzed are given in 
Table 1, Six of these were taken during geomagnetically quiet conditions. The normalized 
intensity measurements obtained on orbit 8884 on 17 August 1970 are shown in Fig. 4a. 
In Figs. 4b and 4c the retrieved molecular oxygen number density and temperature are 
shown and they are compared with CIRA 1965 (model 3, 06: 00 hr.). In general, the re- 
trieved molecular oxygen number density distribution in the 140-180 km region is about 
a factor of two lower than the number densities predicted by the CIRA 1965 model at- 
mosphere at the same local time as the measurements. The composite of all of our quiet 
time data is shown in Fig. 5 along with the results of Krankowsky et al. (1968), Weeks and 
Smith (1968), and Roble and Norton (1972). The envelope encompassing the ozone num- 
ber density profiles obtained from Hays and Roble (1972) and the molecular oxygen number 
density given by the mean CIRA 1965 atmosphere are also shown. The small number of 
stellar occultation profiles makes detailed seasonal variations difficult to discern. How- 
ever, there does appear to be a marked decrease in the molecular oxygen number density 
corresponding to the decrease in solar activity during the period from the beginning of 1970 
to late 1971. The monthly mean of the solar F10-7 flux during this period decreased from 
153 to 120 (x 10 2 ")W cm 2 . The daily solar F10*7 emission for each of the individual 
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Table 1 


Orbit 

Date 

Universal time 
(hr) (min) 

Local time 
(hr) (min) 

Long.* 

(deg) 

Lat. 

(deg) 

Daily solar 
F10-7 emission 
(xlO-* 1 Wcm" ! ) 

5778 


5 

28 

23 

52 

-84 

48 

1729 

8884 


4 

55 

5 

55 

15 

4 

151*6 

8898 

8/18/70 

4 

19 

5 

51 

• 23 

9 

1490 

8913 

8/19/70 

5 

22 

5 

54 

8 

11 

145-4 

8942 

8/21/70 

5 

55 

4 

27 

-22 

24 

142-8 

8943 

8/21/70 

7 

35 

4 

23 

-48 

24 

142-8 

8957 

8/22/70 

6 

58 

4 

22 

-39 

25 

141*7 

8971 

8/23/70 

6 

20 

4 

44 

-24 

27 

146-2 

8985 

8/24/70 

5 

45 

4 

13 

-23 

24 

140*2 

8986 

8/24/70 

7 

25 

4 

13 

-48 

24 

140*2 

9000 

8/25/70 

6 

47 

5 

17 

-32 

27 

J 36-5 

10902 

1/ 4/71 

13 

29 

22 

38 

-223 

-16 

140*2 

11795 

3/ 7/71 

15 

59 

0 

33 

128 

-34 

105-2 

14580 

9/17/71 

6 

22 

1 

50 

-68 

26 

116-5 


* Measured positive eastward from Greenwich. 



Fig. 5. Composite molecular oxygen density results for quiet times. 

.... Krankowsky et a[., 1968. 

scans are given in Table 1. The apparent decrease in the thermospheric molecular oxygen at 
a fixed height probably implies that there is a general decrease in the altitude of the con- 
stant pressure surface at the base of the thermosphere as well as a decrease in the thermo- 
spheric temperature. This suggests a solar cycle influence in the lower thermosphere. 

(b) Molecular oxygen profiles during geomagnetically disturbed periods 
A series of observations was obtained during the magnetically disturbed period at the 
and of August 1970. These data are presented in Fig. 6 where density and temperature at 


4 
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fixed heights are shown as a function of time. For comparison we show the daily sum of 
magnetic K v index on this same figure and give the latitude and longitude of the tangent 
ray point and local time of each occultation scan in Table 1. We note the marked decrease 
in the molecular oxygen number density following the storm at a time when the temperature 
was significantly enhanced. This is obviously in disagreement with the concept of a diffu- 
sively stable atmosphere. The decrease in the molecular oxygen density is apparently due to 



Fig. 6. Molecular oxygen density and temperature variations during the magnetic 

storm op August 1970. 


the existence of a global-scale circulation system which upwells over the polar regions 
and subsides in the equatorial zone. Such a large horizontal scale overturning causes air 
with more atomic than molecular oxygen to move downward over the equator resulting in a 
change in the composition of the lower thermosphere. This result is consistent with the 
observation that energy deposited in the polar thermosphere is redistributed globally 
resulting in a global thermospheric temperature enhancement at low latitudes following a 
large magnetic storm (Hays et a!., 1972b). The daily solar FI 0-7 emission for the days on 
which occultation scans were made is given in Table 1. During the period the daily solar 
F10-7 emission decreased from 151-6 to 136*5 (x IO -22 W cm- 2 ), However, the molecular 
oxygen number density during this period first decreased, but then returned to pre-storm 
levels several days after the storm. This behavior also suggests an atmospheric response 
different from that expected from short term solar e.u.v. changes. 
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4 . SUMMARY 

The stellar occultation technique for measuring molecular oxygen in the Earth’s thermos- 
phere has demonstrated the following general behavior: (a) Thermospheric molecular 
oxygen decreases with decreasing solar activity. This may result from a general cooling of 
the lower thermosphere in addition to the cooling observed in the upper thermosphere. 

(b) Magnetic storms result in strong global circulation systems which upwell in the polar 
regions and sink in equatorial regions. This results in a temporary decrease in the molecular 
constituents in the equatorial thermosphere. 

Acknowledgement — Wc wish to acknowledge the support of Grant NGR 23-005-360, and NASA Contract 
Nas 1 9958. The National Center for Atmospheric Research is sponsored by the National Science 
Foundation. We wish to acknowledge the helpful support of Dr. Houk, and the O.A.O. Wisconsin and 
NASA experimental team. 

REFERENCES 

Brannon, P. J. and Hoffman, J. M. (1971). Molecular oxygen density measurements from 80 to 140 
kilometers. /. geophys. Res. 76, 4630. 

Byram, E. T., Chubb, T. A. and Friedman, H. (1955). Dissociation of oxygen in the upper atmosphere. 
Rhys. Rev. 98, 1594. 

CIRA, Cospar International Reference Atmosphere, (1965). 312 pp. North-Holland, Amsterdam. 
Ditchburn, R. W. and Youno, P. A. (1962). Absorption of molecular oxygen between 1850 and 2500 A 
J. atmos. terr. Phys. 24, 127. 

Friedman, H., Ljchtman, S. W. and Byram, E. T. (1951). Photon counter measurements of solar X-rays 
and extreme ultra-violet light. Phys. Rev. 83, 1025. 

Hall, L. A., Schweizer, W. and Hinteregoer, H. E. (1963). Diurnal variation of the atmosphere around 
190 kilometers derived from solar extreme ultraviolet absorption measurements. /. geophys. Res. 68, 
6413. 

Hays, P. B. and Roble, R. G. (1968a). Stellar spectra and atmospheric composition. /. atmos. Sci. 25, 
1141. 

Hays, P. B. and Roble, R. G. (1968b). Atmospheric properties from the inversion of planetary occultation 
data. Planet. Space Sci. 16, 1197. 

Hays, P. B., Roble, R. G. and Shah, A. N. (1972). Terrestrial atmospheric composition from stellar 
occultations. Science 176, 793. 

Hays, P. B., Jones, R. and Rees, M. H. (1973). Auroral heating and the composition of the neutral 
atmosphere. Planet. Space Sci. in press. 

Hays, P. B. and Roble, R. G. (1973). Observations of mesospheric ozone at low latitudes. Planet. Space 
ScL, 21, 273. 

Hedin, A. E. and Nier, A. O. (1966). A determination of the neutral composition, number density, and 
temperature of the upper atmosphere from 120 to 200 kilometers with rocket-borne mass spectrometers. 
J. geophys. Res. 71, 4121. 

Hinteregoer, H. E., Hall, L. A. and Schmidtke, G. (1965). Solar XUV radiation and neutral particle 
distribution in July 1963 thermosphere. Space Res. V, 1175. 

Hudson, R. D., Carter, V. C. and Breio, E. L. (1969). Photodissociation in the Schumann-Runge Band 
system of O t : Laboratory measurements and atmospheric effects. J. geophys. Res. 74, 4079. 

Jursa, A. S., Nakamura, M. and Tanaka, Y. (1963). Molecular oxygen distribution in the upper 
atmosphere. J. geophys. Res. 68, 6145. 

Krankowsky, D., Kasprzak, W. T, and Nier, A. O. (1968). Mass spectrometric studies of the composition 
of the lower thermosphere during summer 1967. J. geophys. Res. 73, 7291. 

Krepljn, R, W. (1965). NRL solar radiation monitoring satellite: Description of instrumentation and 
preliminary results. Space Res. V, 951. 

Kupperian, J. E., Jr., Byram, E. T. and Friedman, H. (1959). Molecular oxygen densities in the mesos- 
phere at Fort Churchill. J. atmos. terr. Phys. 16, 174. 

Landini, M., Russo, D. and Taoliaferri, G. L. (1965). Atmospheric density in the 120-190 Km region 
derived from the X-ray extinction measured by the U.S. Naval Research Laboratory Satellite 1964-01-D. 
Nature 206, 173. 

Landini, M., Russo, D. and Taoliaferri, G. L. (1967). Atmospheric density measured by the attenuation 
of the solar X-rays monitored on the NRL 1965-16D satellite. Icarus 6, 236. 

Link, R. (1969). Eclipse in Astronomy, 271 pp. Springer, New York. 

Lockey, G. W., Horton, B. H. and Rope, B. (1969). Satellite measurement of upper atmospheric molecular 
oxygen density. Nature, London 223, 387. 



348 


P. B. HAYS and R. G. ROBLE 


May, B. R. (1971). A method of determining the density of thermospheric gases from measurements of 
solar ultra-violet light absorption at grazing-ray and near-vertical incidence. Planet. Space Sci. 19, 27. 

Nier , A. O., Hoffman, J. H., Johnson, C. Y. and Holmes, J. C. (1964), Neutral composition of the 
atmosphere in the 100 to 200 kilometer range. /. geophys. Res. 69, 979. 

Norton, R. B. and Warnock, J. M. (1968). Seasonal variation of molecular oxygen near 100 kilometers. 
J. geophys. Res. 73, 5798. 

Opal, C. B. and Moos, H. W. (1969). Night-time molecular oxygen densities in the 100- to 130-km region 
from Schumann-Runge absorption. J. geophys. Res. 75, 788. 

Qu essette, J. A. (1970). On the measurement of molecular oxygen concentration by absorption spectroscopy. 
J. geophys. Res. 75, 839, 

Reid, R. H. G. and Withbroe, G. L. (1970). The density and vibrational distribution of molecular oxygen 
in the lower thermosphere. Planet. Space Sci. 18, 1255. 

Roble, R. G. and Norton, R. B. (1972). Thermospheric molecular oxygen from solar EUV occultation 
measurements. J. geophys. Res. 77, 3524. 

Roble, R, G. and Hays, P. B. (1972). A technique for recovering the vertical number density profile of 
atmospheric gases from planetary occultation data. Submitted to Planet. Space Sci. 20, 1727. 

Schaefer, E. J, and Nichols, M. H. (1964). Neutral composition obtained from a rocket-borne mass 
spectrometer. Space Res. IV, 205. 

Schaefer, E. J. (1968). Temperature and composition of the lower thermosphere obtained from mass 
spectrometer measurements. Space Res. VIII, 959. 

Stewart, K. H. and Wildman, P. J. L. (1969). Preliminary results of molecular oxygen observations from 
Ariel III satellite. Proc. R. Soc. A3I1, 591. 

Thomas, L. and Norton, R. B. (1967). Absorption of solar X-rays and density changes between 140 
and 160 kilometers. /. geophys. Res. 72, 5552, 

Thomas, L., Venables, F. H. and Williams, K. M. (1965). Measurements of solar X-ray fluxes by the U.S. 
Naval Research Laboratory Satellite 1 965-01 -D. 

Venables, F. H. (1967). Solar X-rays in the wavelength band 44-60A observed by the U.S. NRL satellite 
1965-16D and atmospheric optical densities deduced from these observations. Planet. Space Sci. 15, 681 . 

von Zahn, U. and Gross, J. (1969). Mass spectrometric investigation of the thermosphere at high latitudes. 
/. geophys. Res. 74, 4055. 

Weeks, L. H. and Smith, L. G. (1968). Molecular oxygen concentrations in the upper atmosphere by 
absorption spectroscopy. J. geophys. Res. 73, 4835. 


! 


e? 



APPENDIX VI. 

Observation of Mesospheric Ozone 
at Low Latitudes 


PRECEDING PAGE BLANK NOT FILMED 


113 



Planet. Space Sci. 1973, Vol. 21, pp. 273 lo 279. Pereamon Press. Printed in Northern Ireland 


OBSERVATION OF MESOSPHERIC OZONE 
AT LOW LATITUDES 

P. B. HAYS 

- Departments of Aerospace Engineering, Meteorology and Oceanography, 

University of Michigan, Ann Arbor, Michigan 48104, U.S.A. 

and 

" „ ' R. G. ROBLE 

National Center for Atmospheric Research, Boulder, Colorado 80302, U.S.A, 

(Received in final form 21 August 1972) 

Abstract— Stellar ultraviolet light near 2500 A is attenuated in the Earth’s upper atmosphere 
due to strong absorption in the Hartley continuum of ozone. The intensity of stars in the 
Hartley continuum region has been monitored by the University of Wisconsin stellar photom- 
eters aboard the OAO-2 satellite during occultation of the star by the Earth’s atmosphere. 

These data have been used to determine the ozone number density profile at the occultation 
tangent point. The results of approximately 12 stellar occultations, obtained in low latitudes, 
are presented, giving the nighttime vertical number density profile of ozone in the 60- to 100-km 
region. The nighttime ozone number density has a bulge in its vertical profile with a peak 
of 1 to 2 x 10* cm** at approximately 83 km and a minimum near 75 km. The shape of 
the bulge in the ozone number density profile shows considerable variability with no apparent 
seasonal or solar cycle charge. The ozone profiles obtained during a geomagnetic storm showed 
little variation at low latitudes. 

1. INTRODUCTION 

The first measurement of the ozone concentration in the upper atmosphere was made 
by Johnson et al. (1951) using u.v. absorption spectroscopy. The ozone number density 
distribution was determined from the u.v. absorption measurements made at various 
altitudes by a spectrometer aboard a rocket. Since then daytime measurements of the 
ozone number density distribution have been made from rockets using the Sun as the u.v. 
source (Van Allen and Hopfield, 1952; Johnson et al , 1952; L’vova et al 1964; 
Poloskov et at, 1966; Nagata et al, 1967; Weeks and Smith, 1968; Krueger, 1969) 
and nighttime measurements have been made using the Moon as a source of u.v. light 
(Carver et al, 1967; Carver et al, 1972). The high altitude ozone number density 
distribution has also been determined from solar occultation measurements made from a 
satellite (Rawcliffe.ef al, 1963; Miller and Stewart, 1965). These measurements have 
determined the ozone number density up to an altitude of about 70 km. 

Several other techniques have been used to obtain the ozone number density distribution 
in the upper atmosphere which include (a) satellite eclipse photometry (Venkateswaren 
et al, 1961; Fesenkov, 1967), (b) nighttime airglow spectral photometry (Reed, 1968), 
(c) spectral analysis of backscattered solar radiation as observed from a satellite (Rawcliffe 
and Elliot, 1966; Anderson et al., 1969), (d) chemi-luminescent ozone sondes (Hilsenrath, 
1971) and (e) rocket and ground based observations of the molecular oxygen emission 
Oa^Ag) at 1-27 fi (Evans and Llewellyn, 1970; 1972; Evans et al, 1970). In the latter 
technique, the ozone number density profile is calculated from the measured altitude 
profiles of the Oa^Ag) emission at 1 -27 ft using photochemical theory. Their results give 
the ozone number density profile at twilight up to 100 km and preliminary observations 
indicate a strong seasonal variation of the upper ozone layer at high latitudes. The peak 
ozone number density at 85 km varied between 1-3 X 10 s cm -3 in midwinter to less than 
0*3 X 10 8 cm~ 3 in midsummer. 
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Hays and Roble (1968a) suggested that the nighttime distribution of ozone in the 
upper mesosphere may be obtained from satellite measurements of the intensity of u.v. 
stars during occultation by the Earth’s atmosphere. During the past few years, we have 
used the Orbiting Astronomical Observatory (OAO-2) to obtain u.v. stellar occultation 
data in various spectral regions. The data from the u.v. filters centered at 2390 and 2460 A 
have been used to obtain the nighttime ozone number density distribution from 60 to 
100 km at low latitudes. In this paper we describe the experimental technique and present 
the results which were obtained during quiet and disturbed geomagnetic conditions. 



Fio. 1. Absorption cross sections for molecular oxygen and ozone in the spectral 
REGION EXTENDING FROM 1000 TO 3500 A. 

Dashed curves are the OAO-2 stellar photometer filter transmission curves used in this study, 

2. EXPERIMENTAL TECHNIQUE 

The general details of the stellar occultation technique have been described by Hays 
and Roble (1968a, b; 1972), Hays et al. (1972) and Roble and Hays (1972). Here we 
describe the specifics of the occultation measurements made by the OAO-2 satellite as 
they apply to the determination of the nighttime ozone number density. 

The OAO-2 satellite has one 16 in.-dia u.v. telescope and four 8 in.-dia u.v. telescopes 
and an u.v. spectrometer having a resolution of approximately 5 A. Filter 2-5 and 3-2 
shown in Fig. 1, are the two filters in the University of Wisconsin optical package that 
are used for the ozone stellar occultation measurements. The transmission function of 
these filters is located in the Hartley continuum of ozone ( Fig. 1) with the peak trans- 
mission at 2380 and 2460 A. The detection systems of the u.v. telescope have a maximum 
data acquisition rate of £ sec time integration resulting in a high altitude resolution at 
the tangent point. The intensity data are obtained as a function of time, and by knowing 
the star’s position is orbital elements of the satellite, we can relate the intensity to the 
tangent ray heights during occultation by the Earth’s atmosphere. Figure 2 shows the 
normalized intensity data as a function of tangent ray height for a typical occultation 
scan. Hays and Roble (1972) describe the technique used to relate the normalized intensity 
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data to the tangential column number density of the absorbing species along the ray 
path. This technique allows for the broadband characteristics of the u.v. transmission 
function and can be used as long as absorption is due to a single species. Hays and 
Roble (1968b) have calculated the u.v. transmission of the Earth’s upper atmosphere 
and have shown that the stellar light in the wavelength regions covered by the filter 
transmission functions, shown in Fig. 1, is primarily absorbed by ozone. However, at 



Fig. 2. Normalized intensity profile for a typical ozone occultation. 

altitudes below about 70 km the stellar ultraviolet light is also absorbed by molecular 
oxygen in the Herzberg continuum. The absorption contribution in the Herzberg 
continuum is calculated using the molecular oxygen number density distribution obtained 
from the CIRA 1965 model atmosphere. This absorption contribution is removed from 
the stellar intensity data and the corrected normalized intensity is related to the ozone 
tangential column number density (Hays and Roble, 1972). Once the tangential column 
number density is known as a function of the tangent ray height, the data are inverted, 
using the technique described by Roble and Hays (1972), to obtain the local number 
density at the tangent ray point. 

3. RESULTS 

Stellar occultation measurements using the two filters centered in the Hartley continuum 
of ozone were made during the period extending from January 1970 through August 1971. 
From these data we have obtained 12 orbits on which one or both of the filters could be 
used to determine the ozone density in the mesosphere. These results, corrected for 
molecular oxygen absorption in the Herzberg continuum, and rayleigh scattering are presented 
in Figs. 3-5 where four separate profiles are illustrated in each figure. According to the 
analysis of Roble and Hays (1972), the ozone number density is best retrieved from the 
occultation data between tangent ray heights of 55 to 95 km. The retrieved ozone number 
density data above and below this altitude interval are less reliable. 
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Fio. 3. Altitude profiles for atmospheric ozone obtained using the stellar occultation 

TECHNIQUE. 

Geographic position refers to the tangent ray point of the occultation measurements (time as 
LMT; longitude measured positive eastward from Greenwich). 





Fio. 4. Same caption as in Fio. 3. 
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A careful examination of these ozone profiles in the mesosphere does not show any 
striking systematic seasonal or diurnal pattern in the equatorial regions. There does 
appear to be a systematic increase in the altitude of the high altitude ozone bulge with 
increasing latitude. However, due to the small number of scans, which were made at 
nearly random local time, season, and latitude, it is difficult to place great weight on the 




Fio. 5. Same caption as in Fio. 3. 

slight variations observed. The major conclusion is that ozone varies very little between 
55 and 100 km during the course of the night at low latitudes. The main feature is the 
expected bulge in density which occurs at approximately 85 km and the depletion of 
ozone just below that altitude. This behavior is predicted by most recent theoretical 
studies which incorporate the hydrogen chemistry in their model. 

Numerous theoretical studies of the chemistry of ozone in a moist atmosphere have 
followed the early discussion of Bates and Nicolet (1950) of the influence of hydrogen 
compounds (Hampson, 1964; Hunt, 1966a; Hunt, 1966b; Hesstvedt, 1968; Leovy, 1969; 
Bowman et al., 1970; Hunt, 1971; Shimazaki and Laird, 1972; Strobel, 1972). The 
result of these studies is somewhat confusing due to the large number of choices of possible 
rate coefficients, photodissociation rates, boundary conditions, and eddy mixing rates 
used by these authors. The theoretical results do indicate the general features observed 
in the OAO-2-A2 stellar oceultation measurement. This is illustrated in Fig. 6 where 
three representative theoretical profiles (Hesstvedt, 1968; Shimazaki and Laird, 1970; Hunt, 
1971) are compared with our envelope curve. Hunt (1971) appears to agree best with 
the observations, but it is difficult to assess whether this is fortuitous or the result of a 
correct choice of the multitude of possible variable factors. 
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Fig 6. Comparison of measured ozone envelope with theoretical calculations. 

It should be pointed out that during the period of time in which these observations 
were taken one series of measurements was made while the large magnetic storm of August 
1970 was in progress. There does not appear to be any significant correlation between 
the ozone density in the mesosphere and storm in agreement with the prediction of Maeda 
and Aiken (1968). This is not surprising, but relatively large variation in 0 2 at higher 
altitudes was observed (Hays and Roble, 1972) at the same time. 

SUMMARY 

The stellar occultation measurements of ozone in the mesosphere indicate the following 
conclusions: 

1. Mesospheric ozone varies by as much as a factor of 4 at high altitudes, but does 
not show any clear seasonal or diurnal nighttime pattern. A slight increase in 
the altitude of the 85 km bulge appears to be associated with increasing latitude. 

2. The observations are generally in agreement with results of theoretical predictions 
which utilize a moist atmosphere in which hydrogen compounds are considered 
in the chemistry. 

3. There is no apparent relationship between mesospheric ozone and geomagnetic 
activity at low latitudes. 
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ABSTRACT 


The intensity of stars at wavelengths in the Hartley continuum 
region of ozone has been monitored by the University of Wisconsin 
stellar photometers aboard the 0A0-2 satellite during occupation of 
the star by the earth's atmosphere. These occultation data have been 
used to determine the ozone number density profile at the occultation 
tangent point. The nighttime ozone number density profile has a 
bulge in its vertical profile with a peak of 1 to 3x10^ cm -3 at 
approximately 83 km and a minimum near 75 km. The ozone number 
density at high altitudes varies by as much as a factor of 4, but 
does not show any clear seasonal variation or nighttime variation. 
The retrieved ozone number density profiles define a data envelope 
that is compared with other nighttime observations of the ozone 
number density profile and also the results of theoretical models. 

Calculations are also presented that illustrate the difference 
in retrieving the bulge in the ozone number density profile from 
stellar and solar occultation data. 
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1. Introduction 


Stellar intensity measurements, in certain ultraviolet atmos- 
pheric absorption bands, made by a satellite during occultation of 
the star by the earth* s atmosphere, can be used to obtain the local 
number density of the absorbing species at the occultation tangent 
point (Hays and Roble, 1968a ;b). During the past few years the 
Orbiting Astronomical Observatory (0A0-2) has been used to make 
these occultation intensity measurements using bright ultraviolet 
stars (Hays, et al . , 1972). The occultation technique is based on 
classical ultraviolet absorption spectroscopy. The star is the 
source of ultraviolet light and the stellar photometers aboard the 
QAO-2 satellite are the' detectors . The atmosphere between the star 
and the detectors acts as the absorption cell. During the occulta- 
tion process, the stellar light passes through progressively denser 
regions of the atmosphere and the ultraviolet light is absorbed due 
to the strong absorption features of certain atmospheric gases. The 
basic measurement Is thus the stellar intensity as a function of 
time, but by knowing the star*s position and the satellite orbital 
elements, the intensity is related to the tangent ray height of the 
occulting star. If absorption is due to a single species, the 
normalized stellar intensity is directly related to the tangential 
column number density of the absorbing species through Beer*s Law 
(Hays and Roble , 1972a) . The data giving the tangential column 
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number density as a function of tangent ray height are directly 
inverted, using the technique described by Roble and Hays (1972) 
to obtain the local number density profile of the absorbing species 
at the occultation tangent point. 

Occultation measurements made in the Schumann - Runge continuum 

o 

of molecular oxygen near 1500 A have been used to obtain the 

molecular oxygen number density profile in the lower thermosphere. 

These measurements have been discussed by Hays and Roble (1972a). 

The nighttime number density profile of ozone in the mesosphere has 

also been determined from stellar intensity measurements made in the 

o 

Hartley continuum region of ozone near 2500 A. The results of these 
measurements are discussed by Hays and Roble (1972b). 

In this paper, we compare the results of the QAO-2 ozone 
measurements with theory and other nighttime ozone measurements. We 
also present the results of a single occultation scan and discuss 
the difference in retrieving the ozone number density bulge near 
80 km from solar and stellar occultation data. 
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2. Stellar occultation scan 

The ozone and molecular oxygen absorption cross sections and 
the transmission function of the various ultraviolet filters used 
with the OAQ-2 stellar photometers for the occultation measurements 
are shown in Fig. 1. Filter (4-1) is located in the Schumann - 
Runge continuum region of molecular oxygen and is used to obtain 
the molecular oxygen number density in the lower thermosphere. 
Filters (2-5) and (3-2) are in the Hartley continuum region of 
ozone and are used to obtain the ozone number density profile in 
the mesosphere. Filter (1—4) is used to measure the ozone number 
density in the stratosphere. 

The results presented in this paper were determined from 
occultation data obtained by the 0A0-2 stellar photometers using 
filters (2-5) and (3-2). Hays and Roble (1968b) have shown that 

o 

stellar ultraviolet light near 2500 A is primarily absorbed by 
ozone in the earth’s atmosphere and stellar occultation measure- 
ments within this wavelength region can be used to obtain the 
ozone number density profile in the 60 to 100 km altitude interval. 
In the lower portion of the altitude interval rayleigh scattering 
and molecular oxygen absorption in the Hertzberg continuum have 
small contributions to the total attenuation. These contributions 
are eliminated from the data as discussed by Hays and Roble (1972a) . 

The intensity spectrum of the star above the earth's atmosphere 
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is measured by an ultraviolet spectrometer aboard the 0A0-2 satellite. 
A typical star spectrum, obtained on orbit 12178 on April 3, 1971, is 
shown in Fig. 2. The solar spectrum is also shown in Fig. 2 for 
comparison. The normalized signal of this star, measured during 
occultation using the 0A0-2 stellar photometer with filter (3-2), is 
shown in Fig. 3a. Below 100 km the signal decreases due to absorption 
by ozone. Hie decrease in the signal continues until an altitude of 
about 83 km, then a slight increase in the signal occurs reaching a 
maximum around 75 km. Below 75 km the signal decreases again until 
occultation near 45 km. The ozone number density profile at the 
occultation tangent point is obtained from the intensity data using 
the inversion scheme described by Roble and Hays (1972) and Hays and 
Roble (1972a) and is shown in Fig. 3b. The dip in the normalized 
intensity curve near 80 km, shown in Fig. 3a, is due to a bulge in 
the ozone number density distribution at that altitude. The peak 
ozone number density at the bulge is 3x10^ era* 3 at 83 km. The 
minimum occurs near 73 km, however the magnitude of the number density 
at that altitude is difficult to determine. A complete discussion of 
accuracy of the inversion scheme and an analysis of errors in the 
occultation measurements have been given by Roble and Hays (1972). 
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3 . 


Stellar vs. solar occultation 

The ozone number density profile has also been determined 
from solar occultation measurements made from a satellite 
(Rawcllffe et al . » 1963; Miller and Stewart, 1965). If the intensity 
from the entire solar disk is used during occultation, the effect of 
the finite size of the sun must be considered in analyzing the data. 

In this section we examine whether the bulge in the ozone number 
density near 80 km can be determined from solar occultation measure- 
ments when the entire solar disk is used as the source. For 
simplicity a model of the ozone number density profile is used in 
these calculations. The ozone profile is approximated by a sum of 
gaussian and exponential functions (Roble and Hays, 1972) and the 
profile is shown in Fig. 3b. The normalized intensity calculations for 
a stellar occultation, using GAO-2 filter (3-2), are shown in Fig. 3a. 

The intensity calculations for the broadband filter require a wave- 
length integration over the effective passband of the filter (Hays 
and Roble, 1972a). Therefore, we use the stellar spectrum, shown in 
Fig. 2, the ozone absorption cross-section and the filter transmission 
functions, shown in Fig. 1, to calculate the normalized occultation 
signal for the model ozone number density distribution. The calculated 
normalized intensity profile is shown in Fig. 3a and it has approximately 
the same shape as the 0A0— 2 occultation data. The dip In the intensity 
curve defines the bulge in the ozone number density profile. 
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For a solar occultation, the finite size of the sun must be 
considered because the solar energy emitted from different portions 
of the solar disk is transmitted through the earth*s atmosphere at 
different tangent ray heights. The normalized intensity for a 
solar occultation is calculated using the technique of Roble and 
Norton (1972) and is shown in Fig. 3a. For these calculations^ 
we used the solar energy spectrum, shown in Fig. 2, the ozone 
cross-section and the filter (3-2) transmission function shown in 
Fig. 1. The dip in the intensity curve that occurs in the stellar 
occultation case does not appear in the solar occultation case. 

From normal satellite altitudes the tangent ray heights of the 
upper and lower limb of the solar disk are separated by approxi- 
mately 25 km. Therefore, a smearing of the intensity variation 
caused by the bulge in the ozone number density profile occurs 
when the entire solar disk is used during occultation and a dip 
in the intensity profile is not evident. The calculated normalized 
intensity variation across the solar disk at a tangent ray height 
of 83 km for the center of the solar disk is shown in Fig. 4. 

The tangent ray height of the upper limb is 95 km and the lower 
limb is 70 km. The intensity variation caused by the bulge in 
the ozone number density profile at 83 km occurs near the center 
of the solar disk. The intensity structure across the solar disk 
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is not evident when the overall intensity from the entire solar 
disk is measured during occultation. In theory, the structure 
in the ozone number density profile can be retrieved from solar 
occultation data using the entire solar disk, with an iterative 
technique. However, in practice, the iterative scheme is diffi- 
cult to use especially in the presence of statistical noise. 
Better results would be obtained if a small portion of the solar 
disk is used for the solar occultation measurements . 
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4. Comparison of results 

The ozone number density profiles from 12 stellar occul tat ions 
have been presented by Hays and Roble (1972b) . The measurements 
indicate that mesospheric ozone has a variation of as much as a 
factor of 4 at high altitudes, but does not show any clear seasonal 
or nighttime variation. A slight increase in the altitude of the 
bulge in the ozone number density profile appears to be associated 
with increasing latitude. These measurements plus a few additional 
ones are used to define an envelope representing the ozone number 
density profiles determined from the 0A0-2 stellar occulta tion 
measurements. The envelope is shown in Fig. 5 along with nighttime 
and twilight ozone profiles determined by other experimenters. At 
high altitudes there is good agreement with the ozone number density 
profiles determined by the QAO-2 satellite and the profiles deduced 
by Evans and Llwellyn (1970; 1972). Between 60 and 75 km, the data 
of Evans and Llwellyn (1970; 1972), Evans et al . (1970), Weeks and 
Smith (1968), Miller and Stewart (1965), Tisone (1972) and Rawcliffe 
et al . (1963) are clustered along the upper limb of the QAO-2 ozone 
data envelope. Below 60 km, the 0A0-2 ozone data envelope does not 
agree with the results of Carver, et al . (1967; 1972), Reed (1968) 
and Hilsenrath (1971) * This departure from the other observations 
is probably due to the difficulty of retrieving the ozone number 
density at low signal levels in the presence of statistical noise. 
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Roble and Hays (1972) show that the ozone number density is best 
retrieved between normalized occultation data values of 0.1 - 0.9. 
Therefore, the ozone number density profile determined from the 
GAO- 2 occultation measurements are limited to the 60 - 95 km 
altitude interval. In Fig. 6 the 0A0-2 ozone data envelope is 
also compared with the ozone number density profile determined from 
various theoretical models. In general, there is a large variation 
between the various theoretical models. However, the general features 
of the 0A0-2 ozone measurements appear to be present in the results 
of theoretical predictions which utilize a moist atmosphere in which 
hydrogen compounds are considered in the chemistry. 
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FIGURE LEGENDS 


Fig. 1. Absorption cross-sections for molecular oxygen and ozone. 

The dashed curves give the transmission functions for the 
filters used with the QAO-2 stellar photometers. Filter (4-1) 
is used to obtain the molecular oxygen number density in the 
lower thermosphere and Filters (2-5), (3-2) and (1-4) are used 
to determine the ozone number density in the mesosphere. 

Fig. 2. Solid curve is the stellar intensity as a function of wave- 
length measured by the QAO-2 spectrometer prior to occultation 
on orbit 12178. This star was used to obtain the occultation 
data shown in Fig. 3. The dashed curve is the solar spectrum 
obtained from Ackerman (1970) . 

Fig. 3. a) The dots give the normalized intensity as a function of 

tangent ray height measured by the 0A0-2 stellar photometers 
o 

using the 2460 A filter. The geographic position of the 
occultation tangent point is 35°N, 114°W and time is 22:33 LMT 
on April 3, 1971. The solid and dashed curves are the calculated 
normalized intensities for a stellar and solar occultation 
respectively using the model ozone profile in Fig. 3b. b) The 
dots give the retrieved ozone number density profile from the 
QAO-2 occultation data shown in Fig. 3a. The solid curve ia the 
model ozone number density profile used to calculate the stellar 
and solar intensity profiles in Fig. 3a. 
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Fig. 4. The calculated normalized intensity distribution across 

the solar disk at a tangent ray height of 83 km. The satellite 

altitude at occultation is 715 km. The OAO-2 filter transmission 

o 

function for the 2460 A filter is used in the calculations. The 
tangent ray heights for the upper and lower solar limbs are also 
given in the figure. 

Fig. 5. Comparison of the OAO-2 ozone data envelope with the night- 
time ozone number density profiles determined by various observers. 
Fig. 6. Comparison of the OAO-2 ozone data envelope with the ozone 

number density profile determined from various theoretical models. 
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